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ABSTRACT 



Laboratory measurements are made of the forces and moments 
exerted by water waves on a stationary, submerged, horizontal 
submarine shape and on two-dimensional circular cylinders at 
various relative depths and angles of incidence with respect 
to the wave crest. 

Test conditions include separate cases of drag predomi- 
nance and inertia predominance through variation in the ratio 
u T 
rnsLK 

— g — or its approximate equivalent, the wave height to body 
diameter ratio, H/D. 

Comparisons of the experimental force and moment signatures 
are made with those predicted by two separate theories: 

1. The potential flow solution of Cummins. 

2. A modification of the linear superposition of inertia 
and drag forces proposed by Morrison. In this modification the 
forces and moments are summed in increments of body length using 
drag coefficients determined in steady flow at the local Reynolds 
number and using average inertia coefficients over a full wave 
cycle determined by Keulegan and Carpenter. 

All calculations were performed on an IBM 1620 digital com- 
puter and program statements are included. 

Comparison of theory and experiment shows: 

1. Both theories agree very closely for pre- 
dominance of inertia forces. 
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2. Where inertial forces predominate the predicted 
phase relationship is excellent and the error 
between theory and experiment at maximum force 
has an average value of approximately 19 # for 
the modified Morrison theory and 17# for the 
Potential theory. 

5 . For the case of significant or. predominant 
viscous forces the modified Morrison theory 
gives a very good prediction of phase re- 
lationships and yields an error in maximum 
force having an average value of approxi- 
mately 22# for the modified Morrison theory 
and 32 # for the potential theory. 

4. Instrumentation difficulties are estimated 
to contribute heavily to the apparent error 
particularly with respect to measured wave 
height and vertical moment. 

5. The true error between theory and experiment 
is believed due predominately to the use of 
inexact drag and inertia coefficients. 

6. Additional experiments are needed to define 
the time variation of both the inertia and 
drag coefficients in cases of oscillatory 
flow. 



Adviser: Peter S. E&gieson 

Assistant Professor of 
Hydraulic Engineering 



iii 



ACKNOWLEDGMENT 



The authors wish to express their appreciation to 
Dr, Peter S. Eagleson, Assistant Professor of Hydraulic 
Engineering, their thesis supervisor, for his perceptive 
criticism and valuable guidance during all phases of 
this project. In addition, acknowledgment is due 
Dr. .Robert G. Dean, Assistant Professor of Hydraulic 

i 

Engineering, for his counsel during the preliminary work. 
The authors also wish to thank C. L. Miller, Associate 
Professor of Surveying and the Staff of the Civil Engin- 
eering Computation Laboratory for their willing assistance. 



iv 



TABLE OP CONTENTS 



Page 



I. 



INTRODUCTION 



00000000000000009 



A. Problem Statement .... 



Bo Scope of Investigation 



1 

1 

1 



II . THEORETICAL CONSIDERATIONS ...................... 3 

A. Force Theory ............................. 3 

1. Summary of Previous Work ........... 3 

2. Modified Superposition Theory 4 

3. Potential Flow Theory 17 

Bo Wcive Theory oo.oo.oo»o.....oo......o.o.o.o 19 

C . Computer Program ......................... 21 

1 o General oooo..ooooo.eeeoeoo..o.oo«o. 21 

2. Main Program - Modified Morrison Theory 23 

3. Potential Flow Program ............. 25 

III . EXPERIMENTAL EQUIPMENT AND PROCEDURES .......... 27 

A. Experimental Program ..................... 27 

B. Experimental Equipment 28 

1. Wave Tank and Wave Generator 30 

2 . Test Stand ..... — ................. 30 

3. Instrumentation 30 

a. Horizontal Force Gage 30 

b . Bending Gage ......... — ..... 31 

c . Vertical Force Gage 31 

d. Torsion Gage 32 

e . Wave Gage .................... 33 

f . Recorder 33 

g„ Support Rod Bracket and Shield- 
ing 000000000 . 00 . 00 . 0066 . .OO. 33 



V 



Page 

C . Experimental Procedure ........................ 40 

1. Calibration 40 

a. Force and Moment Transducers ...... 40 

b . Wave Gage 42 

2. Test Procedure 42 

3. Data Analysis 43 

a. Data Records 43 

b. Force and Moment Data Reduction ... 46 

c. Wave Data Reduction 49 

d. Frequency Analysis ................ 55 

e. Error Estimate 58 

IV. PRESENTATION AND DISCUSSION OF RESULTS 62 

A. Presentation of Results 62 

1. Graphical Presentation 62 

2 . Tabular Form ............................ 63 

B. Discussion of Results 63 

1. Model Results - Theory and Experiment ... 63 

2. Prototype Results - Theory Only 66 

V . SUMMARY AND CONCLUSIONS ooo»ooooooooo»oo»o. 0 . 0000 . 9 . 00 . 85 

A . Summary 85 

B. Conclusions 85 

VI . RECOMMENDATIONS FOR FURTHER STUDY 87 

VII . REFERENCES ... 88 

VIII . APPENDICES ... 90 

Appendix A 91 

Appendix B 165 

Appendix C 178 

Appendix D 189 

vi 



LIST OP ILLUSTRATIONS 



Figure 

Number Title Page 

(1) Definition sketch xiii 

(2) Typical incremental body slice 7 

’) 

(5) vs Reynolds number after Schlichting 15 

v u T 

(4) C m vs - after Keulegan and Carpenter 14 

(5) Qualitative variation of particle velocity 

and acceleration in a water wave 22 

( 6 ) Wave tank 54 

(7) Wave generator 54 

( 8 ) Schematic of wave generator 54 

(9) Wave filter 54 

(10 ) Test stand 55 

(11) Experiment models 55 

(12) Torsion gage 55 

(15) Support rod bracket and shielding 55 

(14) Horizontal force gage 56 

(15 ) Horizontal force gage schematic 56 t 

( 16 ) Vertical force gage 57 

(17) Vertical force gage schematic 57 

(18) Torsion gage schematic 58 

(19) Resistance wave gage 59 

(20) Sample calibration curve 4l 

vii 



LIST OF ILLUSTRATIONS 
"[continued j 



Figure 



Number 


Title 


Page 


(21) 


Typical test record 


44 


(22) 


Progressive mean charts 


48 


(25) 


Configuration of gages -torsion gage 


189 


(24) 


Torsion gage circuit diagram 


189 



Table 






Number 


Title 


Page 


I 


Experimental Program 


29 


II 


Magnification Factor and Phase Shift 


58 



viii 



A 

C 

C D 

C m 

D 

E 

F 





LIST OF SYMBOLS 

p 

« cross-sectional area of body, ft 

■ wave celerity, ft/sec 

■ drag coefficient 

■ inertia coefficient 

» diameter of the body (may be a function of 
axial distance), ft. 

p 

» modulus of elasticity, lbs/in 

■ total hydrodynamic force acting on a body 
(vector quantity), lbs. 

» component of total force in the x, y, or z 
(subscripted) direction, lbs. 

= inertial component of F, lbs. 

* drag component of F, lbs . 

« inertial component of force in the subscripted 
direction, lbs. 

» drag component of force in the subscripted 
direction, lbs. 

« maximum value of the forcing function, lbs. 

2 

« modulus of rigidity, lbs/in 

■ wave height, double amplitude, ft. 

4 

« polar moment of inertia, in 

* 2r/L, wave number 

* wavelength of water wave, ft. 

■ mass of fluid displaced by the body, slugs 
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- virtual or added apparent mass, slugs 

« hydrodynamic moment about the y axis of the 
body (vertical moment), ft. -lbs. 

* hydrodynamic moment about the z axis of the 
body (lateral moment ) , ft . -lbs . 

- moment measured by the bending gage, ft. -lbs. 

= wave period, sec. 

■5 

* volume of the body, ft 

* damping coefficient for square low damping 
= equivalent linear coefficient of damping 

- equivalent linear coefficient of damping for 
free vibration. 

= equivalent linear coefficient of damping for 
forced vibration. 

® distance between axis of model and center of 
bending gage, ft. 

■ inside diameter, in. 

* outside diameter, in. 

■ subscript denoting quantity relating to the 
fluid. 

p 

=> gravitational constant, 52.178 ft/sec 
« depth of water from still water to the bottom, ft. 

» spring constant of the system, lbs. /in. 

=* body length, ft. 

* mass of the vibrating system, slugs 

» wave particle true velocity (vector quantity), ft/sec. 

* subscript refers quantity to actual velocity vector. 

■ time, sec. 

* fluid velocity component in the axial (x) 
direction, ft/sec. 



x 



V 


* fluid velocity component in the transverse (y) 
direction, ft ./sec. 


w 


* fluid velocity component in the vertical (z) 
direction, ft ./sec. 


~max 


* refers to the maximum value that quantity has 
attained during the wave cycle. 


X 


* distance along the body axis, ft. 


x forced 


« displacement of the system in forced vibration, ft. 


x free 


» displacement of the system in free vibration, ft. 


y 


** distance in the transverse direction in the 
horizontal plane, perpendicular to the body 
axis , ft . 


z 


« distance perpendicular to the body axis in the 
vertical plane. 


o< 


« angle of attack, measured in the horizontal plane 
between the axis of the body and the direction of 
wave propagation, radians. 


P 


*» angle between vector g and the horizontal plane, 
radians . 


X 


- shear strain at the surface, in/in. 


£ 1 


* principal tensile strain, in/in. 


£ 2 

J* 

J free 


® principal compressive strain, in/in. 
C freeze 


s* 

forced 


* ^forced^c 


A 


» displacement of the free surface from still water, ft. 



» phase shift of the system response, radians. 



M- 


* Poissons ratio. 


^x 


» component of a three dimensional doublet in the 
axial direction, ftvsec. 
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P f 

cr 
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* component of a three dimensional doublet in the 
transverse direction, ft^/sec. 

* component of a three dimensional doublet in the 
vertical direction, ftvsec. 

** 

*» kinematic viscosity, ftvsec: 

- 3.1416 

* density of fluid (water), slugs/ft^. 

- wave angular frequency = 2r/T, sec -1 . 

* shear stress at the outer surface of the tube, 
in/in. 

o 

= velocity potential describing body shape, ft /sec. 

O 

« total velocity potential « 0^ + 0^, ft /sec. 

p 

■ velocity potential of the water wave, ft /sec. 

= wave phase angle, radians. 

= damped natural frequency of the vibrating system, sec -1 . 

* undamped natural frequency of the vibrating system, sec - 
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FIGURE 1. Definition Sketch 
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I . INTRODUCTION 



A . Problem Statement 

The problem of accurate prediction of forces and 
moments on a horizontal submerged body of revolution 
subjected to wave action has recently become quite im- 
portant to several engineering groups . The naval 
architect has become involved with this problem in the 
case of the missile -carrying submarine operating at slow 
speed or hovering near the interface. Another quite un- 
related group, those concerned with the design of off- 
shore oil rigs and such structures on "Texas Towers" 
is also vitally concerned with the wave forces and moments 
on horizontal underwater structural members. 

B. Scope of the Investigation 

The basic objectives of this investigation are as 
follows: 

1. Develop a formulation to predict the force 
and moment signatures of a body of revolution of arbitrary 
shape, under the following general conditions: 

a. Any given wave length, wave height, 
and depth of water (excluding breaking waves). 
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b. Any depth of the object below the 
still water level in which the interaction 
with the free surface is negligible. 

c. Any given angle of attack of the axis 
of the body with respect to direction of wave 
travel . 

d. Any diameter cylinder and any size or 
shape body of revolution submarine. 

Because of the complexities involved in the theoreti- 
cal formulations, high speed digital computer solutions 
are considered a necessity both from the standpoint of 
accuracy and evaluation time. 
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II. THEORETICAL CONSIDERATIONS 



A. Force Theory 

1. Summary of Previous Work 

Considerable experimental and theoretical work has 
been done by other investigators concerning forces and 
moments on vertical pilings and various types of buoyancy 
elements . Outstanding among these studies is the work of 
Harleman (l),(2), Morrison (3)> and Reid and Bretschneider 

(4) . Little information of general utility is available 
for horizontal cylinders however. 

The problem of wave forces on objects has been the 
subject of theoretical studies by Morrison (3), Cummins 

(5) (6), McCamey and Fuchs (7)> and O'Brien (8). The 
first two of these theoretical approaches were adopted 
for application to the case of horizontal bodies: 

(a) The Morrison theory (3) which 
utilizes the concept of a linear super- 
position of inertial and drag components 
to yield total force. Variations of this 
technique are most commonly used to pre- 
dict the forces on vertical piles. 
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(b) A modified potential theory devised 
by Cummins (6) of David Taylor Model Basin 
to predict non-viscous forces and moments 
on a submarine moving beneath a wave train. 

In the opinion of the authors, the Morrison technique 
offered the most accurate approach because of its in- 
clusion (albeit approximate) of viscous effects. Thus the 
Morrison approach, with modifications, was adopted as the 
basic theory for this study. 

The Cummins Potential theory is also applied in a 
parallel development primarily because of its value in 
checking the modified Morrison theory throughout the range 
from negligible viscous effects to cases where viscous 
effects are predominant. 

2. Modified Superposition Theory 

The Morrison theory (J>) operates on the basic assump- 
tion that the total force, F, on a stationary submerged 
object subject to wave action, ffiay be separated into inde- 
pendent inertial and drag components (F^ and F D respectively) 
as described by Eq. (l) 

F - F x + F d (1) 

It should be realized that this separation and assumed 
independence is an expedient and probably a fairly good one, 

but it is not strictly correct. 
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The inertial force is commonly divided explicitly 
into two parts: 

a. The inertia of the accelerating 
fluid particles occuping the body volume 
but calculated in absence of the body. 

b. An additional inertial force due to 
the modification in fluid acceleration 
caused by the shape of the solid body and 
its wake cavity. 

The viscous or drag force as conventionally defined 
also represents the (Implicit) summation of two terms: 

a. The net surface resistance produced 
by the integral of the tangential stresses 
over the body surface. 

b. The form resistance produced by the 
integral of normal stresses over the body 
surface and evaluated under steady state flow 
conditions for the free stream velocity of in- 
terest . 

To the extent that variations in Reynolds number for 
a given acceleration produce appreciable changes in ef- 
fective body shape and volume through variations in the 
boundary layer displacement thickness and the cavity 
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geometry; and to the extent that accelerations for a given 
Reynolds number affect the magnitude of the surface re- 
sistance, the inertial and drag forces as defined will be 
interdependent . 

Any apparent variation in form drag due to the effect 
of acceleration on wake geometry may logically be con- 
sidered an inertial force as all form drag truly is since 
the potential flow pressure distribution for the body plus 
wake cavity will incorporate all resistance save that due 
to tangential stresses. 

Under the above assumption the inertial force may be 
written: 



in which M f is the mass of the displaced fluid, is an 
additional apparent fluid mass, and q is the vector of 
instantaneous fluid velocity. 



F'l » M f dq /dt + dq /dt 



( 2 ) 



Inertia of dis 
placed fluid 



Added inertia 
due to presence 
of body 



Thus 



F-j. = (M + M* ) dq /dt « C m M dq /dt 



( 5 ) 



in which 



C 



m 




( 4 ) 
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and acceleration vectors (and also probably the rate 
of change of acceleration). 



q 




FIGURE 2. Typical Incremental Body Slice 
The total drag force may be written as in equation (5) 




in which the subscript q refers to variability of the 
quantity with the direction (and magnitude in the case of 
Cp ) of the velocity vector , q. 



If the angle of the velocity vector with respect to 
the horizontal plane is defined as 0, then the true z com- 
ponent of the total drag may be written: 



q 




(ID 



If we now assume that: 



q 2 sin 0 = w 2 



( 12 ) 



then 




A * J{D 

q 
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The drag force may be written in terms of the fluid 
velocity as: 

f d *c D p f I |q l q ( 5 ) 

in which the drag coefficient is a function of body shape, 
surface roughness and a body Reynolds number. 

From equations (3) and (5) it can be seen that the 
incremental total force on an incremental slice of a body 
of revolution of diameter D may be described by: 

f 

AF « AF + AFp (6 ) 

AF * AF X + AF y + AF z (7) 



where, for example: 



AF * AF t + AF_ 
Z 1 D 



(8) 



The components of the inertial force may be written: 





C 



m 



/irD 

f Hr 






(9) 



in which: a. x is always oriented along the axis of 

symmetry of the body. 

b. C m is determined by the body geometry 
and its orientation with respect to the total velocity 
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From the above. It can be written: 





P f* p 
2 1 * 





(15) 



The approximation of equation (12) is acknowledged 
to be accurate only for two discrete phase angles of the 
wave motion. It was felt necessary nevertheless, to 
apply this approximation at all phase angles since em- 
pirical drag coefficients for the body shape "seen by 
the vector q" are not available. These coefficients are 
known however for the cylindrical shape presented to the 
component vectors w and v. 

In keeping with this assumption, the C D and C D 

w y 

were obtained from circular cylinder data at Reynolds 
numbers based upon the respective component velocities 
w and y. Similarly the mass coefficients, C m of equation 
(9) will be approximated through use of data from two 
dimensional circular cylinders and will thus be denoted by 
C and C respectively. 

m w "Sr 

In a manner similar to equation (9), the components 
of drag force are now written: 



Pf 

AF « C D -y D |w|w Ax 

z w 

Pf 

AF d » C D -p- D |v | v Ax 

y y 



(14) 
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whereupon the vertical and lateral (perpendicular to 
the body axis in the horizontal plane) components of total 
Incremental force may be written using equations (8), (9), 
and (14): 



AF z 

AF y 




D |w|w + C 
' m. 





D | v| v 




/irD^ n dv 
f 4 dt 



Ax 



(15) 



The vertical and lateral moments (M and M ) about 

z y 

axes through the geometrical center of the body may be ob- 
tained by multiplying equations (15) by the axial position 
measured from the body center: 



AM = AF x 
z z 

(16) 

AMy » AFy x 

Then by summing these incremental forces or inte- 
grating over the length of the body, the total force may 
be obtained. However, these equations cannot be integrated 
directly because every term except p ^ in equations (15) 
and (16) is a function of x, the position along the body 
axis . 

It is further assumed that the average "free stream" 
particle velocity and acceleration to be used in the force 
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equations at a given section is the same as would be 
found at that section on the axis of the body but with 
the body removed from the flow. Note that all three 
of the velocity and acceleration components may be vary- 
ing in magnitude and direction along the axis of the 
body at any particular instant of time. 

The complete equations describing the variations 
of the three components of velocity and acceleration are 
given in Section II-B. 

The coefficients C m and Cp will vary with axial 
position on the body due to the variations in diameter, 
velocity, acceleration, etc., but the quantitative effect 
of these parameters is unknown. For the common case of 
steady flow about a cylinder, experiments have shown a 
well verified relationship between C D and Reynolds 
number (Fig. 2). Likewise for the case of constant ac- 
celeration, theory and experiment have shown the inertia 
coefficient for a circular cylinder to be equal to 2.0. 

A thorough search of available literature revealed no in- 
formation concerning instantaneous values of and C m 
in oscillatory flow. This is due to the fact that it is 
possible to separate the inertial and drag component of 
total force only at four discrete points in a wave cycle - 
points at which the velocity or acceleration go to zero 
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leaving only the inertia or drag components respectively. 

The only information found concerning the variation of C m 

and Cp for circular cylinders in an oscillatory flow was 

that of Keulegan and Carpenter (9). These investigators 

U T 

correlated Cp and C m with an acceleration parameter rP p X ■ 
which can be shown to be proportional to the ratio of 
drag to inertia forces. While both C m and Cp were ob- 
served to vary over the wave cycle, a correlation of only 

U T 

H13.X 

the mean values of these coefficients with — -p — was ob- 
tained. 

Due to the lack of information concerning the in- 
stantaneous values of Cp and C m the following approxi- 
mations were adopted: 

1. Cp is calculated at each point along the body 
from the conventional plot of steady state Cp vs local 
Reynolds number for two dimensional circular cylinders 
(Fig. 3). By this technique, Cp varies over the length 
of the body because both the diameter D, and the velocity 
of the Reynolds number vary with axial position on the 
body . 

2. C and C are determined for respective 

m w m v w T v T 

values of period parameters, -Sp - X - - and ■ - p - 4 x -~, from the 

plot provided in Keulegan and Carpenter (9) and repro- 
duced in Fig. 4. The values of C m in this figure are 
average values of C m over the wave cycle for given values 
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FIGURE 3* Cp vs Reynolds Number* for* Cinculcii* Cylinder’s 
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FIGURE 4. C vs After Keulegan and Carpenter 



of the period parameters . For a particular wave and a 

body of constant diameter along its length, the period 

parameter is constant and the average value of C m is 

determined. For a body of varying diameter along its 

length, such as a submarine shape, the period parameter 

varies along the body and thus the average value of C m 

varies with axial position along the body. As shown in 

Ref, (9), for a given period parameter, the instantaneous 

value of C m does in fact vary. For some period parameters, 

the instantaneous values of C have a large deviation from 

m 

the average value and therefore in these cases the use of 

the average values of C from Fig. 4 can introduce con- 

m 

siderable error. However, it was felt that the average 
values of C m , as shown in Fig. 4, did constitute the best 
information available and that the use of these values 
would provide better accuracy than the use of the constant 
value of 2.0 which is universally used in the case of 
constant acceleration. 

With all the terms of the modified Morrison theory 
defined, the next step is to integrate equations ( 15 ) and 
(16) numerically over the length of the body. A simple 
rectangular numerical integration rule was employed in 
preliminary work. In later work, with an eye towards im- 
proving the accuracy of the integration while at the same 
time reducing the number of numerical integration stations, 
Tchebychef f ' s seven ordinate rule was adopted. The 
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superiority of Tchebycheff ' s seven ordinate rule is 
illustrated by results of a Bureau of Ships comparison 
of Simpson's rule with Tchebycheff ' s rule in which it 
was found that the latter rule was more accurate than 
the twenty station Simpson's rule. Besides being more 
accurate, the Tchebycheff rule is more economical of 
computer time since fewer ordinates are required for a 
given body length. 

By way of summary, the following equations are used 
to represent the forces and moments acting on an elongated 
body of revolution in a train of waves; 



/2 




vertical 








( 17 ) 






F = F 



y lateral 




- 9/2 



r 





vertical 




(18) 



p 



9/2 



M * M. , 
y lateral 
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3° Potential Flow Theory 



The potential flow theory used in this report is an 
adaptation of the theory set forth by Cummins (6). This 
latter theory predicts axial force, lateral force and 
moment, and vertical force and moment on a slender body 
of revolution moving at a constant speed under a sinu- 
soidal wave train. Since this investigation is confined 
to the case of a submarine shape or a cylinder with zero 
forward speed, the potential flow solution was simplified 
considerably. 

As in all potential flow solutions, the assumptions 
must be made that the velocity field is irrotational and 
that a velocity potential 0 ^ does exist . It is further 
assumed that the velocity potential 0 ^ is composed of the 
sum of 0 ^, the velocity potential of the wave without the 
body present, plus 0 ^, the velocity potential of the doublet 
distribution which describes the profile of the body being 
investigated. The diameter of the body is assumed to be 
small compared to the wave length so that the flow which 
would exist in the volume displaced by the body, but with 
the body removed, may be considered approximately uniform. 
The strength of the three dimensional doublets describing 
the body are given as follows: 
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(19) 



^x (x) " - IE d2u 
V y (x) - - g D 2 v 

M z (x) - - H D 2 w 

From these doublet distributions, the forces and 

i 

moments on the body are given as: 



F » -4 ir p. 



+ 9/2 



d^x 

“St 



dx 



F * -4 ir p 



y 



J 



J - i/2 

i/2 






4 - 



J - V 2 



F - -4 IT 



pf 



+ ?/2 

r 



- J ;/2 



d w z 

dt 



dx 



J ?/2 



M =S -4 7r 
8/ 



Pf 



d^y 



x at dx 



^ f /2 



+ y/2 






?/2 



M » -4 ir 



Pf 



^ f'S 



d^z 

x -at dx 



( 20 ) 



( 21 ) 



( 22 ) 



( 25 ) 



( 24 ) 
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The equations are quite a bit simpler than those of 
the modified Morrison theory but, of course, do not in- 
clude any viscous forces. Solution of these equations 
has been carried out using both the Airy (sinusoidal, 
applicable in infinite depth) and Stokian (unsymmetrical , 
applicable in finite depth) wave theories. 

B. Wave Theory 

Early in the investigation, a decision had to be 
made as to which wave theory would best describe the sea 
wave which exerts the force on a hovering submarine or 
an offshore structure. 

It should be pointed out that the sea wave exerting 
the force is not a regular wave but rather a statistical 
array of regular waves. Accepting the notation of super- 
position of harmonics, however, it is of value to study 
the regular wave. For the most practical range of wave 
lengths and water depths this regular wave is described 
by the Stokian theory. 

The following equations as obtained from Ref. (10) 
define the profile, celerity and Internal fluid particle 
velocities and accelerations associated with a Stokian 
wave: 



vV » (K x cos cx. + Ky sinois + ot) 



( 25 ) 
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Wave celerity: 
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Particle velocity: 
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Particle acceleration: 
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The equations describing the sinusoidal wave form 



also used are likewise presented here: 
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A qualitative plot (Figo 5) of the variation of the 
velocity and acceleration for a typical water wave is in- 
cluded at this point as an aid to the reader in understand- 
ing the wave theory and the force theory. 



C . Computer Program 
1 . General 

A few brief comments concerning the computer programming 
problems and techniques should make the program listed in 
the appendix easier to understand. The nature of the equa- 
tions developed to predict the forces and moments on a body 
of revolution invited the use of a high speed digital com- 
puter to perform the complex numerical Integrations. A 



- 21 - 



Direction of Wave Propagation 



«*■ 




Wave Profile 




Vertical Velocity 
(+w } Direction) 




-►-Vertical Acceleration 
(+dw/dt | Direction) 




FIGURE 5, Qualitative Variation of Particle 
Velocity and Acceleration in a 
Water Wave 



well designed computer program also provides added flexi- 
bility or generality to the solution and increases the re- 
liability of the evaluated forces and moments. 
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The IBM 1620 computer of the Civil Engineering Depart- 
ment at M.I,T, was employed exclusively for all the work 
included in this report. It should be noted that this com- 
puter is much slower than the IBM 709 but was used because 
of the ready availability of computer time and because of 
the opportunity to become more familiar with operating 
procedures and programming techniques peculiar to this new, 
popular machine, A listing of the three main progams may 
be found in Appendix C , 

2, Main Program - Modified Morrison Theory 

Because the program was quite long and complex every 
effort was made to compute the various terms in an efficient 
manner. A few of the interesting aspects of the program 
follow; 

a. All eight of the hyperbolic sine and 
hyperbolic cosine terms used in the Stokian 
wave equations are computed by one statement 
(Statement number 1, modified Morrison theory 
program ) . 

b. The vertical and horizontal Stokian 
velocity and acceleration terms (Equations 

28, 29, 30, and 51) are computed in one state- 
ment (Statement number 4). 
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c. The log -log plot of C D vs Reynolds 
number (Fig. 3) was approximated by three 
polynomial equations (each covering a parti- 
cular section of the curve) using a least 
squares technique of curve fitting. The use 
of three separate equations was necessary in 
order to obtain a good approximation of the 
curve with a low order polynomial. The curve 
of Cp vs Reynolds number as derived from 
Schlichting (11 ) and employed on the basis of 
the least square curve fit is shown in Fig. 3« 
Note that this curve was entered for both hori- 
zontal and vertical velocity components to de- 
termine corresponding horizontal and vertical 
C D values . 

d. In a manner similar to that described 

U T 

in o, the plot of C m vs ■ — was approximated 
by three polynomial equations. The original 
curve as derived from Keulegan and Carpenter 
is shown in Fig. 4, 

e. It should be noted that a term called 
"Ratio" was printed out in the results. This 
term is equal to the total vertical drag force 
divided by the total force at that particular 
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time. This was inserted into the program so 
that the magnitude and direction of the verti- 
cal drag force component and the vertical in- 
ertial force component could be calculated if 
desired. 

fo A special test was inserted in the pro- 
gram to check the hyperbolic sine term to insure 
that they were not getting too large for the com- 
puter,, (Largest number allowed is 10 **^). Such 
a condition might arise for a deep water prototype 
wave . 



g. In order to obtain the force and moment 
signature over a wave cycle, all three programs 
were designed to compute these values at eleven 
equally spaced stations corresponding to wave 
phase angles covering 0° to 36O 0 in 36 0 incre- 
ments o The results of each program as presented 
in Appendix A show a listing of these eleven 
groups of values of forces and moments. 

3 . Potential Flow Program 

The potential flow programs have been listed in 
Appendix C and are straightforward and easy to understand. 
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Note that both the sine wave profile potential flow pro- 
gram and the Stokian wave profile potential flow program 
are much shorter and simpler than the modified Morrison 
program. 
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III. EXPERIMENTAL EQUIPMENT AND PROCEDURES 
A, Experimental Program 

An experimental program was selected to test three 
different areas of application of the modified Morrison 
theory : 

1. inertial forces predominant 

2. drag forces predominant 

3. inertial and drag forces about 

the same order of magnitude 

The ratio D/H, was used as a rough measure of the ratio 
of inertial to drag force. Typical values of D/H which cor- 
responded to the three cases above were derived from 
Harleman’s report on vertical cylinders (l). 

In order to obtain more general Information from the 
experimental program, two geometrically similar waves were 
proposed, the first being 10 ft, long and about 8 inches 
high, the second being 5 ft, long and 4 inches high. The 
use of a fairly large wave was deemed necessary in order 
to obtain measurable forces . The use of two wave lengths 
was intended to illustrate the effect of the ratio of body 
length to wave length on the measured wave forces and 
moments . 

Based on the waves selected above and selected values 
of D/H, the following models were constructed for the ex- 
perimental tests (Pig. 11) 
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1" diameter cylinder 

2-7/8" diameter cylinder 

4" diameter cylinder 

5.7” maximum diameter 
submarine shape 



30" long 
30" long 
30" long 
275” long 



aluminum tubing 
aluminum pipe 
plexiglass tubing 
wood 



Since the modified Morrison theory and its associated 
computer program was designed to compute forces and moments 
at any angle of attack and at any depth below the onset of 
surface effects, one model was tested at different depths 
and several models were tested at angles of attack of 45° 
and 90° . 

Table I indicates the model type, orientation, depth, 
etc., for all of the experimental runs performed. 



B. Experimental Equipment 

Fortunately, instrumentation required to measure forces 
arid moments, and laboratory wave making equipment was avail- 
able in the M.I.T. Hydrodynamics Laboratory. This instru- 
mentation was originally designed to measure forces and 
moments on vertical piles and various buoyancy configurations 
proposed for offshore oil rigs. 

A brief description of the available equipment follows. 

A more detailed description of apparatus may be found in 
Ref. (1) and Ref. (2). 
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TABLE I. EXPERIMENTAL PROGRAM 
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ub mar I ne Shape 
-7/8” cylinder 



lo Wave Tank and Wave Generator 



(Figs, 6 through 9) show the 90 foot wave tank with 
the piston type wave-maker at one end. The glass walled 
tank is 37" deep and 30 " wide with a sloping beach wave 
absorber at the far end, 

2. Test Stand (Fig. 10) 

The test stand on which the force and moment gages 
were mounted was constructed of steel pipe and angles 
and was mounted on rubber pads for isolation against ex- 
ternal vibration. It holds the force and moment trans- 
ducers from which the test object is suspended by means 
of two slender, shielded rods passing through the free 
surface. The stand was securely anchored by steel guide 
wires for added rigidity. 

3. Instrumentation 

a. Horizontal Force Gage (Figs. 14 and 15 and 
Ref. l). A highly sensitive Shaevitz linear differential 
transformer is utilized in this "Portal type" gage to yield 
a voltage output proportional to the horizontal force while 
remaining insensitive to any vertical force or bending im- 
posed on the test body. Since this gage was fixed to the 
test stand and could not be rotated, horizontal force in 
the direction of wave propagation only could be measured. 
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b. Bending Gage (Shown mounted directly below the 
horizontal force gage in Figs. 10 and 15, and described 
in Ref. l). A four active arm strain gage bridge mounted 
on a thin walled circular cylinder is the basic con- 
figuration of this gage. The strain gages are so arranged 
that the output of the bridge is four times the output 

of a single gage. While being insensitive to vertical 
force, the gage is sensitive to horizontal force applied 
at a point other than the center of the gage itself. 

Although the modified Morrison theory does not pre- 
dict the axial component of force, it had to be measured 
in the experimental work in order to determine the verti- 
cal or overturning moment. This was necessary because 
the output of the bending gage was the sum of the actual 
overturning moment and that due to the horizontal force 
acting on the body . Overturning moment was found by 
subtracting the axial force times its moment arm from 
the output of the bending gate. Because the gage measuring 
this axial force was fixed in position, overturning moment 
could be measured only for the orientation of the body 
axis perpendicular to a wave crest . 

c. Vertical Force Gage . (Figs. 16 and 17, and 
Ref. l). A Shaevitz linear differential transformer is 
used as the sensing element in this diaphragm-type instru- 
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ment. The gage is constructed so that it is sensitive 
to moment and horizontal force. 

d. Torsion Gage (Pigs. 12 and 18). No suitable 
instrumentation was available to measure the lateral 
or horizontal moment imposed on the body when at an 
angle of attack with the wave motion. A special gage 
was therefore designed and constructed to measure this 
lateral moment. 

The basic requirements or specifications of the 
design were that it be: 

i. sensitive enough to measure the 
torques expected over the range 
of angles of attack set forth in 
the experimental program. 

ii. insensitive to bending, 

iii. Insensitive to vertical force. 

iv. possess adequate stiffness so as 

not to adversely lower the natural 
frequency of the gage assembly. 

v. adaptable with regard to mounting 
to existing gages. 

vi. simple and inexpensive to construct. 

vii. output suitable for the Sanborn re- 
corder. 

Details of the gage are presented in Appendix D and 
Fig. 18. 
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e. Wave Gage , (Figs. 19 a^b^c , and Ref. l). A two 



wire resistance type wave gage was employed for measure- 
ment of the wave profile. Two gages were used, one at 
the location of the body center and another one wavelength 
upstream from this point. 

f. Recorder . (Fig. 6) A Sanborn Model 150 , 

4-channel recording oscillograph was the electronic ampli- 
fying and recording instrument used in conjunction with all 
of the gages employed in this report. 

g. Support Rod Bracket and Shielding . Considerable 
time was spent devising a bracket connection between the 
gages and model support rods that would sustain the forces 
and moments imposed during the .experiments , and at the 
same time permit the body to be easily changed and adjusted 
in position. 

Figure 1J> shows this bracket and the shielding 
(supported independently of the test stand) that was con- 
structed to eliminate the wave forces on the two support 
rods . 

Two important factors in the design of the shielding 
were the alignment of the shielding tubes with respect to 
the support rods and the adjustment of the shielding with 
depth of the model. The use of large diameter shielding 
eases the alignment problem but increases the error in 
measured force and moment due to flow disturbances produced 
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Figure 8. Schematic of wave generator 



Figure 6. Wave Tank 
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Figure 11 Exp. Models 



Figure 12 
Torsion Gage 



Figure 10 Test Stand 
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Figure 13 
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Figure 14. Horizontal Force Gage Figure 15. Horizontal Force Gag© Schematic 
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Figure 16. Vertical Force Gage Figure 17, Vertical Force Gage Schematic 




FIGURE 18 
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a) Photograph 




b) Schematic of Wave Gage 




Figure 19. Resistance Wave Gage 
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by the shielding. With this thought in mind, two sets 
of shielding were constructed (1" O.D. and 1.5" O.D.). 

The shielding size used was the minimum compatible with 
the obtained deflection of the support rods. 

C. Experimental Procedure 
1 . Calibration 

a. Force and Moment Transducers . 

The horizontal force, vertical force, vertical 
moment, and horizontal moment gages were statically cali- 
brated by applying known forces and recording the response 
on the Sanborn recorder. The forces were applied so as 
to obtain the response of the gages in both directions 
from the neutral position. Since there was a large varia- 
tion in the buoyancy of the test models, it was known that 
the zero position of the vertical force gage would be dif- 
ferent with each model. Therefore, this gage was calibrated 
both with no body, i.e., no buoyant force acting, and with 
the very buoyant four inch cylinder attached. This pro- 
cedure shewed that the response of the gage was linear over 
the required range. A sample calibration of the horizontal 
force gage is shown in Fig. 20. 
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Trace Reading, MM 



Figure 20. Typical Calibration Curve 

- 41 - 



£ 



b . Wave Gage . 



The resistance-type wave gages were also 
calibrated "statically" by moving them vertically in equal 
increments above and below the still water line until near 
full scale deflection was achieved on the Sanborn recorder. 
For a discussion of the indirect calibration of the down- 
stream wave gage trace see Section III-C-3-c. 

2. Test Procedure 

Prior to the commencement of the experimental runs, 
a complete calibration of each gage was conducted. Upon 
completion of the calibration procedure, the wave gen- 
erator was started and the controls were adjusted until 
the desired wave length and wave height were achieved. 

The wave generator was then turned off and a few minutes 
later turned on again without adjusting the controls. 
Examination of the wave trace from the Sanborn recorder 
showed that the wave being generated was, within very small 
limits, the same as the originally established wave. The 
procedure was repeated several times with the same re- 
sults; i.e., the wave height, wave length, and frequency 
remained the same once the wave generator controls were 
set . 

The schedule of experiments as given in Appendix A 
was then conducted with each test consisting of the passage 
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of thirty or more waves. In accordance with prior con- 
siderations., the tests involved the use of two different 
waves. The first four runs were made using a wave ten 
feet in length and approximately seven and one half 
inches in height (hereafter referred to as the ten-foot 
wave), and the remaining eleven runs were made using a 
wave five feet in length and approximately four inches 
in height (hereafter referred to as the five-foot wave). 
Upon completion of the test schedule, a complete cali- 
bration procedure was again conducted. 

3° Data Analysis 

a. Data Records 

The experimental data was recorded on 
standard four-channel Sanborn recorder paper. A sample 
of a typical test record is given in Fig. 21. The test 
record shown was obtained from a run in which the 2-7/8 
inch diameter cylinder was positioned with its longi- 
tudinal axis parallel to the direction of wave propagation 
i.e„, with zero degree angle of attack. In this case, 
horizontal force, vertical moment, vertical force, and 
the wave profile were recorded. For those experimental 
runs in which the model was positioned with an angle of 
attack of forty-five degrees, the following items were re- 
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FIGURE 21. SAMPLE TEST RECORD 
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corded; wave profile from the upstream wave gage, lateral 
moment, vertical force, and wave profile from the wave 
gage at the body. These quantities were also measured 
when the model had an angle of attack of ninety degrees 
and in addition lateral force was measured. 

When the body was oriented with an angle of attack 
of forty-five degrees, the horizontal force measured by 
the force gage was not normal or parallel to the body. 
Therefore one could not obtain the overturning moment 
directly from the horizontal force and vertical moment 
gage readings. Because of this and the fact that none of 
the three theories were used to predict horizontal force 
acting at am angle with the body other than zero or ninety 
degrees, neither horizontal force nor vertical moment 
were recorded on these runs. Vertical moment was not re- 
corded for ninety degree angle of attack runs because with 
this orientation there should be no overturning moment ex- 
cept in the case of the submarine shape. 

From the test records one could obtain the force and 
moment signatures. To obtain the values of forces and 
moments at phase angles compatible with those at which 
these quantities had been predicted by theory, the follow- 
ing procedure was followed. 

In accordance with the wave profile equation, the 
phase angle at the crest was defined as ^ = 90° . Since 



the wave length was known, it was possible to locate phase 
angle U/ = 0° at one quarter wave length before the wave 
crest . Having done this, the wave was divided into ten 
equal sections and the forces and moments were read at 
the thirty-six degree intervals for which they had been 
predicted by the computer program. 

It should be noted here that the vertical moment 
recorded on the test records consists of two components. 

One component is due to the horizontal force which acted 
on the model at a known distance from the moment gage. 

The other component is due to the vertical forces which 
acted on the body producing what might be called the over- 
turning moment , It is this overturning moment in which we 
are interested, and thus, the component of vertical moment 
due to the horizontal force had to be subtracted from the 
moment recorded on the test record. The accuracy of this 
subtraction is discussed in Section III-C-3-eo 

b. Force and Moment Data Reduction 

In order to obtain representative force and 
moment values from the test records, it was at first 
thought that an averaging process over approximately ten 
successive waves would be necessary. Therefore, two re- 
presentative test records were investigated to determine 
the difference in values which would be obtained by reading 
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values from one typical wave and from the average of two, 
three, four, and up to ten successive waves. The progressive 
mean charts in Fig. 22 are the result of this investigation. 
The maximum deviation between any one reading and the average 
of the ten successive readings (found in the horizontal force 
reading for the five-foot wave) was 3»97 percent. If one 
carefully selects the wave to be investigated a lesser de- 
viation would be obtained. 

It is felt, therefore, that on the basis of this investi- 
gation one would be justified in obtaining the values of 
forces and moments from one representative wave on each test 
record. The procedure was followed with care being taken 
to select a representative wave rather than an unusual one 
which might differ markedly from succeeding waves. 

The data points were read off the test record in milli- 
meters and were converted to pounds force or foot-pounds 
moment by use of the calibration curves which were obtained 
as discussed in Section III-C-1-a. Since it was known that 
the frequency of the forcing function; i.e., the wave, was 
not greatly different from the natural frequency of the test 
equipment, an investigation was necessarily conducted to de- 
termine whether or not resonance effects were present in the 
recorded values of forces and moments. For a discussion of 
this investigation see Section III -C -3-d. The results of 
the investigation showed that in all cases the phase shift 
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involved was negligible and that there was no significant 
magnification in the records of vertical force and lateral 
moment. The values of horizontal force and vertical moment 
were corrected in accordance with the magnification factors 
determined in the frequency analysis. Having done this, 
the experimental results were tabulated in Appendix A 
and were plotted on the graphs presented in Section IV. 

c - Wave Data Reduction 

The wave data reduction was complicated by the 
fact that the wave gage located at the body did not function- 
properly in recording amplitude during the experimental runs. 
This became evident upon a comparison of the calibration 
traces made for this gage just prior to the commencement of 
the experimental runs and immediately upon conclusion of the 
experiments. These two calibrations differed greatly es- 
pecially in the region below the still water line. As a 
result, one could not feel confident in using either cali- 
bration trace. Furthermore , neither calibration gave wave 
heights consistent with those measured visually during the 
experiments. It is true that it is difficult to measure 
wave heights precisely by marking the levels reached by the 
wave crests and wave troughs at the side of the tanks. How- 
ever, these visual measurements should be accurate within 
one-half inch. The wave heights obtained by use of the 
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calibration curves deviated from the measured values by as 
much as two inches and thus, they were not employed in the 
wave data reduction. 

Some of the test records - in particular, those in 
which horizontal moment was recorded - contain both the 
upstream wave gage trace and the downstream wave gage 
trace. In these cases, it was possible to obtain the 
wave profile from the upstream wave gage trace and cali- 
bration curve. The wave heights obtained from this gage 
did appear to be consistent with the visually measured 
heights, the deviation never beirg greater than one-quarter 
inch . 

Having both wave gage traces on certain records and 
having a reasonable calibration of the upstream gage, it 
was decided to attempt to calibrate the downstream wave 
gage trace on the basis of the upstream trace (assuming 
of course, no change in wave form between the two gages). 

To do this, the trace produced by the upstream gage for a 
given wave was converted from millimeters deflection on 
the test record to feet by means of the calibration curve. 
Thus, at any phase angle the amplitude in feet was known 
from the upstream trace. The deflection in millimeters 
at corresponding phase angles for the same wave on the 
downstream wave gage trace were then known to represent 
a certain number of feet of amplitude. In making this 
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comparison, one had to be certain he was examining the 
traces at the upstream and downstream gages which re- 
presented the same given wave. Since the wave gages 
were located one wave length apart, the corresponding 
wave traces were one length apart on the test record. 

With this in mind, the back calibration was performed 
and the resultant calibration curve was plotted. 

Examination of this curve showed that the portion 
above still water was very nearly the same as that ob- 
tained in the previous direct calibration attempts. The 
portion below the still water level was not so well de- 
fined since the largest trace deflection available for 
back calibration was five millimeters. It was noted that 
the direct calibration curves were nearly parallel from 
about seven millimeters deflection to twenty-five milli- 
meters deflection and that the slope of these lines was 
about the same as the slope at the end of the line which 
was drawn through the back calibration points available. 
Accordingly, the extension of this curve obtained from 
back calibration was drawn parallel to the lines from the 
direct calibration. 

Once this back, or indirect, calibration was com- 
pleted, the amplitudes of the various wave traces were 
checked. The results of this check showed that the peak 
amplitudes obtained by using the back calibration curve 

did correspond closely with those measured visually. 
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One other check was made on this calibration curve. 

Noting that the two direct and the one indirect cali- 
bration attempts resulted in calibration curves which 
showed very little variation above the still water line, 
it was felt reasonable to assume that this portion of 
the calibration was quite accurate. Furthermore, a wave 
profile plotted from the upstream wave gage trace showed 
that the five-foot wave being generated closely approxi- 
mated a Stokian wave. Therefore, assuming first that the 
wave crest obtained from the test record was accurate, 
and second that the wave generated had a Stokian wave 
profile, the wave height (double amplitude) was calcu- 
lated on the basis of Stokian wave theory. The wave 
heights calculated in this manner checked closely with 
those obtained both from the back calibration curve and 
by visual measurement. The following table provides for 
comparison the wave heights obtained by each of these 
methods for several waves. 

Wave Height at Downstream Wave Gage 

By Back Calibration From Stokian Theory 
Wave Visually Using Upstream Gage Using Measured Crest Heigh 

1 7»5 inches 7«4 inches 7-6 inches 

2 7°5 inches 7«5 inches 7»3 inches 

5 4.0 inches 4.0 inches 3.9 inches 
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On the basis of the preceding discussion, it was 
decided to obtain the wave profiles using the back 
calibration curve . These profiles are plotted with the 
results in Section IV and are tabulated in Appendix A. 
Although one could not say that the wave profiles ob- 
tained in this manner are precise, it is felt that they 
are reasonably accurate and that the peak values are 
accurate to within 0.2 inches. 

d. Frequency Analysis 

From the theory of mechanical vibrations, 
it is well known that the response of a mechanical system 
to a periodic forcing function will be subject to in- 
creasing magnification and phase shift as the frequency 
of the forcing function approaches the natural frequency 
of the mechanical system. When working with wave forces 
then, it is desirable to insure that the natural frequency 
of the assembled test equipment is at least seven to ten 
times the frequency of the waves being considered (1) in 
order to avoid the resonance effect. If this condition 
is not met, the resonance effects; i.e., magnification 
factor and phase shift, must be determined for each run. 

The assembled test equipment, including model and 
support rods, was analyzed as a single degree of freedom 



- 53 - 



system sublect to a periodic force and to some unknown 
damping force. Examination of the damped natural fre- 
quency record for the system used revealed that the 
damping force did not satisfy the logarithmic decrement 
condition of linear damping. The damping of bodies os- 
cillating in a fluid is proportional to the hydrodynamic 
drag. At moderate Reynolds number, the steady state drag 
coefficient is constant and the drag is therefore pro- 
portional to the velocity squared. Although it has been 
shown that the steady state drag coefficient does not 
apply exactly to oscillating systems (4), it is reasonable 
to assume that the damping force is proportional to the 
velocity squared. 

For square law damping the equation of motion may 
be written ass 

mx + C.Jc 2 - kx » F sin ui t (37) 

4 & o 

in which the symbols are defined as follows s 

m is the mass of the system; 

Cp is the damping coefficient for 
square law damping; 

k is the spring constant for the system; 

F is the maximum value of forcing function; 
o 

oo is the frequency of the forcing function. 
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Jacobsen (12) has obtained an approximate solution of this 
non-linear differential equation by substituting an equi- 
valent linear coefficient of damping. The equation then 
becomes 



rax + C x + kx = F q sin uj t ( 58 ) 

# 

in which C is the equivalent linear coefficient of 
damping . 

The approximate equivalence of these two systems is 
established by letting the work done per cycle by the 
linear damping force be equal to the work done per cycle 
by the square law damping force. Therefore the following 
equation may be written. 



IT IT 

jsa .0 r^> *. 

AW = 4 [ C 2 x dx = 4 / C x dx (39) 

o o 



Noting that C 2 and C are constant over one cycle and 
assuming the response of the system is periodic and of the 
form: 

X - X Q sin (wt - 9) (40) 

equation (39) may be rewritten as follows: 

T r TT 

cos 2 (wt-O)dt 



C 2 x^ 5 



2oj 



cos^ (wt-O)dt 



* 2 2 
C XTu) 
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Since the values of the integrals do not depend upon 
one may perform the integration and solve for C. 

c *"|r x4y C 2 (41) 

This expression may be used in obtaining the value 
of the ratio of the equivalent linear damping to critical 
damping o If this ratio is known, one may enter a standard 
plot of magnification factor or phase shift versus forced 
to natural frequency for linearly damped systems and there- 
by determine the magnification factor and phase shift. 
Considering both the damped free vibration and damped forced 
vibration, the following expression can be written: 



forced 



free 



ojx 

~cJx 



forced 

free 



(42) 



* v* * 

C ® ^ , o C } C a ® ^ oG 

-forced forced c free free c 



(43) 



Substituting for c forced and C f re e from e S uation W) 
and rewriting equation (42) results in the following: 



T* 

J * ^ x f orced free 

forced “ ^ x free 



(44) 



where is the damped natural frequency of the system. 
In this equation u) and x f orce( j can be read from the test 
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records, and , x f ree J ^ free can be det€rmined 

from the free vibration record. In the case of linear 
damping, the logarithmic decrement can be written as 
follows : 



in ^ 

X B 



2rr J 






free 



)' 



( 45 ) 



where x^ and x^ are the amplitudes of any two successive 
vibrations . 

Therefore, by selecting any two successive vibrations, 
S f ree can te determined for that cycle. If that value of 
i tree and the value of x^ are substituted into equation 
(44), J* forced can be determined. 

The plots available in the literature are made with 
the abscissa equal to , the ratio of the frequency of 
the forcing function to the undamped natural frequency of 



the system. From the test records the ratio 



to_ 

UJa 



is avail- 



able. However, the following expression holds : 






* ) 2 
free 1 



From this it is clear that — is very nearly equal to 77 - 

for small values of For each of the experiments 

* , 

conducted, the value of <J was less than 0.04 and. 



therefore, in these cases — may be considered equal to 

^ 4k 



6 0 
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In accordance with the preceding, the values of 
magnification factor and phase shift were determined 
for each test and were applied to the experimental re- 
sults during the data reduction process, A table of 
magnification factors and phase shifts found is given 
below. 



TABLE II, MAGNIFICATION FACTOR AND PHASE SHIFT 



Body 


Wave 


Magnification 

Factor 


Phase 

Shift 


4 m„ cylinder 


10 


1,03 


0 


4 m, cylinder 


5 


1,06 


0 


2-7/8 cylinder 


10 


lo05 


0 


2-7/8 cylinder 


5 


1,10 


0 


1 in , cylinder 


10 


1,00 


0 


1 in, cylinder 


5 


1,00 


0 


Submarine shape 


10 


1,03 


0 


Submarine shape 


5 


1,05 


0 



e. Error Estimate 

A complete error analysis has not been 
attempted. However, because of the subtraction procedure 
involved in obtaining the value of overturning moment, it 
was considered essential that an estimate of the maximum 
error possible in this quantity be made. As discussed In 
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Section III-C-3-b, the overturning moment was evaluated 
by subtracting the component of vertical moment due to 
horizontal force from the total vertical moment measured 
by the bending gage , Of the two components of vertical 
moment, the component due to horizontal force was by far 
the greater , Therefore, the procedure entailed the sub- 
traction of a factor of nearly the same magnitude as the 
total measured vertical moment from the total measured 
vertical moment. This resulted in a relatively small 
factor for overturning moment. From this it is clear that 
a small percent error in total moment, or in moment due to 
horizontal force, would result in a large error in over- 
turning moment. This effect is most severe when the over- 
turning moment is the least percent of the total vertical 
moment. In particular, the overturning moments determined 
on the experiments using the ten-foot wave would be subject 
to the largest errors. 

An estimate of the maximum error was made using the 
following simplified method: 



jin - 
d(M z ) 



- F , d 

X 

fn(M y - P x .d) 
d(W v - F x .d) 



(45) 
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in which the terms are defined as follows: 



M,. is the overturning moment , 

is the total vertical moment, (measured at gage) 
F is the horizontal force acting on the body, 

«A» 

d is the moment arm; i„e., the distance from 
the center of the test body to the center 
of the bending gage. 

Since we are looking for the maximum error, we are inter- 

d(M -F ,d) 

ested in the maximum value of ■ ^ V TJ - x * — , This will ob- 

vV 

viously be infinite when overturning moment is zero. For 
non-zero overturning moments we will examine the maximum 
value of the numerator d(M v ~F x ,d), To estimate this term, 
the maximum error in reading M m , was assumed to be 1/4 milli- 
meter, which in the case when the attenuation X10 scale was 
used, is equivalent to +0,05 ft, lbs. It was assumed that 
the moment arm could be measured exactly and that the maxi- 
mum error in reading F was 1/4 millimeter, which in the 
case of the attenuation X20 scale is equivalent to + 0,06 

lbs. Since the moment arm was measured to be 55-5 inches, 

$ 

the maximum error in this component of total moment is 
+ O.167 ft. lbs. The maximum error in overturning moment 
then is estimated to be: 

d(^) _ 0.050 + O .167 

M M - F ,d 

z v x 
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Specifically for test run number 1, the estimated error 
in peak value of overturning moment is; 



a( M K> 

H, 



0.217 

OTTO 



48.6 % 



This simplified analysis was made to provide an estimate 
of the error which might obtain in the evaluation of over- 
turning moments by the method used in these experiments . 

The accuracy of the gages and recording instruments has not 
been considered but rather only the accuracies involved in 
reading the recorded data have been considered. However , 
this does show the sensitivity of the resulting overturning 
moment to any inaccuracy in the reading of the Sanborn 
recorder trace. It should be emphasized that this applies 
only to the overturning moment and not to the other quantities 
reported. Furthermore , it should be noted that sensitivity 
results from the method used to obtain overturning moment 
and that this method was dictated by the limitations of 
the instrumentation. 
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IV. PRESENTATION AND DISCUSSION OF RESULTS 



A . Presentation of Results 
1 . Graphical Presentation 

The theoretical and experimental results are pre- 
sented in graphical form in this chapter as an aid to the 
reader in making a qualitative comparison of the results. 

Each graph presents a complete time signature of the re- 
sults obtained from the modified Morrison theory, from the 
Stokian potential theory, and from the experimental runs. 

The following points concerning the graphs should be noted : 

(1) Theoretical axial force is plotted for the 
Stokian potential flow only sire e the modified Morrison 
theory had not been used to predict this quantity. 

( 2 ) The sine wave potential theory results are 
plotted for run number 4 only. 

( 5 ) For the runs with non-zero angles of attack, 
experimental results are not plotted for each of the 
quantities predicted by theory. This is because limitations 
in the instrumentation prevented the experimental measure- 
ment of axial force, horizontal force, and overturning moment 
for the 45° angle of attack runs, and vertical moment in the 
90 ° angle of attack runs. 
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2 . Tabular Presentation 



The complete theoretical and experimental results are 
presented in tabular form in Appendix A. The tables of 
theoretical results are presented directly as they were re- 
ceived from the IBM 1620 computer. The input data, identi- 
fied by the computer program symbols, is printed at the top 
of the page. The output quantities are printed in eleven 
groups corresponding to phase angles from 0 ° to 360 ° in 36 0 
increments. The quantities within each group are identified 
in the format statement included with each table. 

B. Discussion of Results 
1 . Model Results - Theory and E xperiment 

Within the accuracy of the experiments, the measured 
results show general agreement with the results predicted 
by the modified Morrison formulation developed in this re- 
port. In all cases the prediction of phase angle seems very 
good and with the exception of overturning moment, in which 
the experimental values are subject to large errors as dis- 
cussed in Section III-C-3, the quantitative values of forces 
and moments as computed by the modified Morrison formulation 
show fair agreement with the experimental values . Further- 
more, if one assumes that the force distribution on the body 
predicted by theory is as accurate as the magnitude of the 
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total vertical force predicted, then one may conclude that 

the overturning moments predicted by the modified Morrison 

theory are of similar accuracy , 

Although the potential theory and modified Morrison theory 

are completely independent and are based on widely differing 

assumptions, their results show striking agreement for all 

cases except those in which viscous forces become important „ 

One would expect this difference in the theoretical results 

of the two theories in the case where viscous effects are 

relatively large since the potential flow theory does not 

account for viscous effects. This point is illustrated in 

the curves for run number 4 which shows a difference in 

both magnitude and phase angle predicted by the two theories. 

Furthermore, this curve indicates that the modified Morrison 

theory produces closer agreement with experimental results, 

A rough idea of the relative importance of viscous and in- 

u T 

ertia forces for each run is given by the parameter — -g — , 
This parameter can be shown to be proportional to the ratio 
of the maximum values of viscous to inertia forces. 

Irregularities are observed for all runs with the ten- 
foot wave (Huns 1-4) at phase angles of about 270°, These 
irregularities may be attributed to the presence of a sig- 
nificant second harmonic in the profile (Stokian) of this 
wave. This second harmonic produces differences in the 
values of accelerations and velocities encountered in the 
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wave crests and troughs. The results for the five-foot 
wave, which has a negligible second harmonic, do not ex- 
hibit this tendency. This point is illustrated by the 
plot of the moment curve obtained from the sine wave 
potential theory included for run number 4. This curve 
does not have this "hump" since there are no higher 
harmonics and thus the velocities and accelerations en- 
countered in the crests and troughs are of equal magnitude. 

The submarine model experiments were conducted only 
for angles of attack of 0°, 45°, and 90°; however, the 
theoretical forces and moments are included in the tables 
and are plotted on the graphs for the additional cases of 
180°, 235° and 270° to show the effect of body geometry. 
(This was accomplished analytically by simply changing the 
sign of the coefficients for the odd powers of x in the 
fourth degree polynomial describing the body shape). 

It was noted at the time the experiments were con- 
ducted that the torsional gage was not, in fact, completely 
independent of bending for all orientations ( o< ) of the 
gage. For the runs with 45° angle of attack, however, it 
was possible to orient the gage such that bending effects 
were eliminated. For the 90° angle of attack runs, it was 
not possible to accomplish this because of an obstruction 
between the torsion gage and the mounting bracket. Thus, 
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for the submarine shape at 90° (run number 14), bending 
effects are present in the lateral moment values obtained 
experimentally and they must be interpreted with this in 
mind. For runs at 0, of course, the horizontal moment 
is not present. 

2c Prototype Results - Theory Only 

For the purpose of comparing the theories in the 
prototype range, the forces and moments imposed on a 250 
foot long, 52 foot diameter submarine and on a 5 foot 
diameter, 50 foot long cylindrical structural member have 
been calculated for a severe sea state. The wave selected 
(400 ft. long, 25 ft. high) corresponds to a sea wave gen- 
erated by a 58 knot wind (15) « The exact input data and 
results are included in Appendix B. The results of both 
theories show close agreement in both magnitude and phase 
angle for the submarine prototype. In the case of the 
prototype structural member, the theoretical results differ 
markedly in both phase angle and magnitude. This is to be 
expected from the model results since this configuration 
corresponds closely to that of the 1 Inch cylinder in the 
10 foot wave. An examination of the modified Morrison 

theory indicates that the viscous forces are relatively 

« T 

large for this case ( ■ = 11. 85). 
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It is noted that for the prototype structural members 



and for the 1 inch cylinder model (run number 4) the 
potential theory predicted larger forces and moments than 
did the modified Morrison formulation. This was unexpected 
since the fact that the potential theory neglects the large 
viscous effects would lead one to believe that the potential 
theory would result in values of lesser magnitude. This 
difference may be explained by considering the value of C m 
used in the modified Morrison theory. In these particular 
cases , the average C m obtained from Fig, 4 and used in the 
modified Morrison theory is In the range of 0,8. However, 
the "C^" implicit in the potential flow solution for a cir- 
cular cylinder Is 2.0 in all cases. With the exception of 
the 1 inch cylinder (run number 4), the C m determined from 
Fig, 4 is in the range of 2,0 and close agreement between 
the potential theory and the modified Morrison theory Is 
observed. ThuSj, the difference in the theoretical results 
is accounted for by two effects; 

(1) The increased magnitude of viscous 
force accounted for in the modified 
Morrison theory but neglected in 
the potential flow theory. 

(2) The decreased magnitude of inertial 

force due to the reduced C employed 

m 

in the modified Morrison theory as 
compared to constant value of 2.0 im- 
plicit in the potential theory. This 
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difference is present only where 
viscous forces are large and 
likely indicates the previously 
discussed phenomena of flow 
separation on the body shape 
effective in inertia calculations. 

It should be noted at this point that the modified 
Morrison theory did in fact show closer agreement to the 
experimental results in this range. 
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cylinder*, 10 ft* wavelengths angle of attack 0 p axis 10” below still water 
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7 / 8 " cylinder s 5 ft» wavelength* angle of attack 0% axis 6 W below still water 
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7/8” cylinders 5> ft. wavelengths angle of attack 0°, axis 10* 9 below still water 
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7/8" cylinder, 5> ft. wavelength, angle of attack 0®, axis 1 V below still water 
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cylinder’s, 5 ft. wavelengths, angle of attack 0% axis 10" below still water 







cylinder^ 5 fto wavelengths angle of attack axis 10” below still water 
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Submarine shape s $ ft<, wavelengths, angle of attack hZ* » axis 10 M below still water 
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2 7/8 i} cylinder p 5 ft. wavelength^ angle h of attack i;5 @ s> axis 10” below still water 
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Run number l4 (continued) 
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Vo SUMMARY AND CONCLUSIONS 



Ac Summary 

Laboratory measurements of forces and moments exerted 
by water waves on a stationary, submerged, horizontal sub- 
marine shape and two dimensional circular cylinders have 
been mad® at various relative depths and angles of incidence . 
These results have been compared with two theories; 

lo Potential theory of Cummins ( 6 ) 

2o Superposition of inertial and drag 
forces after Morrison ( 5 ), modified by an 
incremental formulation which provides for 
variation in the force coefficients with 
axial body posit ion . 

B. Conclusions 

The comparison of theory and experiment shows; 

1, Both theories agree very closely for predominance of 
inertia forces. 

2. Where inertial forces predominate the predicted phase 
relationship is excellent and the error between theory 
and experiment at maximum force has an average value 
of approximately 19$ for the modified Morrison theory 
and 17 $ for the potential theory. 



3° For the case of significant or predominant viscous forces 
the modified Morrison theory gives a very good prediction 
of phase relationship and yields an error in maximum 
force having an average value of approximately 22 % for 
the modified Morrison theory and 3 2 % for the potential 
flow theory* 

4o Instrumentation difficulties are estimated to contri- 
bute heavily to the apparent error particularly with 
respect to measured wave height and vertical moment* 

5* The true error between theory and experiment is believed 
due predominately to the use of inexact drag and inertia 
coefficients . 

6, Additional experiments are needed to define the time 
variation of both the Inertia and drag coefficients in 
cases of oscillatory flow* 
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VI. RECOMMENDATIONS FOR FURTHER STUDY 



Ac Additional experiments should be run covering a wider 
range of parameters in order to further verify the 
modified Morrison theory. 

Bo Any additional experiments should be run with an 

improved instrument for measuring overturning moments. 
It is recommended that a torsion type gage similar to 
the design presented in this report be employed for 
the direct measurement of overturning moment. 

C. Further study of and C should be made with the 

Dm 

objective of obtaining accurate instantaneous values 
of these coefficients in oscillatory flow. 
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APPENDIX A. TABLES OP RESULTS 



The tables included here contain the experimental and 
theoretical results obtained for each test run. The tables 
are arranged in order of run number starting with run number 
1 and continuing through run number 15. There are four pages 
for each run number. Thus, the first four pages pertain to 
run number 1 and present the results in the following order: 
first page, experimental; second page, modified Morrison 
theory; third page, Stokian potential theory; and fourth 
page, sine wave potential theory. The next four pages per- 
tain to run number 2 and so on through the Appendix. Each 
table is labeled at the top indicating the method and run 
number for quick identification. 

The tables of experimental results are self-explanatory. 
However, as an aid to interpretation of the theoretical 
tables the following example is given. 

Consider the Modified Morrison Solution - Stokian Wave, 

Run Number 1. The wave length (WL) equals 10.0 feet. The 
body length (BL) equals 2.5 feet. The wave height (H) equals 
.650 feet. The angle of attack (ALPHA) equals 0.00001 radians 
^0°. The coefficients of the sixth degree polynomial de- 
scribing the body are as follows: constant term » .5555 feet; 
first power coefficient ■ 0.0; etc. (in this example the 
model was the 4" cylinder). 

* For exceptions, see note on next page. 
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At phase angle 0°, the vertical force * -.65337439 lbs.; 
the vertical moment = .75109979; the ratio of vertical drag 
force to total vertical force (Ratio) * .54955693; 

ETA ® -3.5781317E-02; and celerity - 6 . 8191871 . At phase 

angle 72°, the vertical force - 2.3215859; RATIO - -4.6022410E-02; 

and ETA * . 52853091 ; and celerity » 6 . 8191871 . 

Although there are some slight variations, the format 
of the other tables is essentially the same and should be in- 
terpreted similarly. These tables have been printed directly 
from the computer output and therefore some of the numbers 
are in floating point form. For example: 9 . 66 E-O 6 should 
be interpreted as 9.66 x 10 - ^. 



Note: Run numbers 2, 6 , 12, and 14 have seven pages each 

since the theoretical calculations are included 
for both the stern first and bow first orientations. 
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EXPERIMENTAL RESULTS 



Run Number 1 



Model Angle of 

Attack 


Depth Below 
Surface 


Wave 

Height 


Wave Depth of 

Length Water 


4^ 

o 

M 

o 


0° 


10' 


7 . 56 ” 


10.0' 2.0' 






Vertical 


Vertical 


Axial 


Phase 


ETA 


Force 


Moment 


Force 


Angle 


ft. 


lbs . 


ft .lbs. 


lbs . 


0 


- 0.133 


-1.500 


0.446 


1.150 


36 


0.150 


0.500 


0.550 


1.150 


72 


0.512 


2.500 


0.087 


0.602 


108 


0.522 


2.675 


0.156 


- 0.204 


144 


0.200 


1.475 


0.520 


-1.020 


180 


- 0.055 


- 0.440 


0.048 


- 1.025 


216 


-0.212 


-1 . 920 


- 0.019 


-1.150 


252 


-0.272 


-5.010 


- 0.015 


- 0.527 


288 


-0.278 


-2.900 


0.175 


0.145 


524 


- 0.255 


- 2.590 


0.500 


0.806 


560 


- 0.105 


- 1.275 


0.446 


1.150 
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MODIFIED MORRISON SOLUTION - STOKIAN WAVE 



WL BL 


RUN NUMBER 
H DEPTH 


NO. 1 

L RHO 


Vise ALPHA 


SI S2 


S3 S4 


S5 S6 


S7 


i O o 0 2 « 5 


.630 2.0 


“.8333 1.93 


9.66E-06 .00001 


.3333 0.0 


0.0 0.0 


0.0 0.0 


0.0 


Vertical 


FORMAT 

Lateral 


Vertical 


Lateral 


Force 


Force 


Moment 


Moment 


RATIO 

ETA — CELERITY 

-.65337449 


3. 860437 IE-05 


.75109979 


6. 5102020E-07 


.54955693 
=3. 578 13176=02 

.96905360 


6.8191871 

3.2009790E-05 


.96520920 


-7.5494290E-06 


=.30307991 

.17409580 

2.3215859 


6.8191871 
1.2411 131 E-05 


.48378401 


-1.3067379E-O5 


=4. 6022410E-02 
.32853091 
2.5352703 


6.8191871 
-1. 24086 10E-05 


-.27777548 


-1.3067589E-05 


4.2144374E-02 

.32852988 

1.5564387 


6.8191871 

-3.2007922E-05 


-.75321598 


-7.54991 18E-Q6 


.18870160 

.17409343 

6.4?41400E=02 


6.8191871 

-3.8603406E-05 


-.82226051 


6.5054670E-07 


5.5461631 

-3.5783628E-02 

=1.2178897 


6.8191871 

-3.0453386E-05 


-.52499325 


7.9517624E-06 


=.18619499 

=.19621104 

=1.8881707 


6.8191871 
=1.1 448547 E-05 


-.16734061 


1.20 145 37 E-05 


=3.3927206E=02 

=.27063568 

=2.0162881 


6.8191871 

1.1448860E-05 


4.3490840E-02 


1.201 46 13E-05 


3.1 773741 E=02 
=.27063446 
=1.6713970 


6.8191871 
3.0453931 E-05 


.30233913 


7.9519432E-06 


.13567769 

-.19620692 

=.65334169 


6.8191871 
3. 8604397 E-05 


.75111003 


6.5082500E-07 


.54958608 

-3.5776692E=02 


6.8191871 
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POTENTIAL SOLUTION - STOKIAN WAVE 



RUN NUMBER 1 



WL 


BL 


H 


DEPTH 


z 


RHO 


vise 


ALPHA 


SI 


S2 


S3 


Sk 


S5 


S6 


S7 




10.0 


2.5 


.630 


2.0 


-.8333 


1.93 


9. 66 E -06 


0.0 


.3333 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 





FORMAT 



Vertical 


Vertical 


Lateral 


Lateral 


Force 

AXIAL FORCE 
ETA 


Moment 


Force 


Moment 


-.27774295 

1.7389183 

-3. S 781317 E -02 


.74240273 


.00000000 


.00000000 


1.1916802 
1 .4418401 
.17409580 


,81085038 


.00000000 


.00000000 


2.2957469 

.55899823 

.32853091 


.35934393 


. 00000000 


.00000000 


2.2917395 

-.55901084 

.32852988 


-.35935177 


.00000000 


.00000000 


1 . 1916634 
-1.4418474 
.17409343 


-.81085250 


.00000000 


.00000000 


-.27775891 

-1.7389177 

-3.5783628 E - Q 2 


-.74239937 


,00000000 


.00000000 


-1.3633423 

-1.3717781 

-.19621104 


-.39037691 


.00000000 


.00000000 


-1.8423499 

-.51569103 

-.27063568 


-9.9479440 E -02 


.00000000 


.00000000 


-1.8423424 

.51572614 

-.27063446 


9.9487100 E -02 


.00000000 


.00000000 


-1.3633152 
1.3718012 
-. 19620692 


.39038952 


.00000000 


.00000000 


-.27771099 

1.7389193 

-3.5776692 E - G 2 


.74240911 


.00000000 


.00000000 
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POTENTIAL SOLUTION - SINE WAVE 



RUN NUMBER NO. 1 



WL 


BL 


H 


DEPTH 


z 


RHO 


Vise ALPHA 


SI 


S2 


S3 


Sk 


S5 


S6 


ST 


10.0 


2.5 


.630 


2.0 


-.8333 


1.93 


9.66E-06 0,0 


.3333 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 



FORMAT 



Vertical 

Force 

ETA 



Lateral 

Force 



Vertical Lateral 

Moment Moment 



CELERITY (Last Line Only) 



,00000000 

,00000000 


,00000000 


.97663094 


,00000000 


1,6805604 

.18515273 


,00000000 


.79011016 


.00000000 


2.7192010 

.29958309 


.00000000 


.30179281 


.00000000 


2,7191944 

.29958238 


.00000000 


-.30179960 


,00000000 


1 .6805434 
.18515086 


.00000000 


-.7901 1431 


.00000000 


-2.0960Q00E-Q5 
“2,31 14700E-06 


,00000000 


-.97663080 


,00000000 


”1 .6805774 
“.18515460 


,00000000 


“.79010587 


.00000000 


“2.7192075 

“.29958380 


,00000000 


-.30178599 


« 00000000 


“2,7191880 

“,29958166 


.00000000 


,30180641 


,00000000 


“1,6805264 

“.18514898 


.00000000 


.7901185! 


,00000000 


4.1970000E-05 

4,.62514SQE“06 
7, 1558506 


, 00000000 


.97663081 


,00000000 
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EXPERIMENTAL RESULTS 



Run Number 2 



Model 


Angle of 
Attack, 


Depth Below 
Surface 


Wave 

Height 


Wave 

Length 


Depth of 
Water 


Sub 


0° 


10" 


7.56" 


10.0' 


2.0' 



Phase 

Angle 


ETA 

ft. 


Vertical 
Force 
lbs . 


Vertical 
Moment 
ft. lbs. 


Axial 
Force 
lbs . 


0 


-0.095 


-0.570 


0.517 


1.470 


56 


0.150 


1.400 


0.154 


1.290 


72 


0.508 


2.750 


0.154 


0.462 


108 


0.522 


2.840 


0.212 


-0.557 


144 


0.208 


1.520 


0.202 


-1.545 


180 


- 0.060 


- 1.090 


-0.241 


-1.654 


216 


-0.218 


-2.650 


-0.154 


-1 . 240 


252 


- 0.262 


-5.490 


-0.212 


-0.442 


288 


-0.279 


-5.050 


-0.048 


0.500 


524 


-0.222 


-2.150 


0.165 


1.085 


560 


-0.055 


-0.525 


0.546 


1.460 
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MODIFIED MORRISON SOLUTION - STOKLAN WAVE 



RUN NUMBER NO. 2 



WL 


BL 


SI 


S2 


10.0 


2.2917 


.4725 


-.1319 



Vertical 

Force 

RATIO 

ETA - CELERITY 

- 1 .0656721 

.36256475 

-3.5781317E-02 
.97319690 
-=.34366181 
.17409580 
2.9001994 
“3.6396373E-02 
.32853091 
3.3017310 
2.1 804841 E-02 
.32852988 
2.1029941 
.14030865 
. 1 7409343 
.16700561 
2.3960077 
-3.5783628E-02 
-1.4852996 
-.16403790 
-.19621104 
-2.2956511 
-2.5290965E-02 
=.27063568 
-2.4445476 
1.5987789E-02 
-.27063446 
-2.1540432 
9.4473629E-02 
-.19620692 
-1.0656342 
.36257994 
-3.5776692E-02 



H DEPTH 

S3 S4 

.630 2.0 

-.1529 .0854 

FORMAT 

Lateral 

Force 



4.8203990E-05 

6.8191871 

4.2401360E-05 

6.8191871 

1.9264074E-05 

6.8191871 

-1.2130750E-05 

6.8191871 

-3.8308826E-05 

6.8191871 
-4.8594197E-05 

6.8191871 
-4.0122369E-05 

6.8191871 

-1 .7463842E-05 

6.8191871 

1.0966052E-05 

6.8191871 
3.579051 IE-05 

6.8191871 
4.8204080E-05 

6.8191871 



Z RHO 

S5 s6 

“.8333 1.93 

-.1529 0.0 

Vertical 

Moment 

.45588950 

.36258410 

-8.2453650E-02 

“.51195430 

-.59806700 

-.41479092 

-7.6636613E-02 

.19791034 

.30323009 

.36423571 

.45589055 



Vise ALPHA 

S7 

.66E-06 .00001 

.0 

Lateral 

Moment 

-4.9670603E-06 
-8.4332442E-06 
-8.6852926E-06 
-5.2421597E-06 
4.4381 820E-07 
5.7312159E-Q6 
8.447401 6E-06 
7.9298694E-06 
4.7610783E-06 
1 .4199920E-08 
-4. 9671 660 E -06 
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POTENTIAL SOLUTION - STOKIAN WAVE 



RUN NUMBER NO. 2 



WL BL 


H DEPTH 


Z RHO 


vise 


ALPHA 


SI S2 


S3 S4 


S5 S6 


S7 




10.0 2.2917 


.630 2.0 


-.8333 1.93 


9. 66 E ' 


-06 0,0 


.4725 -.3319 


-.1529 ,0854 -.1529 0.0 

FORMAT 


0.0 




Vertical 


Vertical 


Lateral 




Lateral 


Force 

AXIAL FORCE 
ETA 


Moment 


Force 




Moment 


-.63008841 

2.1713363 

-3.5781317 E -02 


.43358893 


.00000000 




.00000000 


1.2429275 

1.9099326 

.17409580 


.27736616 


.00000000 




.00000000 


2.8367305 

.86769309 

.32853091 


12048792 


.00000000 




. 00000000 


3.0394626 

-.54649164 

.32852988 


-.51652822 


.00000000 




.00000000 


1.6956097 
-1.7256698 
. 1 7409343 


-.60690734 


.00000000 




.00000000 


-.22669670 

-2.1889513 

-3.5783628 E -02 


-.35428654 


.00000000 




.00000000 


-1.6340276 

-1.8073175 

-.19621104 


-5.9909950 E -03 .00000000 




. 00000000 


-2.2216734 

-,78665051 

-.27063568 


,20890660 


.00000000 




.00000000 


-2.2682472 

.49397807 

-.27063446 


.29980174 


.00000000 




.00000000 


-1,8340290 

1.6121879 

-.19620692 


.38454759 


,00000000 




.00000000 


-.63005037 

2,1713403 


.43358829 


.00000000 




.00000000 



-3.5776692 E -02 
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POTENTIAL SOLUTION - SINE WAVE 
RUN NUMBER NO. 2 



vise 

9. S ^6E 

0.0 



FORMAT 



Vertical Lateral 


Vertical 


Force Force 


Moment 


ETA 




CELERITY (Last Line Only) 





WL 


BL 


H 


DEPTH 


z 


RH 0 


SI 

HO . 0 


S 2 

2.2917 


S 3 

.630 


CO • 
CM 


S 5 

-.8333 


S 6 

1.93 


.4725 


•=.1319 


“.1529 


.0854 


-.1529 


0.0 



=.26534472 

.00000000 


.00000000 


.51822792 


1.8923101 

.18515273 


.00000000 


.18638054 


3.3271637 

.29958309 


.00000000 


-.21665824 


3.4911479 

.29958238 


.00000000 


=.53694054 


2.3216263 

.18515086 


.00000000 


-.65212885 


.26531839 

-2.31147 Q 0 E -06 


.00000000 


-.51822497 


=1.8923326 

=.18515460 


.00000000 


-.18637592 


-3.3271738 

=.29958380 


.00000000 


.21666276 


=3.4911417 

=.29958166 


.00000000 


.53694323 


=2.3216062 

=.18514898 


.00000000 


.65212873 


=.26529207 

4.6251450 E =06 

7.1558506 


.00000000 


.51822206 



ALPHA 
06 0.0 



Lateral 

Moment 



.00000000 

.00000000 

.00000000 

.00000000 

.00000000 

.00000000 

.00000000 

.00000000 

.00000000 

.00000000 

.00000000 
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MODIFIED MORRISON SOLUTION - ST OKI AN WAVE 



RUN NUMBER 2 



WL BL 


H DEPTH 


Z RHO VISC 


SI S2 


S3 S4 


S5 S6 S7 


10 . 2.2917 


.630 2.0 -.1 


B333 1.93 9.66E 


.4726 .1319 


-.1529 -.0854 


-.1529 0.0 0 




FORMAT 


Vertical 


Lateral 


Vertical 


Force 


Force 


Moment 


RATIO 


ETA 

A rv oocoooooo 


. CELERITY 




=.63352014 


4. 861 8656E-05 


.31010247 


.63181168 


-3.5781317E-02 


6.8191871 




1.5136521 


3.8329450E-05 


.67792726 


=.19499856 


.17409580 


6.8191871 




3.1592984 


1 .2139486E=05 


.59913571 


=2. 279687 6E=02 


.32853091 


6.8191871 




3.1128327 


-1.9270587E=05 


.21149196 


3 0 3921 681 E=02 


.32852988 


6.8191871 




1.6428603 


=4. 2419961E-05 


-.18511827 


.20363856 


.17409343 


6.8191871 




=.29311979 


=4.8226590E=05 


-.44218832 


=1.3185452 


=3.5783628E~02 


6.8191871 




=1.7479567 


-3.5807553E-05 


“.45458902 


=.11645793 


=.19621104 


6.8191871 




=2.3675716 


=1 . 0971372E-05 


-.35214078 


=1.651 341 5E-02 


=.27063568 


6.8191871 




=2.4129613 


1 .7472348E-05 


“.27155478 


2.4071234E=02 


=.27063446 


6.8191871 




=1.9735026 


4.0142184E-05 


-9.3127519E-02 


.12349923 


=.19620692 


6.8191871 




=.63347626 


4.861 8642E-05 


.3101 1323 


.63185651 


=3.5776692E=02 


6.8191871 





ALPHA 

i 

■06 .00001 
0 

Lateral 

Moment 

5-7335997E-06 
4.4308343E-07 
■5. 245 87 02 E-0 6 
8.6904294E-06 
8.4375908E-06 
4.9689347E-06 
1.5391470E-08 
4.7646544E-06 
7 . 93463 1 6E-06 
8. 45 17525E-06 
5. 7335007E-06 
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POTENTIAL SOLUTION - ST OKI AN WAVE 
RUN NUMBER 2 



WL BL 


H DEPTH 


Z RHO 


Vise ALPHA 


SI S2 


S3 S4 


S5 . S6 


S7 


10. 2.2917 


.630 2.0 


8333 1.93 


9 . 66 E-O 6 0.0 


.4726 .1319 


-.1529 -.0854 -.1529 0. 

FORMAT 


0 0.0 


Vertical 


Vertical 


Lateral 


Lateral 


Force 

AXIAL FORCE 
ETA 


Moment 


Force 


Moment 


“.22680704 
2.1900142 
“3 .5781317 E— -02 


.35453098 


.00000000 


.00000000 


1.6964243 

1.7265114 

.17409580 


.60724320 


.00000000 


.00000000 


3.0409256 

.54676452 

.32853091 


.51675324 


.00000000 


.00000000 


2.8381064 

-.86810571 

.32852988 


.12046094 


.00000000 


.00000000 


1.2435468 

-1.9108547 

.17409343 


-.27759480 


.00000000 


.00000000 


=.63037451 

=2.1723918 

“3o5783628E-02 


-.43385205 


.00000000 


.00000000 


-1.8349311 

-1.6129563 

-.19621104 


-.38473321 


.00000000 


.00000000 


-2.2693604 

-.49419633 

-.27063568 


-.29991471 


.00000000 


.00000000 


-2.2227551 

.78705478 

-.27063446 


-.20895651 


.00000000 


.00000000 


-1.6347997 

1.8082097 

-.1Q6206Q? 


6.0669590E-03 


.00000000 


.00000000 


-.2267651! 

2. 1900135 
-3.5776692E-02 


.35453917 


.00000000 


.00000000 
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POTENTIAL SOLUTION - SINE WAVE 



RUN NUMBER 2 



WL BL 


H DEPTH Z 


RHO Vise 


ALPHA 


SI S 2 


S 3 S 4 S 5 


S 6 S 7 




10 . 2.2917 


.630 2.0 -.8333 


1.93 9 . 66 E - 


-06 0.0 


.4726 .1319 


-.1529 -.0854 


-.1529 0.0 


0.0 




FORMAT 




Vertical 


Lateral 


Vertical 


Lateral 


Force 

ETA 


Force 


Moment 


Moment 


CELERITY (Last Line Only) 






. 26543406 
.,00000000 


.00000000 


.51855884 


.00000000 


2.3227435 


„ .00000000 


.65247326 


.00000000 


.18515273 

3.4928399 


.00000000 


.53716397 


.00000000 


.29958309 

3.3287841 


.00000000 


.00000000 


.21667537 


.29958238 

1.8932403 

.18515086 


.00000000 


.00000000 


-.18657623 


=-.26546036 


. 00000000 


-.51856171 


.00000000 


“ 2 . 3114700 E “ Q 6 

“ 2.3227638 

“.18515460 


.00000000 


-.65247338 


.00000000 


“ 3.4928463 

=.29958380 


.00000000 


-.53716124 


.00000000 


= 3.3287743 

“.29958166 


.00000000 


-.21667086 


.00000000 


“ 1. 8932177 
“.18514898 


.00000000 


.18658081 


.00000000 


.26548670 


.00000000 


.51856463 


.00000000 



4 . 6251450 E -06 

7.11558506 



-103- 



EXPERIMENTAL RESULTS 



Run Number 5 



Model 


Angle of 
Attack 


Depth Below 
Surface 


Wave 

Height 


Wave 

Length 


Depth of 
Water 


2 - 7 / 8 ” 


cyl. 0° 


10" 


7 . 56 ” 


10.0' 


2.0’ 



Phase 

Angle 


ETA 

ft. 


Vertical 
Force 
lbs . 


Vertical 
Moment 
ft. lbs. 


Axial 
Force 
lbs . 


0 


0.0 


-0.290 


1.468 


.762 


56 


.190 


0.280 


1.280 


.618 


72 


.550 


0.720 


0.656 


.276 


108 


.520 


0.750 


- 0.145 


.552 


144 


.210 


0.588 


- 0.825 


-.695 


180 


-.081 


-0.170 


- 1.420 


-.809 


216 


-.220 


-0.670 


-0.800 


-.685 


252 


-.290 


- 0.990 


- 0.254 


-.209 


288 


-.272 


-0.980 


0.540 


.209 


524 


-.177 


-0.720 


1.152 


.666 


560 


.010 


-0.260 


1.448 


.762 
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MODIFIED MORRISON SOLUTION - STOKIAN WAVE 



RUN NUMBER NO. 3 



WL 


BL 


SI 


S2 


10.0 


2.5 


.23958 


0.0 



Vertical 
Force 
RATIO 
ETA - CELERITY 

“.3920360? 
.62068862 
-3.5781317E-02 
.43882926 
-.45392855 
.17409580 
1.1545002 
“6.2800054E-02 
.32853091 
1.2995029 
5.5793746E-02 
.32852988 
.83721586 
.23792960 
.17409343 
9.4619590EHD2 
2.5716860 
“3.5783628E-02 
=.57602820 
-.26718520 
=.19621104 
-.94286363 
-4.6169795E-02 
-.27063568 
-1 .0299263 
4.2269975E-02 
-.27063446 
-.88383029 
.17414023 
-.19620692 
-.39201959 
.62071625 
-3.5776692E-02 



H DEPTH 

S3 S4 

.630 2.0 

0.0 0.0 

FORMAT 

Lateral 

Force 

1.9946700E-05 

6.8191871 

1.6539560E-05 

6.8191871 
6.41 32741 E-06 

6.8191871 

-6.4108394E-06 

6.8191871 
-1.653775 IE-05 

6.8191871 

-1.9945789E-05 

6.8191871 

-1.5734955E-05 

6.8191871 

-5.9154609E-06 

6.8191871 
5. 9 1540 17 E-06 

6.8191871 

1.5735248E-05 

6.8191871 
1.994671 IE-05 

6.8191871 



Z RHO 

S5 s6 

-.8333 1.93 

0.0 0.0 

Vertical 

Moment 

.37339746 

.50582432 

.26227931 

-.12251047 

-.36243717 

-.42156760 

-.28428180 

-9.5429470E-02 

1 . 10961 80E-02 

.13374111 

.37340348 



Vise ALPHA 

S7 

9.66E-06 .00001 

0.0 

Lateral 

Moment 

3.3649450E-07 
-3. 9006067 E-06 
-6.7517396E-06 
-6.7519385E-06 
-3.9010721 E-06 
3.3601 190E-07 
4.1085050E-06 
6. 2077567E-06 
6.2078659E-06 
4. 108781 2E-06 
3. 3639420E-07 



POTENTIAL SOLUTION - STOKIAN WAVE 
RUN NUMBER NO. 3 



WL 


BL 


H 


DEPTH 


z 


RHO 


vise 


ALPHA 


SI 


S2 


S3 


S4 


S5 


S6 


ST 




10.0 


2.5 


.630 


2.0 


-.8333 


1.93 


9.66 E -06 


0.0 


.23958 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 





FORMAT 



Vertical 


Vertical 


Lateral 


Lateral 


Force 

AXIAL FORCE 
ETA 


Moment 


Force 


Moment 


-.14350714 
.89848263 
=3.5781 3 17 E -02 


.38359240 


.00000000 


.00000000 


.61572987 

,74498511 

.17409580 


.41895870 


.00000000 


.00000000 


1.1841239 

.28882908 

.32853091 


.18566960 


.00000000 


.00000000 


1.1841201 

-.28883558 

.32852988 


-.18567363 


.00000000 


.00000000 


.61572125 

-.74498891 

.17409343 


-.41895981 


.00000000 


.00000000 


=. 14351541 

-.89848235 

-3.5783628 E -02 


-.38359071 


.00000000 


.00000000 


-.70442608 
-.70878479 
19621 104 


-.20170405 


.00000000 


.00000000 


-.95192484 

-.26645266 

-.27063568 


-5.1400060 E -02 


.00000000 


.00000000 


-.95192093 

.26647081 

-.27063446 


5.1404050 E -02 


.00000000 


.00000000 


-.70441208 

.70879671 

-.19620692 


,20171057 


.00000000 


.00000000 


-. 14349064 
.89848320 
-3.5776692 E -02 


.38359573 


.00000000 


.00000000 
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POTENTIAL SOLUTION - SINE WAVE 



RUN NUMBER 3 



WL 


BL 


H 


DEPTH 


z 


RHO 


vise 


ALPHA 


SI 


S2 


S3 


S4 


S5 


s6 


S7 




10.0 


2.5 


.630 


2.0 


-.8333 


1.93 


9.66 E -06 


0.0 


.23958 0.0 


0.0 


0.0 


0.0 


0.0 


0.0 







FORMAT 






Vertical 


Lateral 


Vertical 


Lateral 


Force 


Force 


Moment 


Moment 


ETA 








CELERITY (Last Line Only) 






.00000000 


.00000000 


.50461590 


.00000000 


.00000000 

.86832965 


.00000000 


.40824242 


.00000000 


.18515273 

1.4049854 


.00000000 


.15593347 


.00000000 


.29958309 

1.4049820 


.00000000 


-.15593699 


.00000000 


.29958238 

.86832093 

.18515086 


.00000000 


=.40824456 


.00000000 


-1.0840000 E -05 

=2.3114700 E =06 


.00000000 


=.50461586 


.00000000 


=. 8683384? 
=.18515460 


.00000000 


-.40824020 


.00000000 


=1.4049887 

=.29958380 


.00000000 


-.15592996 


.00000000 


=1.4049786 

=.29958166 


.00000000 


.15594051 


.00000000 


=.86831212 

=.18514898 


.00000000 


.40824674 


.00000000 


2.1690000 E =05 


.00000000 


.50461584 


.00000000 



4.6251450 E -06 

7.1558506 



-107- 



EXPERIMENTAL RESULTS 



Run Number 4 



Model 


Angle of 
Attack 


Depth Below 
Surface 


Wave 

Height 


Wave 

Length 


Depth of 
Water 


1" cyl. 


0° 


10" 


9.6" 


10' 


2.0' 



Phase 

Angle 


ETA 

ft. 


Vertical 

Force 

lbs. 


Vertical 
Moment 
ft. lbs. 


Axial 
Force 
lbs . 


0 


0.045 


- 0.080 


0.145 


0.085 


36 


0.255 


-0.009 


0.159 


O.O 65 


72 


0.428 


0.060 


0.100 


0.018 


108 


0.450 


0.090 


-0.024 


- 0 . 058 


144 


0.150 


0.082 


-O.O 69 


-0.097 


180 


- 0.190 


0.055 


-0.144 


-0.099 


216 


- 0.518 


0.010 


- 0.126 


-0.073 


252 


-0.540 


-0.057 


-0.049 


- 0.022 


288 


- 0.290 


-0.110 


0.022 


0.055 


524 


-0.205 


- 0.128 


0.091 


0.075 


560 


0.002 


- 0.094 


0.137 


0.087 



- 108 - 



MODIFIED MORRISON SOLUTION - STOKIAN WAVE 



RUN NUMBER 4 



WL BL 

SI S2 

10.0 2.5 

.0833 0.0 



Vertical 

Force 

RATIO 

ETA - CELERITY 

“.13703102 
.91883421 
-5.7697260E-02 
-6. 9 7 46400 E=0 2 
1.5283408 
.21728531 
3.4180245E-02 
= 1.1703612 
.42710123 
.11418739 
.35033788 
.42709983 
.14344634 
.74311432 
.21728213 
.11478574 
1 .0969000 
=5.77001956=02 
3.3713288E=02 
2.2969364 
=.25294542 
=3.4872173E=02 
=.61123604 
■=.33374512 
=7.7503858E=02 
.27504095 
=.33374389 
=.12116277 
.63913638 
=.25294070 
=.13703047 
.91884118 
=5.7691390E=02 



H DEPTH 

S3 S4 

.800 2.0 

0.0 0.0 

FORMAT 

Lateral 

Force 

3.1822275E-06 

6.9512435 

2.6565185E-06 

6.9512435 

1.0349533E-06 

6.9512435 

-1.0320273E-06 

6.9512435 

-2.6543086E-06 

6.9512435 
-3.1811 167E=06 

6.9512435 
=2.49261 86E-06 

6.9512435 

-9.3317870E-07 

6.9512435 

9.3261365E-07 

6.9512435 

2.4926799E-06 

6.9512435 
3. 1822300E-06 

6.9512435 



Z RHO 

S5 S6 

-.8333 1.93 

0.0 0.0 

Vertical 
Moment 

8.0940480E=03 
6.3173469E-02 
5.1245906E-02 
2.6370342E-02 
8.4539 180E-03 
-3.8730379E-02 
-4.8941810E-02 
-2.2958604E-02 
=1 .8896064E-02 
=2.78981 88E-02 
8.0954960E-03 



Vise ALPHA 

S7 

9.66E-06 .00001 

0.0 

Lateral 
Moment 

6.8420640E=08 
-6.1749196E-07 
-1 .0885991E-06 
-1 .0888208E-06 
-6.1802774E-07 
6.7796940E-08 
6.5960184E-07 
9.7841740E-07 
9.7861 106E-07 
6.60101 12E-07 
6.8404630E-08 
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POTENTIAL SOLUTION - STOKIAN WAVE 



RUN NUMBER NO. k 



WI, 


BL 


H 


DEPTH 


z 


RHO 


vise 


SI 


S2 


S3 


S4 


S5 


S6 


S7 


10.0 


2.5 


.800 


2,0 


-.8333 


1.93 


9.66E 


.0833 


0.0 


0,0 


0.0 


0.0 


0.0 


0.0 



FORMAT 



Vertical 


Vertical 


Lateral 


Force 

AXIAL FORCE 
ETA 


Moment 


Force 


-2.9068397E-02 
. 14332031 
“5.7697260E-02 


6. 1 1 88264E-02 


. 00000000 


9.6308583E-02 

.11961437 

.21728531 


7. 150531 IE-02 


.00000000 


.19388138 

4.6553658E-02 

.42710123 


3.2506532E-02 


.00000000 


.19388071 
“4. 6554701 E-02 
.42709983 


-3.2507236E-02 


.00000000 


9.63071 10E-02 
-. 1 1961497 
.21728213 


-7.1505466E-02 


.00000000 


-2.9069716E-02 
-.14332026 
-5. 77 00 195 E-0 2 


-6. 1 187916E-02 


.00000000 


-.11427434 

-.11228190 

-.25294542 


-2.7498898E-02 


.00000000 


-.14684767 
-4. 202 1420 E-02 
“.33374512 


-5.3094010E-03 


.00000000 


-.14684721 

4.2024287E-02 

-.33374389 


5.3098540E-03 


.00000000 


-.11427236 

.11228381 

-.25294070 


2.7500006E-02 


.00000000 


-2.9065767E-02 


6. 1 1 88935E-02 


.00000000 



.14332043 

=5.7691390E~02 



ALPHA 
06 0.0 

Lateral 

Moment 

.00000000 
.00000000 
.00000000 
.00000000 
.00000000 
. 00000000 
.00000000 
.00000000 
.00000000 
.00000000 
.00000000 
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POTENTIAL SOLUTION - SINE WAVE 
RUN NUMBER NO. 4 



WL 


BL 


H 


DEPTH 


SI 


S2 


S3 


S4 


10.0 


2.5 


.800 


2.0 


.0833 


0.0 


0.0 


0.0 



z 


RHO 


vise 


ALPHA 


S5 


S6 


S7 




8333 


1.93 


9.66E-06 


0.0 


0 


0.0 


0.0 





Vertical 

Force 

ETA 



FORMAT 

Lateral 

Force 



Vertical 

Moment 



Lateral 

Moment 



CELERITY - (Last Line Only) 



.00000000 
. 00000000 


. 00000000 


.13329779 

.23511458 


.00000000 


.21568014 

.38042297 


.00000000 


.21567961 

.38042206 


.00000000 


.13329644 

.23511220 


.00000000 


=1 . 6640000E-06 
-2.935200QE-06 


.00000000 


=. 13329914 
=.23511695 


.00000000 


=.21568064 

=.38042388 


.00000000 


=.21567909 
=.380421 15 


.00000000 


=.13329509 

=.23510982 


.00000000 


3.3290000E-06 

5.8732000E-06 

7.1558506 


.00000000 



7.7463886E-02 


.00000000 


6.2669533E-02 


. 00000000 


2.3937438E-02 


.00000000 


-2.3937983E-02 


.00000000 


“6.2669863E-02 


.00000000 


-7.7463879E-02 


.00000000 


-6.2669194E-02 


. 00000000 


-2.3936898E=02 


.00000000 


2.3938523E-02 


.00000000 


6. 26701 97E-02 


.00000000 


7.7463881E-02 


. 00000000 
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EXPERIMENTAL RESULTS 



Run Number 5 



Model Angle of Depth Below Wave Wave Depth of 
Attack Surface Height Length Water 

4" cyl. 0° 10" 4.0" 5.0' 2.0' 



Phase 

Angle 


ETA 

ft. 


Vertical 

Force 

lbs. 


Vertical 
Moment 
ft . lbs . 


Axial 

Force 

lbs. 


0 


-.070 


-0.175 


1.850 


0.679 


36 


.040 


0.550 


1.560 


0.547 


72 


.158 


0.970 


O.650 


0.264 


108 


.168 


1.120 


-0.284 


-0.0945 


144 


.060 


0.640 


-1.158 


-0.472 


180 


-.055 


-O.065 


-1.495 


-0.604 


216 


-.118 


-0.735 


-1.335 


-0.475 


252 


-.155 


-1.100 


-0.619 


-0.208 


288 


-.155 


-1.100 


0.40j2 


0.160 


324 


-.148 


-0.740 


1.138 


0.472 


560 


-.070 


-0.150 


1.415 


0.688 
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MODIFIED MORRISON SOLUTION - STOKLAN WAVE 



RUN NUMBER NO. 5 



WL BL 

SI S2 

5.0 2.5 

.3333 0.0 



Vertical 

Force 

RATIO 

ETA - CELERITY 



-5.571291 IE-02 
1.0001315 
-2.9928239E-03 
.74230871 
-6.7774962E-02 
9.6941615E-02 
1.2598317 
-1.797 159 IE-02 
.16077232 
1.3051116 
1.7348490E-Q2 
.16077191 
.84292078 
5.9685419E-02 
9.6940590E-02 
5.5717684E-02 
1 .0000466 
-2. 994045 6E-03 
-.74229238 
-6.78Q3181E-02 
-9.8792210E-02 
=1.2598694 
-1.7952980E-02 
=.15593016 
=1.3050985 
1.7331879E-02 
-.15592911 
-.84292842 
5.9708893E-02 
-9.8789370E-02 
-5.5693138E-02 
1.0004865 
-2.9903792E-03 



H DEPTH 

S3 S4 

.333 2.0 

0.0 0.0 



FORMAT 

Lateral 

Force 



1.5377233E-05 

5.1401008 

1.2440683E-05 

5.1401008 

4.7521082E-06 

5.1401008 

-4.7516532E-06 

5.1401008 

-1.2440323E-05 

5.1401008 

-1.5377035E-05 

5.1401008 
-1.24401 7 6E-05 

5.1401008 

-4.7516563E-06 

5.1401008 

4.7518108E-06 

5.1401008 

1.2440353E-05 

5.1401008 

1.5377235E-05 

5.1401008 



Z RHO 

S5 s6 

-.8333 1.93 

0.0 0.0 

Vertical 

Moment 

1 .0730748 
.90095644 
.38122449 
-.28825441 
-.84558410 
-1.0734108 
-.89111583 
-.37470876 
.28241676 
.83543673 
1.0730756 



Vise ALPHA 

S7 

9.66E-06 .00001 

0.0 

Lateral 

Moment 

-8.9955500E-08 
-7.2213753E-06 
-1 . 1566617E-05 
-1.1 566681 E-05 
-7.2215477E-06 
-9.0169010E-08 
7.1657602E-06 
1.17 12441 E-05 
1.1712450E-05 
7.1657824E-06 
-9.0135050E-08 
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POTENTIAL SOLUTION - STOKIAN WAVE 



RUN NUMBER NO. 5 



WL 


BL 


H 


DEPTH 


z 


RHO 


vise 


ALPHA 


SI 


S2 


S3 


S4 


S5 


S6 


S7 




5.0 


2.5 


.333 


2.0 


-.8333 


1.93 


9.66E-06 


0.0 


.3333 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 





FORMAT 



Vertical 


Vertical 


Lateral 


Lateral 


Force 

AXIAL FORCE 

I.TA 


Moment 


Force 


Moment 


6.7400000E-06 
. 69266406 


.99103079 


.00000000 


.00000000 


-2.9928241 E-03 








.73190297 

.56038181 


.80637427 


. 00000000 


.00000000 


9.6941615E-02 

1.1842336 


.30909443 


. 00000000 


.00000000 


.21404637 

.16077232 

1.1842309 


-.30910137 


.00000000 


.00000000 


-.21405121 








.16077192 

.73189555 


-.80637854 


.00000000 


.00000000 


-.56038478 

9.6940590E-02 

-2.3700000E-06 

-.69266407 

-2.9940458E-03 


-.99103072 


.00000000 


.00000000 


-.73190621 

-.56036805 

-9.8792210E-02 


-.79714128 


.00000000 


.00000000 


-1.1842475 

-.21403488 


-.30338394 


.00000000 


.00000000 


-.15593016 

-1.1842390 


.30340452 


.00000000 


.00000000 


.21404938 








-.15592911 

-.73188400 

.56037700 

-9.8789370E-02 


.79715417 


.00000000 


.00000000 


2.5020000E-05 


.99103096 


.00000000 


.00000000 



.69266409 

-2.9903794E-03 
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POTENTIAL SOLUTION - SINE WAVE 



RUN NUMBER NO. 5 



WL 


BL 


H 


DEPTH 


z 


RHO 


vise 


ALPHA 


SI 


S2 


S 3 


S 4 


S 5 


S 6 


S 7 




5.0 


2.5 


.333 


2.0 


-.8333 


1.93 


9.66 E -06 


0.0 


.3333 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 







FORMAT 






Vertical 


Lateral 


Vertical 


Lateral 


Force 


Force 


Moment 


Moment 


ETA 








CELERITY (Last Line Only 






.00000000 

.00000000 


.00000000 


1.0076330 


.00000000 


.74416202 

9.7866442 E -02 


.00000000 


.81519140 


.00000000 


1.2040781 


.00000000 


.31137293 


.00000000 


.15835106 

1.2040753 


.00000000 


-.31137994 


.00000000 


.15835068 

.74415453 


.00000000 


-.81519575 


.00000000 


9.7865453 E - G 2 
=9 . 280 Q 000 E -06 
-1 .2217770 E -06 


.00000000 


-1.0076330 


. 00000000 


-.74416955 
“9. 7867431 E - Q 2 


.00000000 


-.81518706 


.00000000 


“1.2040810 

15835144 


.00000000 


-.31136588 


.00000000 


-=1.2040724 

“.15835030 


. 00000000 


.31138699 


.00000000 


=.74414699 
=9. 7864463 E “02 


.00000000 


,81520009 


.00000000 


1 .8600000 E -05 
2.4447195 E -06 
5.0599504 


.00000000 


1.0076330 


.00000000 
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EXPERIMENTAL RESULTS 



) 

Run Number 6 



Model 


Angle of 
Attack 


Depth Below 
Surface 


Wave 

Height 


Wave 

Length 


Depth of 
Water 


Sub 


0 ° 


10 " 


4.0” 


5.0' 


2 . 0 ' 



Phase 

Angle 


ETA 

ft. 


Vertical 
Force 
lbs . 


Vertical 
Moment 
ft. lbs. 


Axial 

Force 

lbs. 


0 


-.022 


-0.040 


1.592 


1.06 


56 


.040 


1.020 


1.150 


0.802 


72 


.161 


1.640 


0.570 


0.198 


108 


.181 


1.500 


-0.014 


-0.556 


144 


.085 


0.750 


-0.717 


- 1.058 


180 


-.008 


0.545 


-1.185 


-1.114 


216 


-.090 


- 1 . 260 


-0.915 


-0.778 


252 


-.150 


- 1.690 


-0.445 


-0.189 


288 


-.148 


- 1.600 


0.291 


0.587 


524 


-.155 


- 1.080 


1.000 


0.955 


560 


-.060 


- 0.150 


1.480 


1.114 
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MODIFIED MORRISON SOLUTION - STOKIAN WAVE 



RUN NUMBER NO. 6 



WL 


BL 


SI 


S2 


5.0 


2.2917 


.4725 


-.1319 



H DEPTH 

S3 S4 

.333 2.0 

-.1529 .0854 



Z 


RHO 


S5 


S6 


-.8333 


1.93 


-.1529 


0.0 



Vertical 

Moment 

.63397782 

.42649336 

4.0420009E-02 

-.36022288 

-.61203425 

-.63443794 

-.41643671 

-3.3762478E-02 

.35427795 

.60174173 

.63397568 



Vertical 

Force 

RATIO 

ETA - CELERITY 

-.40725375 
.18464792 
-2.9928239E-03 
.89246210 
-7.5641479E-02 
9.6941615E-02 
1.8580283 
-1 .7233900E-02 
.16077232 
2.1113068 
8.0573737E-03 
.16077191 
1.5458420 
3.5188132E-02 
9.6940590E-02 
.39371292 
.19129793 
-2.9940456E-03 
-.90219443 
-7.4545858E-02 
-9.8792210E-02 
-1.8505985 
-1.71 10701 E-02 
-.15593016 
-2.0968439 

8.0280167E-03 
-.15592911 
-1.5444870 
3.5018342E-02 
-9. 8789370E-02 
-.40722337 
.18466219 
-2.99O3792E-03 



FORMAT 

Lateral 

Force 



2.3764159E-05 

5.1401008 
2. 1253247E-05 

5.1401008 

1.0720330E-05 

5.1401008 
3.8011 795 E-06 

5.1401008 

-1.6901116E-05 

5.1401008 

-2.3670300E-05 

5.1401008 
-2, 1445010E-05 

5.1401008 

-1.0932372E-05 

5.1401008 

3.8623144E-06 

5.1401008 
1 .7151 367 E— 05 

5.1401008 
2.376421 2E-05 

5.1401008 



VISC ALPHA 

S7 

1.66E-O6 .00001 
i.O 

Lateral 

Moment 

-2.2867374E-06 
-6.0284617E-06 
-7.4384296E-06 
-6.080881 IE-06 
-2. 47 52897 E-06 
2.1033853E-06 
5.9703338E-06 
7. 5858821 E-06 
6.2301523E-06 
2.4200456E-06 
-2.2868433E-06 
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POTENTIAL SOLUTION - STOKIAN WAVE 







RUN 


NUMBER 


NO. 6 








WL 


BL 


H DEPTH 


Z 


RHO 


vise 


ALPHA 


SI 


S2 


S 3 


S 4 


S 5 


s 6 


S 7 




5.0 


2.2917 


.333 


2.0 


-.8333 


1.93 


9 . 66 E -06 


0.0 


.4725 


=.1319 


=.1529 


.0854 


-.1529 


0.0 


0.0 





FORMAT 



Vertical 


Vertical 


Lateral 


Lateral 


Force 

AXIAL FORCE 
ETA 


Moment 


Force 


Moment 


=.30402612 
1 .0704540 
-2.9928241 E -03 


.58181536 


.00000000 


.00000000 


.88362727 

.95734464 

9.6941615 E -02 


.37093279 


.00000000 


.00000000 


1.7380959 

.48288528 

.16077232 


1 .3682874 E -02 


.00000000 


.00000000 


1 .9250902 
-.17123713 
.16077192 


=.35172981 


.00000000 


.00000000 


1 .3702185 
=.76132099 
9.6940590 E -02 


=.57992482 


.00000000 


.00000000 


.29150092 

-1.0662339 

-2.9940458 E -03 


-.58189824 


.00000000 


.00000000 


“.89230506 

-.96599082 

-9.8792210 E -02 


-.36156328 


.00000000 


.00000000 


=1 .7309466 
=.49244525 
-.15593016 


-7.8057940 E -03 


.00000000 


.00000000 


=1.9120005 

.17398204 

=.15592911 


.34599696 


.00000000 


.00000000 


=1.3692757 

.77258517 

-9.8789370 E -02 


.57050831 


.00000000 


.00000000 


=.30399799 


.58181291 


.00000000 


.00000000 



1 .0704564 
“2.9903794 E -03 



- 118 - 



POTENTIAL SOLUTION - SINE WAVE 
RUN NUMBER NO. 6 



WL 


BL H 


DEPTH 


z 


RHO 


vise 


ALPHA 


SI 


S2 S3 


S4 


S 5 


s6 


S7 




5.0 


2.2917 .333 


2.0 


-.8333 


1.93 


9.66 E -06 


0.0 


.4725 


■=.1319 -.1529 


.0854 


-.1529 


0.0 


0.0 







FORMAT 






Vertical 


Lateral 


Vertical 


Lateral 


Force 


Force 


Moment 


Moment 


ETA 








CELERITY (Last Line Only) 






=.30275899 


.00000000 


.59160504 


. 00000000 


.00000000 








.90283537 
9. 7866442 E =02 


.00000000 


.37238541 


.00000000 


1.7635755 


. 00000000 


1 .0926558 E -02 


.00000000 


.15835106 

1.9506870 


.00000000 


-.35470589 


.00000000 


.15835068 
1 .3926989 
9.7865453 E -02 


.00000000 


-.58485213 


.00000000 


.30274462 


.00000000 


-.59160373 


.00000000 


-1.2217770 E -06 

=.90284827 


.00000000 


-.3 7238179 


.00000000 


=9. 7867431 E =02 
=1.7635821 


.00000000 


-1 .0922017 E -02 


.00000000 


=.15835144 

=1.9506846 


.00000000 


.35470961 


.00000000 


=.15835030 

=1.3926885 


.00000000 


.58485362 


.00000000 


=9. 7864463 E =02 
=.30273029 


,00000000 


.59160239 


.00000000 



2.44471! 95 E -0 6 
5.0599504 



- 119 - 



MODIFIED MORRISON SOLUTION - ST OKI AN WAVE 
RUN NUMBER 6 



WL 


BL 


H DEPTH Z 


RHO 


vise 


ALPHA 


SI 


S2 


S3 S4 S5 


S6 


S7 




5.0 


2.2917 


.333 2.0 -.8333 


1.93 


9.66E-06 


.00001 


.4726 


.1319 


-.1529 -.0854 -. 


1529 


0.0 0.0 





FORMAT 



Vertical 

Force 

RATIO 

ETA 

" O o O o o « o 

.243179116 
“.30979921 
-2.9928239E-03 
1; 4377068 
“3;7846560E-02 
9.6941615E-02 
2; 0782457 
“8.1 886751 E-03 
.16077232 
1.9229715 
1 .6656603E“02 
.16077191 
1; 0279684 
6; 56884 17E“02 
9.6940590E-02 
~. 25696082 
=.29272429 
“2;9940456E»O3 
“1.4369942 
«3;7649103E«02 
“9;8792210E“02 
-2; 0641 397 
“8.1578538E-03 
-.15593016 
-1; 9 1481 36 
1. 65421 20E-02 
-.15592911 
-1; 0371724 
6'.4862716E-Q2 
-9.8789370E“02 
-.24321012 
“.30975883 
”2. 990379 2E-03 



Lateral 

Force 

. CELERITY 
2.3681461E-05 

5.1401008 

1.6909574E“05 

5.1401008 

3.8038297E-06 

5.1401008 

-1.0724437E-05 

5. 1401008 
-2. 1262513E-05 

5.1401008 

-2.3774998E-05 

5.1401008 
-1.7 15929 IE-05 

5.1401008 

-3.8642422E-06 

5. 1401008 
1.0937246E-05 

5. 1401008 
2. 1454866E-05 

5. 1401008 
2.3681405E-05 

5.1401008 



Vertical 

Moment 

.62698556 

.64723345 

.40566327 

1.4 14485 IE-02 

-.37553658 

-.62728462 

-.63707673 

-.39913365 

-2.0227219E-02 

.36524532 

.62698920 



Lateral 

Moment 

2.1041582E-06 
-2.4772545E-06 
-6. 084839 1E-06 
-7.4428914E-06 
-6.0317703E-06 
-2.2876035E=06 
2.4221004E-06 
6. 234241 IE-06 
7.5904264E-06 
5.9735358E-06 
2. 104053 IE-06 



- 120 - 





POTENTIAL SOLUTION 
RUN NUMBER 


- ST OKI AN 
6 


WAVE 




WL BL 


H DEPTH Z 


RHO' 


vise 


ALPHA 


SI S2 


S3 S4 S5 


S6 


S7 




5.0 2.2917 


.333 2.0 =.8333 


1.93 


9 . 66E 


-06 0.0 


.4726 .1319 


=.1529 -.0854 

FORMAT 


-.1529 ( 


D.O 


0.0 


Vertical 


Vertical 


Lateral 




Lateral 


Force 

AXIAL FORCE 
ETA 


Moment 


Force 




Moment 


.29160968 
1.0667290 
=2.9928241 E =03 


.58226004 


.00000000 




.00000000 


1.3708319 

.76168365 

9.6941615 E -02 


.58024713 


.00000000 




.00000000 


1.9259741 

.17133102 

.16077232 


.35188697 


.00000000 




.00000000 


1.7389105 

=.48309554 

.16077192 


=1 .3752260 E -02 


.00000000 




.00000000 


.88405976 

“.95778034 

9.6940590 E -02 


=.37120049 


.00000000 




.00000000 


=.30413929 

=1.0709503 

-2.9940458 E -03 


=.58217668 


.00000000 




.00000000 


=1.3699091 

=.77294061 

=9.8792210 E =02 


-.57082253 


.00000000 




.00000000 


=1.9128844 

=.17406119 

=.15593016 


-.34614419 


.00000000 




.00000000 


=1.7317457 

.49267300 

=.15592911 


7.8805147 E -03 


.00000000 




.00000000 


=.89271518 

.96643675 

=9.8789370 E =02 


.36183327 


.00000000 




.00000000 


.29163802 

1.0667263 

=2.9903794 E =03 


.58226279 


.00000000 




.00000000 



- 121 - 



POTENTIAL SOLUTION - SINE WAVE 
RUN NUMBER 6 

WL BL H DEPTH Z RHO VISC 

SI S 2 S 3 S 4 S 5 S 6 S 7 

5.0 2.2917 .333 2.0 -.8333 1.93 9.66 E -06 

.4726 .1319 -.1529 -.0854 -.1529 0.0 0.0 



ALPHA 



0.0 



FORMAT 



Vertical Lateral 

Force Force 

ETA 

CELERITY (Last Line Only) 



.30285746 

.00000000 


.00000000 


1.3933225 

9.7866442 E -02 


.00000000 


1.9515833 

.15835106 


.00000000 


1.7644023 

.15835068 


.00000000 


.90327658 

9.78654536=02 


.00000000 


=. 30287 180 
-1.2217770 E -06 


.00000000 


-1.3933328 
=9. 786743 IE -02 


.00000000 


”1.9515858 

-.15835144 


.00000000 


”1.7643958 

-.15835030 


.00000000 


-.90326363 

-9.7864463 E -02 


.00000000 


.30288615 

2.4447195 E -06 

5.0599504 


.00000000 



Vertical 

Moment 


Lateral 

Moment 


59197134 


.00000000 


58517722 


. 00000000 


35486429 


.00000000 


.0995362 E -02 


.00000000 


37265512 


. 00000000 


59197267 


. 00000000 


58517575 


.00000000 


35486054 


.00000000 


.0999907 E -02 


.00000000 


37265877 


.00000000 


59197399 


.00000000 



-122 



EXPERIMENTAL RESULTS 







Run Number 7 






Model 


Angle of 
Attack 


Depth Below Wave 

Surface Height 


Wave 

Length 


Depth of 
Water 


2-7/8" 


cyl. 0° 


6" 4.0" 


5.0' 


2.0' 



Phase 

Angle 


ETA 

ft. 


Vertical 

Force 

lbs. 


Vertical 
Moment 
ft. lbs. 


Axial 

Force 

lbs. 


0 


-.045 


-0.180 


0.970 


0.564 


36 


.050 


0.565 


1.058 


0.256 


72 


.160 


0.860 


0.525 


0.145 


108 


.178 


0.940 


-0.125 


-0.0454 


144 


.082 


0.620 


-0.694 


-0.291 


l8o 


-.011 


-0.0 


-0.940 


-0.564 


216 


-.099 


-0.670 


-0.852 


-O .509 


252 


-.155 


-1.110 


-0.402 


-0.147 


288 


-.157 


-1.170 


0.206 


O.O 856 


524 


-.158 


-0.850 


0.820 


0.262 


560 


-.095 


-0.200 


1.040 


0.546 



- 125 - 



MODIFIED MORRISON SOLUTION - STOKIAN WAVE 
RUN NUMBER NO. 7 



WL BL 


H DEPTH 


Z RHO 


Vise ALPHA 


SI S2 


S3 S4 


S5 s6 


S7 


5.0 2.5 


.333 2.0 


-.50 1.93 


9.66E-06 .00001 


.23958 0.0 


0.0 0.0 


0.0 0.0 


0.0 




FORMAT 






Vertical 


Lateral 


Vertical 


Lateral 


Force 


Force 


Moment 


Moment 


RATIO 








ETA- CELERITY 






-.10045818 


1.1731923E-05 


.85684437 


-9.7913000E-09 


1.0000865 








-2.9928239E-03 


5.1401008 






.54235583 


9.4915410E-06 


. 74933124 


-5.491 1659 E-06 


=.16737838 








9.6941615E-02 


5.1401008 






.98350970 


3.6257727E-06 


.35258253 


-8.8720026E-06 


=4.1596958E=02 








.16077232 


5.1401008 






1.0653301 


-3.6248424E-06 


-.18454019 


-8 .872 1240E-06 


3.84Q3023E-02 








.16077191 


5.1401008 






.72390724 


-9. 4907862 E-06 


-.64958780 


-5. 49 1487 IE-06 


.12540162 








9.6940590E-Q2 


5.1401008 






.10046756 


-1 . 1731 45 1 E— 05 


-.85774119 


-1.0177070E-08 


.99999253 








-2.9940456E-03 


5.1401008 






-.54230327 


-9.4909970E-06 


-.73796280 


5.4850365E-06 


-.16749355 








-9.8792210E-02 


5.1401008 






=.98358730 


-3.6252022E-06 


-.34440498 


8.8882383E-06 


-4.1531486E-02 








-.15593016 


5.1401008 






=1 .0652823 


3.6253802E-06 


.17814737 


8.8883073E-06 


3.8348995E-02 








-.15592911 


5.1401008 






-.72394949 


9.4912768E-06 


.63738151 


5.4852326E-06 


.12546893 








-9.8789370E-02 


5.1401008 






-.10044247 


1.1731925E-05 


. 85684566 


-9.9280300E-09 


1.0002437 








-2.9903792E-03 


5.1401008 







- 124 - 



POTENTIAL SOLUTION - STOKIAN WAVE 



RUN NUMBER NO. 7 



WL 


BL 


H 


DEPTH 


z 


RHO 


vise 


ALPHA 


SI 


S2 


S3 


S4 


S5 


S6 


S7 




5.0 


2.5 


.333 


2.0 


O 

un 

• 

i 


1.93 


9.66 E -06 


0.0 


.23958 0.0 


0.0 


0.0 


0.0 


0.0 


0.0 





FORMAT 



Vertical 


Vertical 


Lateral 


Lateral 


Force 

AXIAL FORCE 
ETA 


Moment 


Force 


Moment 


8.1 300000 E -06 
.52845497 
-2.3928241 E -03 


.80327558 


. 00000000 


.00000000 


.59324146 

.42752920 

9.6941615 E -02 


.65538552 


.00000000 


.00000000 


.95987326 

.16330047 

.16077232 


.25163669 


.00000000 


.00000000 


.95987095 

-.16330415 

.16077192 


-.25164238 


.00000000 


.00000000 


.59323546 

=.42753148 

9.694090 E -02 


-.65538900 


.00000000 


.00000000 


6.9000000 E -07 

=.52845497 

=2.9940458 E “03 


-.80327549 


.00000000 


.00000000 


=.59324246 

=.42752572 

=9.8792210 E “02 


-.64433682 


.00000000 


.00000000 


=.95988866 

=.16329602 

“.15593016 


-.24480489 


.00000000 


,00000000 


=.95988179 

.16330710 

-.15592911 


.24482150 


.00000000 


.00000000 


“.59322446 

.42753258 

“9.8789370 E =02 


.64434728 


.00000000 


.00000000 


2.2940000 E -05 

.52845498 

“2.9903794 E -03 


.80327575 


.00000000 


.00000000 



- 125 - 



POTENTIAL SOLUTION - SINE WAVE 



RUN NUMBER NO. 7 



WL 


BL 


H 


DEPTH 


z 


RHO 


vise 


ALPHA 


SI 


S2 


S3 


S4 


S5 


S6 


S7 




5.0 


2.5 


.333 


2.0 


50 


1.93 


9.66 E -06 


0.0 


.23958 0.0 


0.0 


0.0 


0.0 


0.0 


0.0 







FORMAT 






Vertical 


Lateral 


Vertical 


Lateral 


Force 


Force 


Moment 


Moment 


ETA 








CELERITY (Last Line Only) 






. 00000000 
.00000000 


.00000000 


.79146542 


.00000000 


.58451686 

9.7866442 E -02 


.00000000 


.64030830 


.00000000 


.94576717 


.00000000 


.24457406 


.00000000 


.15835106 

.94576490 


.00000000 


-.24457957 


.00000000 


.15835068 

.58451098 


.00000000 


-.64031174 


.00000000 


9. 7865453 E “02 
-7.29 GQ 000 E -06 
-1.2217770 E -06 


.00000000 


-.79146543 


.00000000 


=>.58452278 


.00000000 


-.64030489 


.00000000 


“9. 7867431 E -02 
». 94576943 
-. 15835144 


.00000000 


-.24456853 


.00000000 


94576266 

=*.15835030 


.00000000 


.24458509 


.00000000 


=*.58450507 

=9.7864463 E -02 


.00000000 


.64031514 


.00000000 


1 .4590000 E -05 
2.4447195 E -06 
5.0599504 


.00000000 


.79146544 


.00000000 



-126 



EXPERIMENTAL RESULTS 



Run Number 8 



Model 


Angle of 
Attack 


Depth Below 
Surface 


Wave 

Height 


Wave 

Length 


Depth of 
Water 


2-7/8" 


cyl. 0° 


10" 


4.0" 


5.0' 


2.0' 



Phase 

Angle 


ETA 

ft. 


Vertical 

Force 

lbs. 


Vertical 
Moment 
ft. lbs. 


Axial 
Force 
lbs . 


0 


-.045 


-0.115 


0.705 


0.575 


56 


.058 


0.210 


0.621 


0.518 


72 


.170 


0.450 


0.562 


0.147 


108 


.170 


0.460 


-0.0572 


-0.0772 


144 


.720 


0.510 


-0.424 


-0.262 


180 


0.000 


-0.0 


-0.500 


-0.564 


216 


-.105 


-0.54 


-0.460 


-0.241 


252 


-.155 


-0.570 


-0.225 


-0.147 


288 


-.157 


-0.650 


0.1555 


0.059 


524 


-.158 


-0.490 


0.542 


0.262 


560 


-.068 


-0 . 200 


0.697 


0.550 



- 127 - 



MODIFIED MORRISON SOLUTION - STOKIAN WAVE 
RUN NUMBER NO. 8 



WL 


BL 


H 


DEPTH 


Z 


RHO 


vise 


ALPHA 


SI 


S2 


S3 


S4 


S5 


S6 


S7 




5.0 


2.5 


.333 


2.0 


-.8333 


1.93 


9.66E-06 


.00001 


.23958 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 





FORMAT 



Vertical 


Lateral 


Vertical 


Lateral 


Force 

RATIO 


Force 


Moment 


Moment 


ETA-CELERITY 






-3.8186266E-02 


7.9453148E-06 


.55121012 


-4.6406470E-08 


1 .0000982 
-2.9928239E-03 


5.1401008 






.37272200 

=9.24l3739E-02 


6.4280829E-06 


.46759073 


-3-731 1588E-06 


9.6941615E-02 


5.1401008 






.64331278 

=2.4080146E°02 


2.4555025E=06 


.20377983 


-5.9763409E-06 


.16077232 


5.1401008 






.67429364 

2.2974254E=02 


-2.4550004E-06 


-.14013014 


-5.9764138E-06 


.16077191 


5.1401008 






.44160723 

7.7998510E-02 


-6.4276956E-06 


-.42960376 


-3.7313612E-06 


9.6940590E“02 


5.1401008 






3.8188677E-02 

1.0000345 


-7.9451227E-06 


-.55143539 


-4.6662200E-08 


-2.994G456E-03 


5.1401008 






=.37271072 
=9.245171 IE-02 


=6.4277262E=06 


-.46256107 


3.7024194E-06 


=9.8792210E=02 


5.1401008 






-.64333724 

-2.4053078E-02 


-2.4551783E=06 


-.20037928 


6.0516824E-06 


-.15593016 


5.1401008 






=.67428197 


2.455 17 18E-06 


.13718487 


6.0517352E-06 


2.2950594E=02 

-.15592911 


5.1401008 






=.44161413 


6.4278051 E-06 


.42436576 


3.7025515E-06 


7.8027478E-02 

-9.8789370E=02 


5.1401008 






=3.8176093E=02 


7.9453158E-06 


.55121060 


=4.64991 30E-08 


1.0003646 

=2.99037925=03 


5.1401008 







-128- 



POTENTIAL SOLUTION - STOKIAN WAVE 



RUN NUMBER NO. 8 



WL 


BL 


H 


DEPTH 


z 


RHO 


vise 


ALPHA 


SI 


S2 


S 3 


S 4 


S 5 


S 6 


S 7 




5.0 


2.5 


.333 


2.0 


-.8333 


1.93 


9.66 E -06 


0.0 


.23958 0.0 


0.0 


0.0 


0.0 


0.0 


0.0 





FORMAT 



Vertical 


Vertical 


Lateral 


Lateral 


Force 

AXIAL FORCE 
ETA 


Moment 


Force 


Moment 


3.5100000 E -06 
.35789295 
-2. 9928241 E -03 


.51205618 


.00000000 


.00000000 


.37816729 

.28954396 

9.6941615 E -02 


.41664592 


. 00000000 


.00000000 


.61188234 

.11059572 

.1607/232 


.15970615 


.00000000 


.00000000 


.61188085 

=.11059823 

.160/7192 


-.15970975 


.00000000 


.00000000 


.37816347 
-. 28954552 
9. 6940590 E =02 


-.41664813 


.00000000 


.00000000 


=1 . 2100000 E -06 
-.35789294 
=2. 9 9 4045 8 E -0 3 


-.51205618 


.00000000 


. 00000000 


=.37816898 
-.28953685 
-9 . 8/92210 E -02 


-.41187533 


.00000000 


.00000000 


-.61188945 

-.11058978 

-.15593016 


-.15675560 


.00000000 


. 00000000 


-.61188510 

.11059729 

-.15592911 


.15676623 


.00000000 


. 00000000 


-.37815752 

.28954148 

-9.8 J 89370 E -02 


.41188198 


.00000000 


. 00000000 


1 .2930000 E -05 

.35789297 

-2.9903794 E -03 


.51205628 


.00000000 


. 00000000 



- 129 - 



POTENTIAL SOLUTION - SINE WAVE 
RUN NUMBER NO. 8 



WL 


BL 


H 


DEPTH 


z 


RHO 


vise 


ALPHA 


SI 


S2 


S 3 


S 4 


S 5 


S 6 


S 7 




5.0 


2.5 


.333 


2.0 


-.8333 


1.93 


9.66 E -06 


0.0 


.23958 0.0 


0.0 


0.0 


0.0 


0.0 


0.0 







FORMAT 






Vertical 


Lateral 


Vertical 


Lateral 


Force 


Force 


Moment 


Moment 


ETA 








CELERITY (Last Line Only) 






.00000000 

.00000000 


.00000000 


.52063444 


.00000000 


.38450145 
9. 7866442 E =02 


.00000000 


.42120165 


.00000000 


.62213577 

.15835106 


.00000000 


.16088344 


.00000000 


.62213428 

.15835068 


. 00000000 


-.16088707 


. 00000000 


.38449757 
9. 7 8 6545 3 E -0 2 


.00000000 


-.42120390 


.00000000 


-4.7900000 E -06 

-1.2217770 E -06 


.00000000 


-.52063443 


.00000000 


=.38450534 
-9.7867431 E -02 


.00000000 


-.42119942 


.00000000 


=.62213725 

=.15835144 


.00000000 


-.16087980 


.00000000 


=.62213279 

=.15835030 


.00000000 


.16089070 


. 00000000 


=.38449369 
=9. 7864463 E =02 


.00000000 


.42120614 


.00000000 


9.6000000 E -06 

2.4447195 E -06 

5.0599504 


.00000000 


.52063443 


.00000000 



-I5O- 



EXPERIMENTAL RESULTS 



Run Number 9 



Model 


Angle of 
Attack 


Depth Below 
Surface 


Wave 

Height 


Wave 

Length 


Depth of 
Water 


2 - 7 / 8 " 


cyl. 0 ° 


14" 


4.0” 


5.0' 


2 . 0 ' 



Phase 

Angle 


ETA 

ft. 


Vertical 

Force 

lbs. 


Vertical 
Moment 
ft. lbs. 


Axial 
Force 
lbs . 


0 


-.044 


- 0.12 


0.497 


0.246 


36 


.050 


0.220 


0.434 


0.205 


72 


.171 


0.450 


0.211 


0.0955 


108 


.170 


0,475 


-0.0764 


-0.0454 


144 


.075 


0.305 


-0.348 


- 0.200 


180 


-.008 


- 0.0 


-0.364 


-0.291 


216 


-.091 


- 0.300 


- 0.322 


- 0.232 


252 


-.150 


- 0.500 


-0.159 


-0.109 


288 


-.150 


-0.525 


0.120 


0.0227 


324 


-.125 


-0.350 


0.400 


0.1545 


360 


-.044 


- 0.090 


0.440 


0.232 



- 131 - 



MODIFIED MORRISON SOLUTION - STOKE AN WAVE 



RUN NUMBER NO. 9 



WL 


BL 


H 


DEPTH 


z 


RHO 


vise 


ALPHA 


SI 


S2 


S3 


S4 


S5 


S6 


S7 




5.0 


2.5 


.3333 


2.0 


-1.1667 


1.93 


9.66E-06 


.00001 


.23958 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 





FORMAT 

Vertical Lateral 

Force Force 

RATIO 

ETA-CELERITY 



-1. 36297 10E-02 
1 .00011198 
-2.9982189E-03 
.23819585 
-5.1261350E-02 
9.70281 20E-02 
.39976126 
=1 .35 12446E-02 
.16091934 
.41056391 
1 .3157189E-02 
.16091893 
.26261383 
4.6495156E-02 
9.7027090E-O2 
1 .3629824E-=02 
1 .0001109 
=2.999441 7 E-03 
=.23819433 
=5. 1274583 E=02 
=9.88820505=02 
-.39976762 
-1. 3505355 E-02 
=.15606845 
=.41056307 
1.3151108E-02 
-.15606740 
=.26261347 
4.6507267E-02 
=9.8879205 E =02 
-1.3623464E-02 
1.0005784 
-2.9957720E-03 



5.5779876E-06 

5.1403027 

4.5128549E-06 

5.1403027 
1 .7238892E-06 

5.1403027 
-1.723581 IE-06 

5.1403027 

-4.5126304E-06 

5.1403027 

-5.5779025E-06 

5.1403027 

-4.5125741E-06 

5.1403027 

-1.7236630E-06 

5.1403027 
1 .7236025E-06 

5.1403027 
4.51 25797 E-06 

5.1403027 

5.5779886E-06 

5.1403027 



Vertical 

Moment 

.33903409 

.28218506 

.11674675 

-9.4175370E-02 

-.26864597 

-.33908964 

-.27999279 

-.11532011 

9.2862308E-02 

.26640541 

.33903427 



Lateral 

Moment 

-6.2296580E-08 
-2. 6286462E-06 
-4.1716517E-06 
-4. 17 17003 E-06 
-2.6287874E-06 
-6.2479620E-08 
2. 590099 IE-06 
4.2726280E-06 
4.2726685E-06 
2.5901974E-06 
-6. 2361610E-08 



- 132 - 



POTENTIAL SOLUTION - STOKIAN WAVE 



RUN NUMBER NO. 9 



WL BL 


H DEPTH Z 


RHO 


vise 


SI S2 


S3 S4 S5 


s6 


S7 


5.0 2.5 


.3333 2.0 -1. 


1667 1.93 


9.66E 


.23958 0.0 1 


D.O 0.0 0.0 

FORMAT 


0.0 


0.0 


Vertical 


Vertical 


Lateral 




Force 

AXIAL FORCE 
ETA 


Moment 


Force 




1 .4920000E-06 
.25125886 
=2.9982191 E-03 


.31224070 


.00000000 


•* 


.23059796 

.20327622 

9.7028120E-02 


.25362835 


.00000000 




.37311299 

7.7644880E-02 

.16091934 


9.71 17562E-02 


.00000000 




.37311209 

-7.7646626E-02 

.16091893 


=9.71 19740E-02 


.00000000 




.23059564 

=.20327731 

9.7027090E-02 


-.25362969 


.00000000 




=1 .4010000E-06 
=.25125886 
=2. 999441 9E=03 


-.31224067 


.00000000 




=*.23059939 

=.20326765 

=9.8882050E-02 


-.25158583 


.00000000 




=*.37311633 

=7.7638500E=02 

=.15606845 


-9.5853898E-02 


.00000000 




=.37311365 
7. 764377 IE-02 
-.15606740 


9. 586041 4E-02 


.00000000 




-.23059241 

.20327091 

-9.8879205E-02 


.25158989 


.00000000 




7 . 2440000E-06 


.31224073 


.00000000 





. 25125886 
-2.9957722E-03 



ALPHA 
06 0.0 

Lateral 

Moment 

.00000000 

.00000000 

.00000000 

.00000000 

.00000000 

.00000000 

.00000000 

.00000000 

.00000000 

.00000000 

.00000000 



-153- 



POTENTIAL SOLUTION - SINE WAVE 
RUN NUMBER NO. 9 



WL 


BL 


H 


DEPTH 


z 


RHO 


vise 


ALPHA 


SI 


S2 


S 3 


S 4 


S 5 


S 6 


S 7 




5.0 


2.5 


.333 


2.0 


-1.1667 


1.93 


9.66 E -06 


0.0 


.23958 0.0 


0.0 


0.0 


0.0 


0.0 


0.0 







FORMAT 






Vertical 


Lateral 


Vertical 


Lateral 


Force 


Force 


Moment 


Moment 


ETA 








CELERITY (Last Line Only) 






1 .0000000 E -09 


.00000000 


.34243582 


.00000000 


.00000000 

.25289736 


.00000000 


.2770361 1 


.00000000 


9.7866442 E - Q 2 

.40919610 

.15835106 


.00000000 


.10581753 


.00000000 


.40919511 

.15835068 


.00000000 


-.10581992 


.00000000 


.25289481 

9.7865453 E -02 


.00000000 


-.27703760 


.00000000 


-3. 1520000 E -06 
-1 .2217770 E -06 


.00000000 


-.34243583 


.00000000 


-.25289992 
-9. 786743 IE -02 


.00000000 


-.27703464 


.00000000 


-.40919707 

-.15835144 


.00000000 


-. 10581514 


.00000000 


-.40919413 

-.15835030 


.00000000 


.10582231 


.00000000 


-.25289224 

-9.7864463 E -02 


.00000000 


.27703905 


.00000000 


6.3210000 E -06 

2.4447195 E -06 

5.0599504 


.00000000 


.34243581 


.00000000 



- 154 - 



EXPERIMENTAL RESULTS 



Run Number 10 



Model 


Angle of 
Attack 


Depth Below 
Surface 


Wave 

Height 


Wave 

Length 


Depth of 
Water 


1" cyl. 


0° 


10" 


4" 


5.0' 


2.0' 



Phase 

Angle 


ETA 

ft. 


Vertical 

Force 

lbs. 


Vertical 
Moment 
ft. lbs. 


Axial 

Force 

lbs. 


0 


-.045 


-0.020 


0.050 


0.045 


36 


.048 


0.025 


0.073 


0.051 


72 


.160 


0.057 


0.001 


0.007 


108 


.168 


0.060 


-0.018 


-0.015 


144 


.075 


0.054 


-0.055 


-O.O 56 


180 


-.006 


0.025 


- 0.070 


-0.054 


216 


-.124 


-0.010 


- 0.080 


-0.045 


252 


-.160 


-0.054 


-0.040 


-0.025 


288 


-.165 


-0.055 


40.019 


0.007 


524 


-.155 


-0.040 


40.046 


0.051 


560 


-.045 


-0.009 


0.050 


0.044 



- 135 - 



MODIFIED MORRISON SOLUTION - STOKIAN WAVE 



RUN NUMBER NO. 10 



WL 


BL 


H 


DEPTH 


z 


RHO 


vise 


ALPHA 


SI 


S2 


S3 


S4 


S5 


S6 


ST 




5.0 


2.5 


.333 


2.0 


-.8333 


1.93 


9 . 66 E -06 


.00001 


.0833 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 





FORMAT 



Vertical 

Force 

RATIO 

ETA-CELERITY 

H.2810525E-02 
1 .0000341 
-2.9928239E-G3 
3. 5845034 E-02 
•=.317511491 
9.6941615E-02 
7.1318274E-02 
-7.1439488E-02 
.16077232 
8.15O80?OE-O2 
6.2509647E-02 
.16077191 
5.8607244E-02 
.19419706 
9.6940590E-02 
1. 281082 IE-02 
1.0000119 
-2.994Q456E-03 
“3.5841463E-02 
=>.31765204 
“9.8792210E-02 
=7. 1325904E-02 
-7.1337350E-02 
*.15593016 
“8. 1502066E-Q2 
6.2434105E=02 
-.15592911 
-5. 86 10359 E-02 
.19425258 
=9.8789370E“O2 
=1.2809363E°02 
1.0001259 
°2. 9903792 E-03 



Lateral 

Force 



9.6059051 E-07 

5.1401008 
7. 77 29 245 E-07 

5.1401008 

2.9712008E-07 

5.1401008 

-2.9650780E-07 

5.1401008 

-7.7683522E-07 

5.1401008 

-9.6039534E-07 

5.1401008 

-7.7706175E-07 

5.1401008 

-2.9690855E-07 

5.1401008 

2.9670088E-07 

5.1401008 

7.7703604E-07 

5.1401008 
9. 6059066 E-07 

5.1401008 



Vertical 

Moment 

6.3910638E-02 
5.8332292E-02 
3.0427627E-02 
-9.461 7670E-03 
-4.5731380E-02 
-6.3982845E-02 
-5.7770185E-02 
-2.9991294E-02 
9.1620840E-03 
4.5102641E-02 
6.3910796E-02 



Lateral 

Moment 

-5.45661 80E°09 
-4. 5094491 E-07 
-7. 224405 6E-07 
-7.2252108E-07 
=4.51 19486E-0? 
-5.7943800E-09 
4.4745134E-07 
7.3152894E-07 
7.3164625E-07 
4.4773135E-07 
-5.4678250E-09 



-136- 



POTENTIAL SOLUTION - STOKXAN WAVE 



RUN NUMBER NO. 10 



WL BL 


H DEPTH Z 


RHO 


vise 


ALPHA 


SI S2 


S3 S4 S5 


S6 


S7 




5.0 2.5 


.333 2.0 -.8333 1.93 


9 . 66 E-O 6 0.0 


.0833 0.0 0.0 0.0 0.0 

FORMAT 


0.0 


0.0 




Vert iced 


Vertical 


Lateral 




Lateral 


Force 

AXIAL FORCE 
ETA 


Moment 


Force 




Moment 


4.2300000E-Q7 
4.3265530E-02 
“2.9928241 E-03 


6.1902265E-02 


.00000000 




. 00000000 


4.571 6489E-02 
3.5002849E-02 
9.6941615E-02 


5.0368153E-02 


.00000000 




.00000000 


7 • 3970203E-02 
1 .3369871 E-02 
.16077232 


1 .9306810E-02 


.00000000 




.00000000 


7. 3970023 E-02 
-1.3370174E-02 
.16077192 


-1 .9307244E-02 


.00000000 




. 00000000 


4.5716027E-02 

-3.5003036E-02 

9.6940590E°02 


-5.0368419E-02 


.00000000 




.00000000 


“1 .4700000E-07 
“4.3265530E-02 
“2.9940458E-03 


-6. 1902258E-02 


.00000000 




.00000000 


“4. 5716691 E-02 
-3.5001988E-02 
“9. 8792210E-02 


-4. 979 1437 E-02 


.00000000 




.00000000 


“7.3971062E-02 
-1.3369 153 E-02 
-.15593016 


-1 . 89501 18E-02 


.00000000 




.00000000 


-7.3970534E-02 

1.3370060E-02 

-.15592911 


1 .8951404E-02 


.00000000 




.00000000 


-4.5715306E-02 

3.5002548E-02 

-9.8789370E-02 


4.9792242E-02 


.00000000 




.00000000 


1 .5630000E-06 
4.3265531 E-02 
“2.9903794E-03 


6.1902273E-02 


.00000000 




.00000000 



- 157 - 



POTENTIAL SOLUTION - SINE WAVE 



RUN NUMBER NO. 10 



VL BL 


H DEPTH 


Z RHO Vise 


ALPHA 


SI S2 


S3 S4 


S5 S6 S7 




5.0 2.5 


333 2.0 


.8333 1.93 9.66E- 


-06 0.0 


.0833 0.0 0 


1.0 0.0 0 

FORMAT 


.0 0.0 0.0 




Vertical 


Lateral 


Vertical 


Lateral 


Force Force 

ETA 

CELERITY (Last Line Only) 


Moment 


Moment 


.00000000 

.00000000 


.00000000 


6.2939282E-02 


.00000000 


4.6482220E-02 
9. 7866442 E-02 


.00000000 


5.0918893E-02 


.00000000 


7.5209735E-02 

.15835106 


.00000000 


1 .94491 3 IE-02 


.00000000 


7.5209557E-02 

.15835068 


.00000000 


-1.944957 IE-02 


.00000000 


4.6481753E-02 
9. 7865453 E=02 


.00000000 


-5.0919165E-02 


.00000000 


=5 . 8000000E-07 
=1 .2217770E-06 


.00000000 


-6.2939279E-02 


.00000000 


-4.6482692E-02 
=9.7867431 E-02 


.00000000 


-5.0918622E-02 


.00000000 


-7.5209916E-Q2 

=.15835144 


.00000000 


=1 .9448692E-02 


.00000000 


=7.5209377E=02 

=.15835030 


.00000000 


1 .9450009E-02 


.00000000 


=4.6481281 E=Q2 
=9. 7864463 E=02 


.00000000 


5.0919437E-02 


.00000000 


1.1600000E-06 
2. 44471 95 E-0 6 
5.0599504 


.00000000 


6. 2939281 E-02 


.00000000 



-I38- 



EXPERIMENTAL RESULTS 







Run Number 11 








Model 


Angle of 
Attack 


Depth Below 
Surface 


Wave 

Height 


Wave 

Length 


Depth of 
Water 


4" cyl. 


45° 


10" 


.333' 


5.0' 


2.0' 



Phase 

Angle 


ETA 

ft. 


Vertical 

Force 

lbs. 


Lateral 
Moment 
ft. lbs. 


0 


- .066 


-0.40 


O .165 


36 


.050 


0.45 


- 0.150 


72 


.155 


1.15 


-0.555 


108 


.176 


1.375 


- 0.508 


144 


.086 


0.900 


-0.555 


180 


.000 


0.170 


- 0.152 


216 


-.081 


-0.895 


0.220 


252 


-.150 


-1.425 


0.412 


288 


-.140 


- 1.50 


0.510 


524 


-.120 


-1.13 


O .565 


560 


00 

Lf\ 

O 

0 

1 


- 0.565 


0.158 



- 139 - 



MODIFIED MORRISON SOLUTION - STOKIAN WAVE 
RUN NUMBER NO. 11 



WL BL 


H DEPTH 


Z RHO 


Vise ALPHA 


SI S2 


S3 S4 


S5 s6 


S7 


5.0 2.5 


.333 2.0 


-.8333 1.93 


9.66E-06 .7854 


.3333 0.0 


0.0 0.0 


0.0 0.0 


0.0 




FORMAT 






Vertical 


Lateral 


Vertical 


Lateral 


Force 


Force 


Moment 


Moment 


RATIO 








ETA-CELERITY 






=8.071 1259E-02 


1.3556066 


. 86440597 


4.3132300E-03 


.94162735 








-2.9928239E-03 


5.1401008 






.94046140 


1.1129530 


.73598710 


-.38441282 


“6.6442780E=02 








9.6941615E=02 


5.1401008 






1.6036812 


.45788096 


.30424316 


-.63039264 


=1.5769Q66E=02 








.16077232 


5.1401008 






1.6542549 


-.36568620 


-.23855109 


-.64874333 


1.5287344E-02 








.16077191 


5.1401008 






1.0654249 


-1.0497193 


-.67667498 


-.41978674 


5.8649902E-02 








9.6940590EHD2 


5.1401008 






7.1276057E-02 


-1,3306928 


-.86495544 


-2. 1183350E-02 


1.0662756 








-2.9940456E-03 


5.1401008 






“.94367970 


-1.0945443 


-.72256447 


.39012173 


»6.5901916E=O2 








=9.8792210E=02 


5.1401008 






-1.5963416 


-.42636886 


-.29541750 


.65065809 


“1.5666322E-02 








“. 15593016 


5.1401008 






-1.6463493 


.41029211 


.23059435 


.65578251 


1.5191465E-02 








-.15592911 


5.1401008 






-1.0680311 


1.0893145 


.66294345 


.40364841 


5. 82 29 859 E-0 2 








-9.8789370E-02 


5.1401008 






-8.0685962E-02 


1.3556068 


, 86440703 


4. 3033800E-03 


.94192038 








“2.9903792E-03 


5.1401008 







140 - 



POTENTIAL SOLUTION - STOKIAN WAVE 



RUN NUMBER NO. 11 



WL 


BL 


H 


DEPTH 


z 


RHO 


vise 


SI 


S2 


S3 


S4 


S5 


S6 


S7 


5.0 


2.5 


.333 


2.0 


-.8333 


1.93 


9.66E 


.3333 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 



FORMAT 



Vertical Vertical Lateral 

Force Moment Force 

AXIAL FORCE 
ETA 



-4.3504400E-03 


.79844474 


1.2413204 


.62065791 
-2.9928241 E-03 
.92612099 


.65222329 


.99935320 


.49967476 
9.6941615E-02 
1 .5041883 


.25060494 


.38056036 


.19027949 

.16077232 

1.5041847 


-.25061060 


-.38056897 


“.19028379 

.16077192 

.92611150 


-.65222668 


-.99935861 


49967745 

9.6940590E-02 

“4.3620200E-03 


-.79844462 


-1.2413204 


“.62065795 

“2.9940458E-03 

“.92881906 


-.63968263 


-1.0091380 


-. 50456710 
“9.8792210E-02 
“1. 4971527 


-.24285113 


-.38660216 


-.19330037 

-.15593016 

-1.4971421 


.24286761 


.38662839 


.19331347 

-.15592911 

-.92879100 


.63969308 


1 .0091539 


.50457511 
~9. 8789370E-02 
-4.3272400E-03 


.79844491 


1 .2413202 



.62065786 

-2.9903794E-03 



ALPHA 
*06 .7854 

Lateral 

Moment 

7.7979500E-03 

-.37161692 

-.59108048 

-.59107913 

-.37161330 

-7.793310QE-03 

.36680138 

.60369940 

.60369490 

.36678987 

-7.8071900E-03 



-l4l- 



POTENTIAL SOLUTION - SINE WAVE 
RUN NUMBER NO. 11 



WL 


BL 


H 


DEPTH 


Z 


RHO 


vise 


ALPHA 


SI 


S2 


S3 


S4 


S5 


s6 


S7 




5.0 


2.5 


.333 


2.0 


-.8333 


1.93 


9.66 E -06 


.7854 


.3333 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 












FORMAT 








Vertical 




Lateral 




Vertical 




Lateral 


Force 






Force 




Moment 




Moment 



ETA 

CELERITY (Last Line Only) 



.00000000 

.00000000 


1.1344337 


.81182070 


.00000000 


.94300269 
9. 7866442 E -02 


.91777527 


.65677606 


-.33741583 


1. 5 258088 
.15835106 


.35055616 


.25086413 


-.54594971 


1.5258050 

.15835068 


-.35056408 


-.25086978 


-.54594840 


.94299310 

9.7865453 E -02 


“.91778015 


-.65677954 


-.33741244 


“1.1760 G 00 E -05 

-1.2217770 E -06 


-1.1344337 


-.81182068 


4.2100000 E 


=.94301224 
=9. 7867431 E -02 


“.91777036 


-.65677254 


.33741926 


-1.5258124 

“.15835144 


-.35054825 


-.25085846 


.54595099 


“1.5258013 

-.15835030 


.35057200 


.25087545 


.54594710 


“.94298360 

“9.7864463 E -02 


.91778507 


.65678305 


.33740902 


2.3540000 E -05 

2.4447195 E -06 

5.0599504 


1.1344338 


.81182067 


-8.4300000 E 



- 142 - 



EXPERIMENTAL RESULTS 



Run Number 12 



Model Angle of Depth Below Wave Wave Depth of 

Attack Surface Height Length Water 

Sub 45° 10" 4" 5-0' 2.0' 



Phase 

Angle 


ETA 

ft. 


Vertical 
Force 
lbs . 


Lateral 
Moment 
ft . lbs 


0 


0.0 


-0.570 


0.575 


56 


.096 


0.800 


0.275 


72 


.168 


1.785 


0.040 


108 


.160 


1.740 


-0.412 


144 


.078 


0.880 


-O.658 


180 


-.018 


-0.220 


-0.228 


216 


-.104 


-1.250 


-0.012 


252 


-.154 


-1.80 


0.518 


288 


-.162 


-1.840 


0.158 


524 


-.095 


-1 . 250 


0.525 


560 


0.0 


-0.160 


0.550 



- 145 - 



MODIFIED MORRISON SOLUTION - STOKIAN WAVE 



RUN NUMBER NO. 12 



Z RHO 

S5 S6 

-.8333 1.93 

.1529 0.0 



Vise ALPHA 

S7 

9.66E-06 .7854 

0.0 



WL BL 

SI S2 

5.0 2.2917 



.4725 =.1319 



Vertical 
Force 
RATIO 
ETA-CELERITY 



=.34686607 

.26385562 

=2.9928239E=03 

1.0711499 

-6.8259710E-02 

9;6941615E“02 

2.0856479 

=1.2872018E=02 

.16077232 

2.2835412 

7o386Q835E=03 

,16077191 

1.6056247 

3.7786663E-02 

9.6940590E=02 

.32566528 

.28144983 

=2.9940456E=03 

=1.0831157 

=6.6980028E=02 

=9;8792210E=02 

=2.0718152 

=1.2742726E=02 

=.15593016 

=2.2630461 

7.3288308E-03 

=.15592911 

=1.6068091 

3o7386962E=02 

=9.878937QE=02 

=.34683275 

.26388177 

=2.9903792E=03 



H DEPTH 
S3 S4 
.333 2.0 



-.1529 .0854 



FORMAT 

Lateral 

Force 



1.8130627 

5.1401008 

1.5880276 

5.1401008 

.78059646 

5. 1401008 
-.31500120 

5. 1401008 

- 1.2978882 

5.1401008 

=1.7883219 

5.1401008 

=1.5849077 

5. 1401008 
=.75885390 

5. 1401008 
.36880574 

5.1401008 

1.3537978 

5.1401008 

1.8130655 

5. 1401008 



Vertical 

Moment 

.48054794 
.27208392 
=6. 1388744E-02 
-.36519791 
-.51482938 
-.47957981 
=.26238120 
6.6296023E-02 
.35855059 
.50590800 
.48054515 



Lateral 

Moment 

-.19023335 
-.36521406 
-.40573238 
-.30128432 
-8.1956835E-02 
.17630616 
.36837013 
.41677843 
.30262819 
6.6600537E-02 
=. 19023872 
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POTENTIAL SOLUTION - STOKIAN WAVE 



RUN NUMBER NO. 12 



WL 


BL H DEPTH 


z 


RHO 


vise 


ALPHA 


SI 


S2 S3 


S4 


S5 


s6 


S7 




5.0 


2.2917 .333 


2.0 


-.8333 


1.93 


9.66E-06 


.7854 


.4725 


=.1319 =.1529 


.0854 


-.1529 


0.0 


0.0 







FORMAT 






Vertical 


Vertical 


Lateral 


Lateral 


Force 

AXIAL FORCE 
ETA 


Moment 


Force 


Moment 


-.23387342 
.83097342 
=2. 992 824 IE-03 


.43835437 


1.6619528 


-.17910207 


1.0532023 

.71788710 

9.6941615E-02 


.22877252 


1.4357795 


-.34391059 


1.9432085 
.33567919 
. 16077232 


=7. 2679808E-02 


.67136089 


-.37642837 


2.0844044 

=.17023849 

.16077192 


-.34850318 


-.34047823 


-.26992435 


1.4201464 

-.61343203 

9.6940590E-02 


-.48804073 


-1.2268685 


-6.41887 24E-02 


.21427129 

-.82824556 

-2.9940458E-03 


-.43707164 


-1 . 6564970 


.16843436 


-1.0636552 
-.72806019 
=9. 8792210E-02 


-.21979586 


-1.4561257 


.34205610 


-1.9300801 

-.34469126 

-.15593016 


7.6948100E-02 


-.68938506 


.38595126 


-2.0658441 

. 1 7484745 
-.15592911 


.34216804 


.34969621 


.27766383 


-1.4218009 

.62529852 

-9.8789370E-02 


.47985970 


1.2506016 


5.9447037E-02 


-.23384252 

.83097455 

-2.9903794E-03 


.43835125 


1.6619550 


-. 17910714 
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POTENTIAL SOLUTION - WAVE SINE 



RUN NUMBER NO. 12 



WL BL 


H DEPTH 


Z RHO VISC 


ALPHA 


SI S2 


S3 S4 


S5 s6 S7 




5.0 2.2917 


.333 2.0 


-.8333 1.93 9.66E 


-06 .7854 


.4725 “.1319 


=.1529 .0854 
FORMAT 


-.1529 0.0 0.0 




Vertical 


Lateral 


Vertical 


Lateral 


*°rca Force 

ETA 

CELERITY (Last Line Only) 


Moment 


Moment 


= . 22783398 
. 00000000 


1.5163531 


.44504662 


=. 15880438 


1.0761532 

9.7866442E-02 


1.3214484 


.22804316 


-.31344939 


1.9690847 

.15835106 


.62179317 


=7. 6065490E-02 


=.34836683 


2.1098892 

.15835068 


-.31536719 


-.35111957 


-.25021939 


1 .4447842 
9.7865453E-02 


=1.1320673 


-.49205727 


-5.6496230E-02 


.22781824 

-1.2217770E-06 


-1.5163542 


-.44504497 


. 15880668 


=1.0761670 
=9. 786743 1E=02 


=1.3214427 


-.22803989 


.31345058 


=1,9690912 

=.15835144 


-.62178291 


7.6069086E-02 


.34836641 


=2. 1098860 
=.15835030 


.31537815 


.35112216 


.25021757 


=1.4447724 

=9.7864463E=02 


1.1320748 


.49205786 


5.6493670E-02 


=.22780250 

2.4447195E=06 


1.5163555 


.44504332 


-.15880901 



5.0599504 
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MODIFIED MORRISON SOLUTION - STOKIAN WAVE 
RUN NUMBER 12 



WL 


BL 


H DEPTH 


Z 


Rho VIS 


ALPHA 


SI 


S2 


S3 S4 


S5 


S6 S7 




5.0 


2.2917 


.333 2.0 


.8333 


1,93 9.66E-06 


.7854 


.4726 


.1319 


-.1529 -.0854 


-.1529 0.0 0.0 





FORMAT 



Vertical 


Lateral 


Vertical 


Lateral 


Force 


Force 


Moment 


Moment 


RATIO 

ETA 

^ J. n 09ooo»oo 


. CELERITY 






.14241650 


1.8118475 


.46759916 


.19169182 


-.64377758 


' 






-2;9928239E-03 


5.1401008 






1.4849918 

-4.0869344E-02 


1.3675264 


.54521879 


-5.4364837E-02 


9.6941615E-02 


5. 1401008 






2.2508878 
-7.4961981 E-03 


.42337731 


.39394133 


-.28213441 


.16077232 


5. 1401008 






2. 1403730 


-.68010281 


9.9991 7 18E-02 


-.41031499 


1. 25469 15E-02 
.16077191 


5. 1401008 






1.2179713 


-1.5334495 


-.22132178 


-.38208777 


6; OO48844E-02 
9.694O59OE-02 


5. 1401008 






16389811 


-1.7990193 


-.46944718 


-.20088372 


-.55857185 

-2.9940456E-03 


5. 1401008 






-1.4873984 


-1.3632438 


-.53640372 


5.9416145E-02 


-4; 0400702E-02 

-9.8792210E-O2 


5.1401008 






-2.2309538 

-7.4367725E-03 


-.38955015 


-.38681404 


. 29754649 


-. 15593016 


5. 1401008 






-2. 1256296 


.74113958 


10430827 


.41993030 


1.2424588E-02 
-. 15592911 


5. 1401008 






-1.2287752 


1 . 5808456 


.21154973 


.37494512 


5.9058263E-02 

-9.8789370E-02 


5.1401008 






.14245032 


1.8118447 


.46760341 


.19168646 


-.64362417 

-2.9303792E-03 


5.1401008 
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POTENTIAL SOLUTION - ST OKI AN WAVE 
RUN NUMBER 12 

WL BL H DEPTH Z RHO VISC ALPHA 



SI S2 S3 S4 S5 S6 S7 



5.0 2.2917 


.333 2.0 


8333 1.93 9. 


66E-06 .7854 


.4726 .1319 


=.1529 “.0854 -.1529 0.0 

FORMAT 


0.0 


Vertical 


Vertical 


Lateral 


Lateral 


Force 

AXIAL FORCE 
ETA 


Moment 


Force 


Moment 


.21435714 

.82864214 
=2. 9928241 E-03 


.43734990 


1.6572902 


.16848355 


1.4208121 

.61372989 

9.6941615E-02 


.48830756 


1.2274642 


-6.4276127E-02 


2.0853962 

.17032646 

.16077232 


.34865410 


.34065414 


-.27011388 


1.9441436 
-. 33583372 
.16077192 


7.2655570E-02 


-.67166993 


-.37665020 


1,0537206 

“.71822699 

9.6940590E-02 


=.22896111 


-1 .4364592 


-.34408256 


“.23396772 

“.83137099 

=2,9940458E“03 


=.43863275 


-1.6627479 


-.17915761 


=1.4224905 

=.62559327 

=9.8792210E»02 


=.48012026 


-1.2511911 


5.9538064E-02 


=2.0668344 

“.17492537 

=.15593016 


=.34231067 


-.34985202 


.27786250 


=1.9309969 

.34486083 

=.15592911 


-7.6919930E-02 


.68972420 


.38617909 


=1.0641510 

.72841051 

“9.8789370E=02 


.21998577 


1.4568263 


.34222349 


.21438818 

.82864066 

“2.9903794E=03 


.43735326 


1.6572873 


.16847846 
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POTENTIAL. SOLUTION - SINE WAVE 
RUN NUMBER 12 



WL BL 


H 


DEPTH 


Z 


RHO 


vise 


ALPHA 


SI S2 


S3 


S4 


S5 


S6 


S7 




5.0 2.2917 


• 333 


2.0 -.8333 


1.93 


9. 


66E— 06 


.7854 


.4726 .1319 


-.1529 


-.0854 


-.1529 


0.0 


0.0 





FORMAT 



Vertical Lateral 

Force Force 

ETA 

CELERITY (Last Line Only) 



Vertical 

Moment 



Lateral 

Moment 



.22790985 

.00000000 


1.5170801 


.44532798 


.15885454 


1.4454616 

9.7866442E-02 


1.1326164 


.49232599 


“5.6575 1 22E-02 


2.1108938 

.15835106 


.31552987 


.35127136 


-.25039490 


1.9700325 

.15835068 


-.62207899 


7. 604241 IE-02 


-.34857192 


1.0766825 

9.7865453E-02 


-1.3220736 


-.22823230 


-.31360571 


“.22792559 

“1.2217770E-O6 


-1.5170789 


-.44532962 


-.15885223 


“1 .4454734 
“9.7867431 E-02 


“1.1326089 


-.49232541 


5.6577673E-02 


“2.1108971 

“.15835144 


-.31551891 


-.35126876 


.25039674 


“1.9700261 

“.15835030 


.62208922 


-7.6038782E-02 


.34857232 


“1.0766688 

“9.7864463E-02 


1.3220793 


.22823556 


.31360453 


.22794135 

2.4447195E“06 

5.0599504 


1.5170776 


.44533126 


.15884992 
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EXPERIMENTAL RESULTS 



Run Number 13 



Model Angle of Depth Below Wave Wave Depth of 

Attack Surface Height Length Water 

2-7/8" cyl. 45° 10" 4" 5 . O' 2.0' 



Phase 

Angle 


ETA 

ft. 


Vertical 
Force 
lbs . 


Lateral 
Moment 
ft . lbs 


0 


-.060 


-0.175 


0.062 


36 


.050 


0.220 


- 0.070 


72 


.158 


0.500 


- 0.158 


108 


.171 


0.550 


- 0.222 


144 


.087 


0.340 


- 0.128 


180 


-.020 


-0.025 


-0.015 


216 


-.045 


-0.430 


0.095 


252 


-.148 


- 0.680 


0.251 


288 


-.151 


-0.730 


0.248 


324 


-.125 


-0.525 


0.190 


360 


-.060 


-0.175 


0.060 
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MODIFIED MORRISON SOLUTION - STOKIAN WAVE 



RUN NUMBER NO. 13 



WL 


BL 


H 


DEPTH 


z 


RHO 


Vise ALPHA 


SI 


S2 


S3 


S4 


S5 


S6 


S7 


5.0 


2.5 


.333 


2.0 


-.8333 


1.93 


9.66E-06 .7854 


.23958 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 



FORMAT 



Vertical 



Late red. 



Vertical 



Lateral 



Force 

RATIO 

ETA-CELERITY 
-5.4521970E-02 
.95561037 
-2.9928239E-03 
.47238862 
-9.0653981 E-02 
9.69416158=02 
.81948759 
-2.1 1241 61 E-02 
.16077232 
.85410791 
2.Q268351E-02 
.16077191 
.55803130 
7. 6741 283 E-02 
9.6940590E-02 
4.9675167E-02 
1 .0488501 
-2.9940456E-03 
-.47409111 
-8.9903309E-02 
-9.8792210E-02 
-.81576235 
-2.0989261 E-02 
-.15593016 
-.85000207 
2.01451 19E-02 
=.15592911 
-.55932092 
7.6204914E-02 
“9.8789370E-02 
-5.4509106E-02 
.95583692 
-2.9903792E-03 



Force 



.69850443 

5.1401008 

.57711209 

5.1401008 

.24304751 

5.1401008 

-.17999241 

5.1401008 

-.53378292 

5.1401008 

-.68136093 

5.1401008 

-.56271369 

5.1401008 

-.22052879 

5.1401008 

.20923694 

5.1401008 

.55906414 

5.1401008 

.69850469 

5.1401008 



Moment 


Moment 


.44399876 


4.4500330E-O3 


.38302166 


-.19452101 


.16200587 


-.32232331 


-.11682648 


-.33472755 


-.34266005 


-.21857038 


-.44437083 


-1.3115675E-02 


-.37615274 


.19911275 


-.15740188 


.33371347 


.11280934 


.33736293 


.33558003 


.20862138 


.44399947 


4.4450320E-03 
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POTENTIAL SOLUTION - STOK1AN WAVE 



RUN NUMBER NO. 13 



WL BL 


H DEPTH 


Z RHO 


Vise ALPHA 


SI S2 


S3 s4 


S5 s6 


S7 


5.0 2,5 


.333 2.0 


-.8333 1.93 


9.66E-06 .7854 


.23958 0.0 ( 


3.0 0.0 0.0 0.0 

FORMAT 


0.0 


Vertical 


Vertical 


Lateral 


Lateral. 


Force 

AXIAL FORCE 
ETA 


Moment 


Force 


Moment 


=2. 2478100E-03 
.32068804 
“2.9928241 E-03 


.41254883 


.64137845 


-4.0291 280E-03 


.47851791 

.25817721 

9.6941615E-02 


.33699758 


.51635633 


-. 19201093 


.77719991 

9.8315610E-02 

.16077232 


.12948521 


.19663193 


-.30540568 


.77719801 
“9.8317831 E-02 
.16077192 


“.12948814 


-.19663636 


-.30540498 


.47851305 

“.25817861 

9.6940590E-02 


-.33699935 


-.51635911 


-.19200906 


°2.2538100E“03 

“.32068807 

-2.9940458E-03 


-.41254875 


-.64137851 


-4.0267310E-03 


“.47991197 

“.26070504 

“9.8792210E-02 


-.33051795 


-.52141 198 


.18952277 


“.77356464 

9.98?6470E“02 

15593016 


-.12547887 


-.19975367 


.31192573 


“.77355916 

9.9883238E“02 

“.15592911 


.12548739 


.19976721 


.31192343 


".47989755 

.26070917 

“9.8789370E-02 


.33052333 


.52142027 


.18951685 


“2. 2358300E-03 
.32068802 
“2.9903794E-03 


.41254890 


.64137838 


-4.0338S20E-03 
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POTENTIAL SOLUTION - SINE WAVE 



RUN NUMBER NO. 13 



WL BL 


H DEPTH 


Z RHO 


Vise ALPHA 


SI S2 


S3 S4 


S5 s6 


S7 




5.0 2.5 


.333 2.0 


-.8333 1.93 


9.66E 


-06 .7854 


.23958 0.0 


0.0 0.0 


0.0 0.0 


0.0 






FORMAT 








Vertical 


Lateral 


Vertical 




Lateral 


Force 


Force 


Moment 




Moment 


ETA 










CELERITY (Last 


Line Only) 








.00000000 
. 00000000 


.58615118 


.41946006 




7. 0000000 E-09 


.48724054 

9.7866442E-02 


.47420577 


.33934995 




-.17433956 


. 78837093 


.18112906 


.12961911 




-.28208701 


.15835106 

.78836904 


-.18113314 


-.12962204 




-.28208636 


.15835068 

.48723562 


“.47420830 


“.33935175 




-.17433780 


9.7865453E-02 

-6.0600000E-06 


-.58615120 


-.41946005 




2 . 1 830000E-06 


-1.2217770E-06 

“.48724547 


-.47420324 


-.33934814 




.17434131 


=9. 7867431 E“02 
“.78837280 


-.18112495 


-.12961618 




.28208769 


“.15835144 

“.78836716 


.18113723 


.12962497 




.28208569 


=.15835030 
“.48723070 
=9. 7864463 E“02 


.47421085 


.33935357 




.17433604 


1 .2150000E“05 
2.4447195E“06 
5.0599504 


.58615118 


.41946002 




-4.3530000E-06 
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EXPERIMENTAL RESULTS 



Run Number 14 



Model 


Angle of 
Attack 


Depth Below 
Surface 


Wave 

Height 


Wave 

Length 


Depth of 
Water 


Sub 


90° 


10" 


4" 


5.0' 


2.0' 



Phase 

Angle 


ETA 

ft. 


Vertical 

Force 

lbs. 


Lateral 

Force 

lbs. 


Lateral 
Moment 
ft . lbs . 


0 


-.040 


-0 . 000 


2.970 


0 . 668 . 


56 


.018 


1.125 


2.460 


0.450 


72 


.167 


1.875 


0.524 


0.150 


108 


.180 


1.680 


-1.555 


-0.550 


144 


.090 


0.600 


-2.740 


-0.525 


180 


.006 


-0.525 


-2.810 


-0.574 


216 


-.082 


-1.400 


-1.950 


-0.518 


252 


-.128 


-1.875 


-0.581 


-0 . 158 


288 


-.150 


-1.825 


0,647 


0.175 


524 


-.150 


-1 . 280 


2.050 


O.58O 


560 


-.082 


-0.100 


2.900 


0.605 
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MODIFIED MORRISON SOLUTION - STOKIAN WAVE 



RUN NUMBER NO. 14 



WL 


BL 


SI 


S2 


5.0 


2.2917 


.4725 


=.1319 



Vertical 

Force 

RATIO 

ETA-CELERITY 



=.13035610 
.88040191 
-2.9928239E-03 
1 .3976820 
=5.2303486E=O2 
9.6941615E=02 
2.3905812 
=4.0181?87E=03 
.16077232 
2.4097864 
3.9863515E-03 
.16077191 
1.5438726 
4.7351346E=02 
9.6940590EH32 
9.9154320E-02 
1.1574458 
=2.9940456E=03 
=1.4084771 
=5.1089904E»02 
=9.8792210E=02 
=2.3655876 
=3.9639019E=03 
=.15593016 
=2.3843253 
3.9333142E=03 
=.15592911 
=1.5523515 
4. 63 5649 3 E=0 2 
=9.8789370E=02 
=.13031927 
.88065072 
=2.9903792E=03 



H DEPTH 
S3 S4 

.333 2.0 

=.1529 .0854 



FORMAT 

Lateral 

Force 



2.7810507 

5.1401008 

2.2825804 

5.1401008 

.97067518 

5.1401008 

-.71302742 

5.1401008 

= 2.1558606 

5.1401008 

=2.7660623 

5.1401008 

=2.2724493 

5.1401008 

=.89879094 

5.1401008 

.84581929 

5.1401008 

2.2645414 

5.1401008 

2.7810505 

5.1401008 



Z RHO 

S5 S6 

-.8333 1.93 

-.1529 0.0 

Vertical 

Moment 

1.1714574E-02 
-.16533858 
-.27920571 
-.28079203 
-.17781536 
-8.0747 160E-03 
.16656457 
.27628350 
.27782829 
.17883774 
1.1710319E-02 



Vise ALPHA 

S7 

J.66E-06 1.5708 

).0 

Lateral 

Moment 

=.32546625 

-.26537561 

-.10979980 

8.7388972E-02 

.25443751 

.32419009 

.26599728 

.10439476 

-9.9927040E-02 

-.26535734 

-.32546621 
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POTENTIAL SOLUTION - STOKIAN WAVE 



RUN NUMBER NO. 14 



WL BL 


H DEPTH Z 


RHO 


Vise ALPHA 


SI S2 


S3 S4 S5 


s6 


S7 


5.0 2.2917 


.333 2.0 -.8333 1.93 


9.66E-06 1.5708 


.4725 -=.1319 ■ 


-.1529 .0854 -. 

FORMAT 


1529 0.0 


0.0 


Vertical 


Vertical 


Lateral 


Lateral 


Force 

AXIAL FORCE 
ETA 


Moment 


Force 


Moment 


“1 .4343673E-02 
=4. 7140246E-06 
=2.9928241 E-03 


1 .6637308E-03 


2.5696509 


-.29846216 


1.3531741 

“3.7718456E-06 

9.6941615E-02 


-.15717159 


2.0560618 


-.23880787 


2.2082557 
=1 . 43081 82E-06 
.16077232 


-.25648726 


.77995005 


-9.0588050E-02 


2.2082495 
1 .4308553E=06 

.16077192 


25648513 


-.77997022 


9.0595010E-02 


1.3531582 

3.7718689E-06 

9.6940590E-02 


15716606 


“2.0560748 


.23881233 


“1.4363073E-02 
4.7 14025 2 E~06 
“2.9940458E=03 


1 .6704736E-03 


-2.5696513 


.29846226 


“1.3620530 
3 . 8555752E-06 
-9 . 87922 10E-02 


.15820277 


-2.1017036 


.24410920 


“2.1850503 
1 .4825446E-06 
15593016 


.25379193 


-.80814642 


9.3862940E-02 


“2,1850342 
“1 .4826496E-06 
=. 1559291 1 


.25378870 


.80820372 


-9.3874370E-02 


“1 .3620105 

“3.8556383E-06 

“9.8789370E-02 


.15819431 


2. 1017380 


-.24411604 


“1 .4309783E-02 
4. 7 1 40232E“06 
-2.9903794E-03 


1 .6598486E-03 


2.5696501 


-.29846208 
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POTENTIAL SOLUTION - SINE WAVE 
RUN NUMBER NO. 14 



WL BL 


H DEPTH 


Z RHO VISC 


ALPHA 


SI S2 


S 3 s4 


S5 S6 S7 




5.0 2. 29117 . 


3333 2.0 


-.8333 1.93 9.66E 


-06 1.5708 


.4725 -=.1319 *=.1529 .0854 


-.1529 0.0 0.0 






FORMAT 




Vertical. 


Lateral 


Vertical 


Lateral 


Force 


Force 


Moment 


Moment 


ETA 








CELERITY (Last Line Only) 






1 .258J347E-06 


2.350501 1 


-2.45791 17E-06 


-.27300818 


.00000000 
1 .3815936 


1 .9015927 


-.16047248 


-.22086658 


9.795461QE-02 

2.2354619 


.72633710 


-.25964720 


-8.4361 112E-02 


.15849372 

2.2354559 


=.72635584 


-.25964510 


8.4367592E-02 


.15849334 

1.3815778 


=1.9016042 


-.16046694 


.22087059 


9.7953620E-02 
-1. 848391 9E-05 
~1.2228777E~06 


=2.350501 1 


4.4042005E-06 


.27300818 


=1.3816076 

-9.7955600E-02 


=1.9015823 


.16047406 


.22086543 


=2.2354672 

=.15849410 


=.72632068 


.25964781 


8.4359212E-02 


=2.2354505 


.72637227 


.25964445 


-8.4369492E-02 


=.15849296 

=1.3815637 

=9.7952629E=02 


1 .9016144 


.16046531 


-.22087179 


3.5747789E=05 

2.4469220E=06 


2.3505011 


-6.408877 IE-06 


-.27300818 


5.0599504 
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MODIFIED MORRISON SOLUTION - ST OKI AN WAVE 



RUN NUMBER 14 



WL BL 


H DEPTH, 


Z RHO VISC 


ALPHA 


SI S2 


S3 S4 


S5 S6 S7 




5.0 2,2917 


.333 2.0 


8333 1.93 9.66E 


-06 1.5708 


,4726 ,1319 


-,1529 -.0854 


-.1529 0.0 0 


.0 




FORMAT 




Vertical 


Lateral 


Vertical 


Lateral 


Force 


Force 


Moment 


Moment 


RATIO 








ETA 

" «*• 000000000 


CELERITY 






“,13040307 


2.7824348 


-1.1721 280E-02 


•32557516 


,88036485 








“2.9928239E-03 


5. 1401008 






1 ,3983871 


2.2837083 


.16539219 


. 26546633 


“5- 2 29 3609 E -02 








9.6941615E-02 


5. 1401008 






2,3917694 


.97113828 


.27929969 


. 10983879 


-4,0173267E-03 








,16077232 


5. 1401008 






2; 4 1098 18 


“.71340309 


.28088945 


-8.7415700E-02 


3.9854807E-03 








,16077191 


5.1401008 






1.5446282 


-2. 1569441 


.17787858 


-.25452046 


4.7343124E-02 








9:6940590E-02 


5.1401008 






9'. 91 85720E-02 


“2.7674421 


8.0800517E-03 


-.32429866 


1 .1574476 








-2.9940456E-03 


5.1401008 






-I. 4091875 


-2.2735824 


-.16661884 


-.26608952 


=5. 1080306E-02 








“9» 87922 10E-02 


5.1401008 






-2.3667636 


-.89923661 


-.27637657 


-. 10443455 


“3.9630652E-03 








-.15593016 


5.1401008 






-2.3855079 


.84624299 


-.27792466 


9. 9955973E-02 


3.9324517E“03 








-.15592911 


5.1401008 






-1-.55311 13 


2.2656690 


-.17890115 


.26544327 


4.6348326E-02 








-9.8789370E-02 


5.1401008 






“'.13036613 


2.7824346 


-1.1716937E-02 


.32557516 


.88061415 








-2.9903792E-03 


5. 1401008 







- 158 - 



POTENTIAL SOLUTION - ST OKI AN WAVE 
RUN NUMBER 14 



WL BL 


H DEPTH Z 


RHO 


vise 


ALPHA 


SI S2 


S3 S4 S5 


S 6 


S7 




5.0 2.2917 


.333 2.0 -.8333 


1.93 


9.66E 


-06 1.5708 


.4726 .1319 


-.1529 -.0854 -.1529 0.0 

FORMAT 


0.0 


Vertical 


Vertical 


Lateral 




Lateral 


Force 

AXIAL FORCE 
ETA 


Moment 


Force 




Moment 


-1.4353223E-02 
-4.7163498E-06 
=2.9928241 E-03 


-1.6691295E-03 


2.5709183 




.29856379 


1.3538396 

“3.7737089E-06 

9.6941615E-02 


.15722112 


2.0570776 




.23889235 


2.2093442 

-1.4315287E-06 

.16077232 


.25657303 


.78033725 




9. 0623791 E-02 


2.2093395 
1 .4315564E-06 
.16077192 


.25657389 


-.78035240 




-9.0620931 E-02 


1.3538276 

3.7737266E-06 

9.6940590E=02 


.15722341 


=2.0570872 




-.23889047 


=1.4367740E=02 

4.7163503E=06 

-2.9940458E-03 


-1.6663234E-03 


=2.5709186 




-.29856376 


=1.3627229 

3.8574796E=06 

=9.8792210E=02 


=. 15825284 


-2.1027417 




-.24419526 


=2.1861273 
1 ,4832805E=06 

=.15593016 


=.25387688 


-.80854764 




-9.3899964E-02 


=2.1861126 

-1.4833762E-06 

=.15592911 


-.25387652 


.80859976 




9. 3901 248E-02 


=1.3626843 

=3.8575372E=06 

=9.8?89370E=02 


=, 15825189 


2. 1027730 




.24419606 


-1.4319258E-02 

=4.7163486E=06 

=2.9903794E=03 


-1.6651312E-03 


2.5709177 




.29856371 
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POTENTIAL SOLUTION - SINE WAVE 



RUN NUMBER 14 



WL BL 


H DEPTH 


z 


RHO VISC 


ALPHA 


SI S2 


S3 S4 


S5 


S6 S7 




5.0 2.2917 


.333 2.0 ■ 


-.8333 1.93 9.66E 


-06 1.5708 


.4726 .1319 


-.1529 “.0854 


-.1529 0.0 


0.0 




FORMAT 






Vertical 


Lateral 




Vertical 


Lateral 


Force 

ETA 


Force 




Moment 


Moment 


CELERITY (Last Line Only) 








-1.2580297E-06 


2.3495438 




-2.4572832E-06 


.27285526 


.00000000 

1.3810288 

9.7866442E-02 


1 .9008195 




.16037862 


.22074575 


2.2345507 


.72604359 




.25950027 


8.4318433E-02 


.15835106 

2.2345462 


-.72605764 




.25950115 


-8.4315783E-02 


.15835068 

1.3810170 


-I .9008282 




.16038093 


-.22074405 


9.7865453E-02 
-1.6018482E-05 
~1 . 221 7770E-06 


-2.3495438 




3. 961 3342E-07 


-.27285526 


=1.3810429 


-1.9008094 




-.16038031 


-.22074454 


“9. 786743 IE-02 
“2.2345560 


-.72602716 




-.25950091 


-8.4316529E-02 


-.15835144 

“2.2345408 

”.15835030 


.72607407 




-.25950054 


8.4317687E-02 


“1. 3810030 


1.9008383 




-.16037932 


.22074525 


“9.7864463E-02 










3.3275712E-05 

2.4447195E-06 

5.0599504 


2.3495438 




1 .6081609E-06 


.27285526 
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EXPERIMENTAL RESULTS 



Run Number 15 



Model Angle of Depth Below Wave Wave Depth of 

Attack Surface Height Length Water 



2 - 7 / 8 " cyl. 90 ° 


10" 


4 " 


5. O' 2.0 ’ 


Phase 

Angle 


ETA 

ft. 


Vertical 
Force 
lbs . 


Lateral 
Force 
lbs . 


0 


-.04 


-0 . 280 


1.344 


56 


.018 


0.175 


1.725 


72 


.167 


0.480 


1.270 


108 


.180 


0.530 


0.000 


144 


.090 


0.345 


-1.190 


180 


-.006 


0.005 


-1.672 


216 


-.082 


- 0.345 


-1.780 


252 


-.128 


-0.610 


- 1.225 


288 


-.150 


- 0.805 


- 0.291 


324 


-.130 


-0.620 


O.618 


360 


-.082 


-0.300 


1.620 
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MODIFIED MORRISON SOLUTION - STOKIAN WAVE 
RUN NUMBER NO. 15 



WL BL 

SI S2 

5.0 2.5 

.23958 0.0 



Vertical 

Force 

RATIO 

ETA-CELERITY 

-8.4924284E-02 
.92041975 
-2.9928239E-03 
.58738901 
-8. 5342747 E~02 
9.6941615E-02 
1 .0333142 
=6. 830347 IE-03 
.16077232 
1.0474277 
6.7386139E-03 
.16077191 
.68764178 

7.2901237E-02 
9.6940590E-02 
7.1399702E-02 
1.0947662 
-2.9940456E-03 
-.59234245 
-8.3328616E-02 
-9.8792210E-02 
-1.0225424 
-6.7448044E-03 
-.15593016 
-1 .0363297 
6.6559595E-03 
-.15592911 
-.69104321 
7.1429343E-02 
=9.8789370E-02 
=8.49083 18E-02 
.92059300 
-2.9903792E-03 



H DEPTH 

S3 S4 

.333 2.0 

0.0 0.0 



FORMAT 

Lateral 

Force 



1.2068576 

5.1401008 

1.0044417 

5.1401008 

.45211919 

5.1401008 

-.27722200 

5.1401008 

-.91818883 

5.1401008 

-1.1960239 

5.1401008 

-.98560006 

5.1401008 

<=.39617175 

5.1401008 

.35955338 

5.1401008 

.97972001 

5.1401008 

1.2068577 

5.1401008 



Z RHO 

S5 s6 

-.8333 1.93 

0.0 0.0 

Vertical 
Moment 

-2.6120000E-06 
-2.3280000E-06 
-9.1000000E-07 
7.1000000E-07 
1 .9400000E-06 
2.5630000E-06 
2.2090000E-06 
9. 1000000E-07 
=6. 7 000000 E -07 
-1.8400000E-06 
-2.5610000E-06 



Vise ALPHA 

S7 

9.66E-06 1.5708 

0.0 

Lateral 
Moment 

-6. 0000000 E -08 
1 .5200000E-06 
2.5700000E-06 
2.7860000E-06 
1 . 8800000E-06 
1.1000000E-07 
-1 . 6100000E-06 
-2. 69 30000 E -06 
-2.8250000E-06 
-1.7300000E-06 
-1 .0000000E-O8 
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POTENTIAL SOLUTION - STOKIAN WAVE 



RUN NUMBER NO. 15 



WL 


BL 


H 


DEPTH 


z 


RHO 


vise 


ALPHA 


SI 

1 


S2 


S3 


S4 


S5 


S6 


S7 




5.0 


2.5 


.333 


2.0 


“.8333 


1 .93 


9.66E-06 


1 .5708 


.23958 0.0 


0.0 


0.0 


0.0 


0.0 


0.0 







FORMAT 






Vertical 


Vertical 


Lateral 


Lateral 


Force 

AXIAL FORCE 
ETA 


Moment 


Force 


Moment 


“6.2?69549E“03 
»2.0627328E“06 
=2.9928241 E-03 


“2.4270400E-06 


1.1244115 


4. 0000000E-08 


.5921 1291 
-1.6504608E-06 
9.6941615E=02 


-1.9860000E“06 


.89967887 


1.6000000E-06 


.96627438 

-6.2608937E-07 

.16077232 


“7.5000000E-07 


.34128612 


2.4880000E-06 


.96627204 

6.261Q355E-Q7 

.16077192 


7.5000000E=07 


=.34129384 


2.4760000E-06 


.59210682 
1 .6504697E“06 
9.6940590E“02 


1 .9520000E-06 


-.89968374 


1.5800000E-06 


=6. 2843737 E-03 
2.0627330E-O6 
=2. 9940458 E=03 


2.3788500E-06 


-1.1244116 


4.0000000E-08 


“.59599813 
1 .6870986EHD6 
=9.8792210E“02 


1 ,9000000E“06 


-.91965041 


-1.5300000E-06 


=.95612028 
6.4872346E-07 
“. 15593016 


7.0000000E=07 


-.35362413 


-2.5540000E-06 


“.95611358 

“6.4876735E=07 

“.15592911 


“7. 1000000E-07 


.35364806 


-2.5550000E-06 


=.59598038 
“1 .6871250E-06 
“9.8789370E-02 


“1 .9070000E-06 


.91966481 


-1 .5500000E-06 


=6.2621 124E=03 
=2.0627322E=O6 
“2.9903794E“03 


“2.3788400E“06 


1.1244112 


5.0000000E-08 
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POTENTIAL SOLUTION - SINE WAVE 



RUN NUMBER NO. 15 



WL 


BL 


H 


DEPTH 


z 


RHO 


vise 


SI 


S2 


S3 


S4 


S5 


S6 


S7 


5.0 


2.5 


.333 


2.0 


.8333 


1.93 


9.66E 


.23958 0.0 


0.0 


0.0 


0.0 


0.0 


0.0 



FORMAT 



Vertical 


Lateral 


Vertical 


Force 

ETA 


Force 


Moment 


CELERITY (Last Line Only) 




2.0000000E-13 

.00000000 


8.3439272 


-2.0038592E-05 


4.9044471 

9.7866442EH32 


6.7503716 


-1 .6060000E-05 


7.9355537 

.15835106 


2.5783919 


-6.0000000E-06 


7.9355348 

.15835068 


-2.5784501 


6. 1000000E-06 


4.9043975 

9.78654535=02 


=6.7504074 


1 .5920000E-05 


=6.1252766E“05 

~1.2217770E=06 


-.3439272 


1 .9640178E“05 


=4.9044966 
=9. 7867431 E=02 


=6.7503354 


1 .5980000E-05 


=7.9355726 

=.15835144 


“2.5783335 


6.0000000E-06 


=7.9355159 

“.15835030 


2.5785082 


-6. 1000000E-06 


=4.9043480 
=9. 7864463 E=02 


6.7504436 


-1 .6010000E-05 


1 .2253890E=04 
2.4447195E=Q6 
5.0599504 


8.3439272 


-1 .96401 77 E-05 



ALPHA 
06 1.5708 



Lateral 

Moment 



.00000000 
1 . 1700000E-05 
1 .8690000E-05 
1 .8690000E=05 
1 .1500000E-05 
.00000000 
1 . 1700000E-05 
1 .8690000E-Q5 
1 .8690000E-05 
K1500000E-05 
.00000000 
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APPENDIX B. PROTOTYPE RUNS - THEORY ONLY 



This appendix contains the tables of results obtained 
by each of the three theories for the two prototype bodies 
investigated. These tables should be interpreted in the 
same manner as the tables in Appendix A, 

The tables included are; 



Run 

No. 


Body 


Angle of 
Attack 


Wave 

Length 


Wave 

Height 


PI 


Sub 


0° 


400 


25 


P2 


Sub 


45° 


400 


25 


P3 


Sub 


90° 


400 


25 


P4 


5 ft. Cyl. 


0° 


400 


25 


P5 


5 ft. Cyl. 


45° 


400 


25 


P6 


3 ft. Gyl. 


90° 


400 


25 
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MODIFIED MORRISON SOLUTION - ST OKI AN WAVE 
RUN NUMBER PI 



WL BL 

400. 250. 

SI 

31.1558 

S5 

8.8705E-08 



H DEPTH Z RHO VISC ALPHA 

25. 900. -80. 1.93 9.66E-06 .00001 



S2 S3 

4.0329E-02 -1.6597E-03 



S6 

=5. 18Q5E-10 



S7 

-7.3899E-12 



S4 

3.0559E-06 



FORMAT 



Vertical 


Lateral 


Force 


Force 


RATIO 


ETA ......... 




74538.300 


6,5026223 


1068301 1 


=1,7792336E~06 


46.121821 


435068.94 


4.7691744 


=1 .4848854E-02 


7.3473300 


46. 121821 


629486.25 


1,2140561 


3. 8461091 E-03 


11,888219 


46,121821 


584478,32 


=2.8047922 


6, 0928432E-03 


1 1.888191 


46. 121821 


315352,20 


-=5.7522996 


2.3201865E-02 


7,3472560 


46.121821 


=74542.960 


=6,5026138 


=. 10682333 


=9.3504230E=05 


46.121821 


=435072.43 


-4.7691390 


=1 . 484866QE-02 


=7,3474055 


46.121821 


=629487.17 


=1.2140064 


“3. 8459941 E =03 


=11.888245 


46. 121821 


=584476.37 


2.8048376 


6.0929667E-03 


=11.888160 


46.121821 


-315348.08 


5.7523250 


2.3202257E-02 


=7.3471825 


46. 121821 


74547,674 


6.5026104 


=. 10681649 


1 .8175827E-04 


46. 121821 



Vertical 

Moment 


Lateral 

Moment 


29896485. 


36.580903 


26619339. 


-150.52323 


13171736. 


-280. 13237 


-5288808. 4 


-302.74000 


-21737195. 


-209. 71068 


-29896515. 


-36.578644 


-26619236. 


150.52526 


-13171538. 


280.13332 


5289027.2 


302.73953 


21737350. 


209,70899 


29896545. 


36.576416 
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POTENTIAL SOLUTION - STOKIAN WAVE 



RUN NUMBER PI 



WL BL 


H DEPTH 


Z RHO 


VISC ALPHA 


400. 250. 


25. 900. 


-80. 1.93 


9. 66E-06 0.0 


SI 


S2 


S3 


S4 


31.1558 


4.0329E-02 


=1.6597E=03 


3.0559E-06 


S5 


S6 


S7 




8.87Q5E-08 


=5. 1805E-10 


”7.3899E=12 





FORMAT 



Vertical 


Vertical 


Lateral 


Lateral 


Force 

AXIAL FORCE 
ETA 


Moment 


Force 


Moment 


75340.444 

292910.92 

-1.7792336E-06 


27606684. 


.00000000 


.00000000 


405289.76 

214827.62 

7.3473300 


24271339. 


. 00000000 


. 00000000 


580431.46 

54687.091 

11,888219 


1 1665127- 


. 00000000 


.00000000 


533867 . 10 
=126342. 14 
11.888191 


=5396787.6 


. 00000000 


,00000000 


283382,68 

=259112.74 

7.3472560 


=20397303. 


. 00000000 


. 00000000 


=75344.747 

=292910,64 

»9,3504230E=05 


=27606708, 


. 00000000 


.00000000 


=405292.94 

=214826,14 

=■7.3474055 


=24271240. 


. 00000000 


. 00000000 


=580432.26 

=54684.958 

=11,888245 


=11664942. 


. 00000000 


.00000000 


=533865.24 
126344. 10 
=11.888160 


5396988.8 


. 00000000 


. 00000000 


=283378.88 

259113.78 

=7.3471825 


20397443. 


. 00000000 


. 00000000 


75349.048 

292910.36 

1.8175827E=04 


27606732. 


. 00000000 


. 00000000 



- 167 - 



MODIFIED MORRISON SOLUTION - ST OKI AN WAVE 



RUN NUMBER P2 



WL BL 

400. 250, 

SI 

31.1558 

S5 

8.8705E-O8 



Vertical 
Force 
RATIO 
ETA .... 
52398.946 
=.20312320 
=1.7792336E-06 
490911.25 
=1.6581949E=02 
7.3473300 
741381.64 
=4. 1239994E-03 
11.888219 
709285.35 
5.3722337E=03 
11.888191 
406039.31 
2.2158463E=Q2 
7.3472560 
=52404.544 
=.20310139 
=9.350423OE=O5 
=490915.57 
=1.6581720E=02 
=7.3474055 
=741382.99 
=4, 1238890E=03 
=11,888245 
=709283.28 
5.3723510E=03 
= 11.888160 
=406034.55 
2.2158834E=02 
=7.3471825 
52410. 149 
=.20307938 
1.8175827E=04 



H DEPTH 

25. 900. 

S2 

4.0329E=02 

S6 



. . . » CELERITY 
540860.11 

46.121821 

413640.46 

46.121821 

129083.47 

46. 121821 
=204576.22 

46. 121821 
“459754.37 

46. 121821 
=538490,31 

46.121821 
=410808. 17 

46. 121821 
=125447.79 

46. 121821 
208106.34 

46. 121821 
462362.04 

46. 121821 
540859.52 

46.121821 



Z RHO 

- 80 . 1.93 

S3 

-1.6597E-03 

S7 

-7.3899E-12 



23791644. 
22396650. 
12302205. 
-2392214. 1 
=16094864. 
-23791683. 
-22396577. 
= 12302050 . 
2392391.6 
16094997. 
23791723. 



=5. 1805E-10 

FORMAT 

Vertical 
Moment 



Lateral 

Force 



VISC ALPHA 

9.66E-06 .785 

S4 

3.0559E-06 



Lateral 

Moment 

3519782.3 
-7004072.0 
-14901050. 
=17150768. 
-12801209. 
-3509499.7 

7092794.7 

14959752. 

17116605. 

12731170. 

3519536.3 
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POTENTIAL SOLUTION - STOKIAN WAVE 



RUN NUMBER P2 



WL BL 


H DEPTH 


Z RHO 


.vise 


ALPHA 


400, 250. 


25. 900. 


“80. 1.93 


9.66E-06 . 


785 


SI 


S2 


S5 


S4 




31.1558 


4.0329E-02 


-1.6597E-03 


3.0559E^06 




- S5 


S6 


S7 






8.87Q5E-08 


=5. 1805E-10 


-7.3899E-12 







FORMAT 



Vertical 


Vertical 


Lateral 


Lateral 


Force 

AXIAL FORCE 
ETA 


Moment 


Force 


Moment 


57572.878 

247897.01 

-1.7792336E-06 


21999681. 


495399.38 


3102645. 1 


458544.91 

188583.47 

7.3473300 


20378215. 


376866.73 


=6629946.5 


684367,62 

57237,141 

11.888219 


10972928. 


114383.16 


-13830112. 


648783.94 

=95971,949 

11.888191 


-2623661,3 


-191791.09 


-15747622. 


365385.75 

=212522.83 

7.3472560 


-15218097. 


=424707.31 


-11650047. 


=57578.005 

=247896.87 

=9.350423OE=O5 


=21999712. 


-495399.08 


-3102531.0 


=458548.87 

=188582.28 

=7.3474055 


=20378146. 


=376864.37 


6630053.3 


=684368.79 

=57235.360 

=11.888245 


=10972785. 


=114379.60 


13830169. 


=648781.97 

95973.633 

=11.888160 


2623826.5 


191794.46 


15747607. 


=365381.34 

212523,79 

=7.3471825 


15218218, 


424709.23 


11649969. 


57583, 162 
247896,72 
1.8175827E-04 


21999744. 


495398.79 


3102417.2 
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MODIFIED MORRISON SOLUTION 
RUN NUMBER ?J> 



WL BL 

400. 250, 

SI 

31.11558 

S5 

8, 8765E-Q8 



Vertical 

Force 

RATIO 

eta .oooo.oo 

-15837.002 

.99997777 

-1.7792336E-06 

520359.89 

-1.9695024E-Q2 

7.3473300 

857165.51 

-1.6057055E-03 

11,888219 

859916.44 

I . 6006424E-03 

II, 888191 
540852,19 

1. 89490 10E -02 
7.3472560 
15830.375 
1.0003963 
~9. 35042 30E-Q5 
-520365.37 
- 1.9 69459 QE-02 
=7.3474055 
-857167.63 
-1.6056246E-03 
-11.888245 
-859914.46 
1. 60072 14E-03 
- 11,888160 
=540846.97 
1.894939 IE-02 
-7.3471825 
-15823.744 
1.0008156 
1.8175827E-04 



H DEPTH 
25. 900. 

S2 

4.0329E-02 



. . o CELERITY 
899610.69 

46.121821 

733516.35 

46.121821 

290451.46 

46.121821 

-265169.60 

46.121821 

-721117.47 

46.121821 

-898412.64 

46.121821 
=727996. 19 

46.121821 

-281360.54 

46.121821 

274260,50 

46. 121821 

726637.61 

46.121821 

899610.77 

46.121821 



Z RHO 

-80. 1,93 

S3 

=1 .6597E-03 
S 7 

=7. 3899E-1 2 



-56209.530 
3559510.8 
5813676.2 
5824150.7 
3633037.0 
56163.450 
-3559548.0 
—5813690. 8 
-5824137.9 
-3633000.5 
- 56118.520 



ALPHA 

1.5708 



0559E-06 



Lateral 

Moment 

6099689.8 

4955091.8 

1928641.7 
=1839775.3 
-4910747.2 
-6095139.5 
=4935552.2 
-1897139.8 

1871270.9 

4930283.7 
6099692.3 



S6 

-5.1805E-10 



FORMAT 



Lateral Vertical 

Force Moment 



ST OKI AN WAVE 



vise 

9. 66E-06 
S4 

3. 
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POTENTIAL SOLUTION - STOKIAN WAVE 
RUN NUMBER P3 



WL BL 


H DEPTH 


z RHO , 


VISC ALPHA 


400, 250, 


25. 900. 


-80. 1.93 


9.66E-06 1.5708 


SI 


S 2 


S3 


S4 


311 o 1)558 


4.0329E=02 


-1.6597E-03 


3.0559E-06 


S5 


S6 


S7 




8.87O5E-08 


=5.1805E=10 -7.3899E=12 

FORMAT 




Vertical 


Vertical 


Lateral 


Lateral 


Force 

AXIAL FORCE 
ETA 


Moment 


Force 


Moment 


<=,3220380 
Hi; 5237554 
=1.7792336E=06 


=128.36620 


830610.71 


5587767.5 


48822046 

=1,2327430 

70473300 


3284310. 1 


671977.65 


4520668. 8 


789958,39 

=,47086261 

11,888219 


5314249.5 


256670.79 


1726820.2 


789956,70 

,47087210 

11,888191 


5314313.9 


-256675.99 


-1726617.0 


488217.03 

1,2327488 

7.3472560 


3284482.9 


-671980.85 


-4520543.9 


-5; 7766547 
1.5237554 
=9.350423OE=O5 


86.243699 


=830610.71 


“5587767.5 


=488226,40 

1,2327364 

=7.3474055 


=3284343.9 


-671974.07 


-4520644.0 


=789960,25 

.47085194 

=11.888245 


=5314261.6 


-256665.00 


-1726781.5 


=789954.82 

=.47088275 

=11.888160 


=5314301.8 


256681.80 


1726656.5 


=488212. 10 
=1.2327554 
=7.3471825 


=3284450.3 


671984.46 


4520568.2 


11.877490 

=1.5237554 


=45.201560 


830610,71 


5587767.5 



1.8175827E-04 
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MODIFIED MORRISON SOLUTION - ST OKI AN WAVE 



RUN NUMBER P4 



WL 


BL 


H 


DEPTH 


Z 


RHO 


vise 


ALPHA 


SI 


S2 


S3 


S4 


S5 


S6 


S7 




400. 


30. 


25. 


200. 


=50. 


1.93 


9. 66E-06 


.00001 


3.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 





FORMAT 



Vertical 


Lateral 


Vertical 


Lateral 


Force 


Force 


Moment 


Moment 


RATIO 


ETA .00000.090 


CELERITY 




' 


^32.3^097 


2.71 30823E=02 


1138.6820 


-7.51 17000E-05 


.99641 D 65 


=.10785793 


46.045438 


- 




226.05225 


2. 1919489E-02 


1547.5368 


=1 . 8881820E-02 


=1.5069499 


7.3140010 

847.39942 


46.045438 

8.3656598E-03 


608.97260 


-3. 04531 36E-02 


=8. «0 1624 IE-02 


11.975477 

991.48486 


46.045438 

=8.3650646E=03 


=101.79570 


-3. 045 3407 E-02 


7.2663654E-Q2 


11.975448 

907.34958 


46.045438 
=2. 19191 31 E-02 


-307.91550 


-1.8882354E-02 


.37543818 

7.3139255 


46.045438 






528.52299 


=2. 7 130 69 6E =02 


=1 144.0439 


-7.5401000E-0C 


1.0036277 


=.10794965 

=228.91758 


46.045438 
=2. 1979087E-02 


-153ft. 5595 


1.8835693E-02 


=1.4807528 


=7.3807330 

=844.70425 


46.045438 

-8.4021924E-03 


-601.61430 


3.0575267E-02 


=8.4821649E=02 


=11.800986 

=988.00271 


46.045438 

8.4026244E-03 


98.407200 


3.0575219E-02 


7.2525922E-02 


=1 1 . 800904 


46.045438 


- 




=906.85234 


2. 19794156=02 


299.51630 


1. 883531 3E-02 


.37380010 


=7.3805145 


46.045438 


- 




=532.33647 


2.71 30821 E=02 


1138.6967 


=7.5597000E=05 


.99643817 


=.10767439 


46.045438 
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POTENTIAL SOLUTION - ST OKI AN WAVE 
RUN NUMBER P4 



WL 


BL 


H 


DEPTH 


z 


RHO 


vise 


ALPHA 


SI 


S 2 


S 5 


S 4 


S 5 


S6 


S 7 




400 , 


30 . 


25 . 


200 , 


- 50 . 


1.93 


9 . 66 E -06 


0.0 


3.0 


0,0 


0,0 


0.0 


0.0 


0.0 


0.0 





FORMAT 



Vertical 


Vertical 


Lateral 


Lateral 


Force 

AXIAL FORCE 
ETA 


Moment 


Force 


Moment 


= 4.7516730 
1222. 1132 
=.10785793 


2839.0509 


.00000000 


.00000000 


1409.6328 

987.35929 

7.3140010 


2307.6479 


.00000000 


.00000000 


2287.0513 

376,81598 

11.975477 


883.98800 


.00000000 


.00000000 


2287.0459 

= 376.82450 

11.975448 


= 884.00800 


.00000000 


.00000000 


1409.6185 

= 987.36455 

7.3139255 


= 2307.6598 


.00000000 


.00000000 


= 4.7692880 
= 1222. 1132 
=. 10794965 


= 2839.0505 


.00000000 


.00000000 


= 1412.5836 

= 990.05408 

= 7.3807330 


- 2286.0161 


.00000000 


.00000000 


= 2279.3685 

= 378.47613 

= 11.800986 


- 870.60640 


.00000000 


.00000000 


= 2279.3522 
378.50180 
= 1 1 , 800904 


870.66540 


.00000000 


.00000000 


= 1412.5410 

990.06990 

= 7.3805145 


2286.0527 


.00000000 


.00000000 


= 4.7163930 
1222. 1131 
=,10767439 


2839.0510 


.00000000 


.00000000 
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MODIFIED MORRISON SOLUTION - ST OKI AN WAVE 



RUN NUMBER P5 



WL 


BL 


H 


DEPTH 


z 


RH0 


vise 


ALPHA 


SI 


S2 


S3 


S4 


S5 


S6 


S7 




400, 


o 

o 

rr\ 


25. 


200. 


=60. 


U93 


9. 66E-06 


. ?8S 


3.0 


o.o 


0.0 


O 

o 

o 


0.0 


0.0 


0.0 





FORMAT 



Vertical 




Lateral 


Vertical 


Lateral 


Force 




Force 


Moment 


Moment 


RATIO 

Si A 090.0 


O O P 0 


CELERITY 






=637.66629 




719.63697 


807 . 69670 


49.124700 


. 9963800^ 

=. 10786793 




46.046438 


- 




227.29197 
=1 . 6048260 




698.46102 


1101.9172 


=124.17170 


7.3140010 




46.046438 


- 




849. 16468 
=8.78l8059E- 


-02 


493.76306 


411.36610 


=426.09400 


11.976477 




46.046438 






998.30831 

7.4J00303E- 


=02 


76.382443 


=92.865200 


=649.32655 


11.976448 




46.046438 






911.36986 




=397.81182 


=214.30870 


=642.64049 


. 37630623 
7.3139266 




46.046438 






633.76637 

1.0036697 




=637.40698 


=811,66440 


=61.804800 


=.10794966 




46.045438 


- 




=230.22869 

=1.4781238 




=669.46373 


=1096.4742 


266.21960 


=7.3807330 




46.045438 






=846.33640 




. =256.11165 


=406.69090 


522.62789 


=8.7734782E 

=11.800986 


=02 


46.045438 






=994.84164 




164.14776 


89.948200 


567.23478 


7.4643684E=02 








=11.800904 




46.045438 






=910.83630 




539.84917 


208,32270 


398.93070 


.37363278 

=7.3806145 




46.045438 






= 537.65125 

.95640642 




719.53842 


807.70720 


49 ,118100 


=. 10767439 




46.046438 
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POTENTIAL SOLUTION - ST OKI AN WAVE 
RUN NUMBER P5 



WL 


BL 


H 


DEPTH 


Z 


RHO 


vise 


ALPHA 


SI 


S2 


S3 


S4 


S5 


S6 


S7 




400. 


30. 


25. 


200. 


1 

cn 

O 

• 


1.93 


9.66E-06 


.785 


3*0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 





FORMAT 



Vertical 


Vertical 


Lateral 


Lateral 


Force 

AXIAL FORCE 
ETA 


Moment 


Force 


Moment 


=4.8413650 

868.52576 

=.10785793 


2013.8965 


1735.6688 


=3.4335000 


1416.1622 

701.67832 

7.3140010 


1637.0074 


1402.2395 


=852.92800 


2297.7335 

267.78543 

11.975477 


627. 10130 


535.14457 


=1375.5690 


2297.7279 

=267.79148 

11.975448 


=627.11540 


-535.15664 


-1375.5659 


1416.1478 
=701.68207 
7; 3 139255 


=1637.0154 


-1402.2470 


=852.91890 


=4.8590710 

=868.52575 

=.10794965 


=2013.8952 


-1735.6687 


-3.4230000 


=1 41 9. 1 686 
=703.62060 
=7.3807330 


=1621.5329 


=1406.1210 


850.81360 


=2289.9054 
=268.98205 
=1 1 . 800986 


•617. 52880 


=537.53592 


1381.1282 


=2289.8892 

269.00030 

=11.800904 


617. 57130 


537.57235 


1381.1181 


=1419.1257 

703.63184 

=7.3805145 


1621.5587 


1406.1434 


850.78790 


=4.8059270 

868.52572 

=.10767439 


2013.8955 


1735.6686 


-3.4546000 
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MODIFIED MORRISON SOLUTION - STOKIAN WAVE 
RUN NUMBER P6 



WL 


BL 


H 


DEPTH 


z 


RHO 


vise 


ALPHA 


SI 


S2 


S3 


S4 


S5 


S6 


S7 




400 . 


0 

O 


0 

un 

CM 


200. 


»50. 


1.93 


9.66E-06 


1.5708 


3.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 





Vertical 


FORMAT 

Lateral 


Vertical 


Lateral 


Force 


Force 


Moment 


Moment 


RATIO 

ETA 

A ol 00000000 


CELERITY 






=543.05184 


1070.6491 


-4.2000000E-03 


=1. 1000000E=03 


.99634869 

=.10785793 

228.52365 


46.045438 

1140.3940 


=2.7000000E=03 


2,0000000E=04 


=1.5028611 

7.3140010 


46.045438 






848.27305 


914.79429 


=2.9000000E=03 


2. 8000000E=03 


=9.4046414E=02 

11.975477 

1007.8263 


46.045438 

286.66058 


1 . 0000000E=04 


4.6700000E=03 


7.9158946E=02 

11.975448 

915.40059 


46.045438 

=505.51289 


1 .0000000E=03 


4. 1000000E=03 


.37518376 

7.3139255 


46.045438 


- 


539.07899 


=966.65891 


3.7000000E-03 


8.0000000E=04 


1.0036914 

=.10794965 

=231.53307 


46.045438 

=840.77703 


5.0400000E=03 


=1.8000000E=03 


=1.4756470 

=7.3807330 


46.045438 


■ 




=845.33244 


=379.68211 


9. 0000000E=04 


=3. 7800000E=03 


=9.4059494E=02 

=11.800986 

=1004.3534 


46.045438 

251.31371 


4. 0000000 E=04 


=3.3400000E“03 


7.9171370E=02 
=1 1 . 800904 


46.045438 


- 




=914.84950 


809.75845 


=4.0000000E~03 


-2. 1000000E=03 


.37347413 

=7.3805145 

=543.03754 


46.045438 

1070.6527 


-3.7000000E=03 


=1.0000000E=03 


.99637494 

=.10767439 


46.045438 
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POTENTIAL SOLUTION - ST OKI AN WAVE 



RUN NUMBER P6 



WL 


BL 


H 


DEPTH 


z 


RHO 


vise 


SI 


S2 


S5 


S4 


S5 


S6 


ST 


400. 


30. 


25. 


200. 


“50. 


1.93 


9.66E 


3.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 



FORMAT 



Vertical 


Vertical 


Lateral 


Force 

AXIAL FORCE 
ETA 


Moment 


Force 


“4,9322170 

-4.5256967E-03 

“.10785793 


=1.0474700E“02 


2466.9918 


1422.7199 
=3. 6562199E=03 

7.3140010 


“8.7000000E-03 


1993.0335 


2308.4630 

=1.3953268E=03 

11.975477 


“2.9000000E-03 


760.60336 


2308.4575 

I. 3953583E-03 

II. 975448 


2.8000000E-03 


=760.62053 


1422.7054 

3.6562394E-03 

7.3139255 


7.6000000E=03 


“1993.0441 


=4.9500041 

4.5256967E=03 

=.10794965 


9.2906000E“03 


“2466.9918 


=1425.7825 

3.6664833E=03 

“7.3807330 


7.7000000E-03 


=1998.6281 


“2300.4881 
1 .4016504E-03 
=11,800986 


2.6000000E-03 


“764.05036 


=2300.4716 
-1.4017454E-03 
=1 1.800904 


“2.8000000E-03 


764. 10214 


=1425.7393 

“3.6665419E=03 

“7.3805145 


“7.7000000E-03 


1998.6600 


=4.8966323 
=4.5256967E=03 
=. 10767439 


“9.2877000E-03 


2466.9918 



ALPHA 

•06 1.5708 



Lateral 

Moment 

.00000000 

6.3000000E-03 

8.9000000E-03 

9. 0000000E-=03 

5.9000000E=03 

.00000000 

5.5000000E-03 

•9.0000000E-03 

•8.8000000E=03 

•5 . 9000000E-03 

.00000000 
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APPENDIX C. COMPUTER PROGRAMS 



This section contains the IBM 1620 computer programs 
utilized in the calculations for each of the three theories 
discussed. The programs are written in the standard IBM 
Fortran language as adapted for the IBM 1620 computer. 

As an aid to the reader in the event that a close in- 
spection of the mechanics of a program is desired, a list 
and description of the symbols used in the modified Mor- 
rison theory is given below. The explanation does not 
cover every symbol but does cover the most used symbols 
and those which might be difficult to understand. Although 
the Stokian potential flow program is not explained, its 
structure is similar to the modified Morrison program. The 
sine wave potential program simply involved the programming 
of equations derived directly from Cummins (6) report, and 
for this reason the symbols involved are not discussed. 



Program 

Symbol 


Meaning in Theory Symbols 
(Explanation) 


WL 


L 


(wavelength) 


BL 


$ 


(body length) 


H 


H 


(wave height -double amplitude) 


Depth 


h 


(depth of water) 


z 


z 


(depth of the body axis below 
still water) 
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Program 

Symbol 

RHO 

vise 

ALPHA 
SI, S2 . . 

B 

CEL 

SIGMA 

TIME 

DIAM 

D(J) 

X(J) 

PIHOL 

PIHSQ 

ALSIN 

ALCOS 

Y 

T 



Meaning in Theory Symbols 
(Explanation) 

(fluid density) 

V (fluid viscosity) 

c* (angle of attack) 

o S7 Coefficients of a 6th order parabola 

(Sl+S2x+S3x 2 ... S7x°) describing 
the variation of body diameter over 
the length. x=0 is located at mid- 
ships , +x towards the bow 

K (wavenumber « 2 jt/L) 

C (wave celerity) 

cf“ (wave angular frequency) 

T (wave period) 

Part of 6th order parabola describing 
the diameter 

D(x) Diameter of the body at a parti- 
cular axial position 

Axial position along the axis of 
symmetry . Midship section is the 
origin. x(l), x(2) ... x(7) locate 
the seven Tchebycheff stations . 

rH 




sin c* 
cos c* 

ifj = (Kx cos + Ky sine/. + cr t ) 
Note that y * 0 for all cases 

Stepping parameter simulating the 
passage of time. It was stepped from 
0 to 10 inclusive giving eleven sets 
of forces and moments over the wave 
cycle. 
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Program 

Symbol 


Meaning in Theory Symbols 
(Explanation) 


HH (1 ) 


V 

max 


HH (2 ) 


w 

max 


CMCN 


u T w T 

DlxT - ° r “ 5 ^y-( de P endin S on index L) 


CM (<J,L) 


C (vertical and transverse for 
m each x(J) at a particular time) 


RNCON 


l. v J P ( x ) or -Lw]D_ (depending on 
« v index K) 


CD (K, K) 


C D (corresponds to vertical and 
transverse Reynolds number at 
each x(J)) 


ETA 


y\ (displacement of free surface from 
still water) 


VPI 


AF_ (incremental vertical component 
1 z of vertical force) 


VFD 


AF d (incremental, drag component of 
D z vertical force) 


HFI 


AF, (incremental inertial component 
x y of transverse force) 


HFD 


AF d (incremental drag component of 
y transverse force) 


SMVFI 


F_ (inertial component of vertical 
L z force) 


SMVFD 


F D (drag component of vertical force) 
z 


SMVF 


F (total vertical force) 
z 


SMHF 


F (total transverse force) 


SMVM 


M (moment about y axis ) 
z 


SMHM 


My (moment about z axis ) 


RATIO 


(drag component of vertical force) 
total vertical force 



-i8o- 



The hyperbolic sine and cosine terms required in the 
wave equations were computed by statement number 1 in 
which 

Program Meaning in Theory Symbols 

Symbols (Explanation) 



A (1 ) * Kh 


C(l,l) 


- sinh 


Kh 


C (2,1) * cosh Kh 


A (2 ) * 2Kh 


C(l,2) 


» sinh 


2Kh 


C (2,2) * cosh 2Kh 


A(3) « K(h+z) 


C(l,5) 


» sinh 


K(h+s ) 


C (2,3) * cosh K (h-H& ) 


A (4) - 2K (h+z ) 


C (1,4) 


« sinh 


2K (h+z) 


C (2,4) - sinh 2K(h+z) 



The most illustrative way to show how symbols E(l,l), 

E(l,2) .... etc. are used is by rewriting the transverse 
and vertical velocity and acceleration equations with these 
symbols inserted. 

Thus 

v - E (l,l)sin < / / + F(l,l)cos 2 H' + 0(1,1) (see equation 28) 

w = F(2,2)cos^/' + E(2,2)sin 2^ + G(2,2) (see equation 29) 

dv/dt « F(l,2)cos(/' + E(l,2)sin 2 V' + G(l,2) (see equation 30), 
dw/dt = E(2,l)sin^ + F(2,l)cos 2V 7 + G(2,l) (see equation 31) 

By noting the E(M,N) terms are always multiplied by sine 
terms and the F(M,N) terms are always multiplied by cosine 
terms, I believe that the reader can work out statement 
number 4 (VIP(M,N) .... ) to see that the vertical and trans- 
verse velocity and acceleration terms at a particular x(J) 
and a particular time (I) are as follows; 
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VIP (1,1) * v 
VIP (1,2) - dv/dt 
VIP (2,1) - dw/dt 
VTP (2,2) « w 

It should be noted that the above explanation of the 
symbols in the program should not be considered a complete 
analysis of the program. It is offered only as a guide 
to help those interested in following the programming 
procedure involved in the modified Morrison program. 
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R.J. ENGLERT J,R. THUNE THESIS PROJECT 
FORCES AND MOMENTS ON A BODY OF REVOLUTION IN WATER WAVES 
THEORY BASED UPON THE SUPERPOSITION OF INERTIAL AND 
DRAG COMPONENTS TO OBTAIN TOTAL FORCE 

DIMENSION C(4,2),X(7)iVIP(2,2),D(7)iRN(2,2),CD(2,2),CM(7,2) 
D I MENS ION A(4),E(2,2),F(2,2),G(2,2),HH(2), DSQ(7 ) 

10 PAUSE 

ACCEPT , WL 9 BL , H p DE PTH , Z , RHO , V I SC , ALPHA , S 1 ,S2,S3,S4 f S5,S6,S7 
B=6. 2832/WL 
A( I )=B*DEPTH 
A(2)=2.*A( I ) 

A( 3 )=A( 1 )+B*Z 
A(4)=2.*A(3 ) 

DO 1 J-1,4 
DO I J=1 , 2 

1 C(l p J)=,5*(EXP(A( I ))+((-!. )**J)*1./EXP(A( I ))) 

TANH=C( 1 0 1 )/C(1 „2) 

PiH0L=3. l4l6 Vr H/WL 
PSHSQ=P!HOL*PIHOL 
CSQ=C(1J)**2 
ALS 1 N=S 1 N( ALPHA) 

ALCOS® COS (ALPHA) 

IF(C(1,1M0.0E10)8,9,9 
9 CEL«SQR(32.174*TANH/B) 

F(1,1)-0.0 
E( 2 „2 )=0«0 
ETACN^O. 0 
GO TO 6 

8 SS1=PIHSQ*(5.+2.*C(2,2)+2.*C(2,2)**2)/(8.*CSQ**2) 
CEL-SQR(32. 174*TANH/B*(1 .+SS1 )) 

F(1 ,1 )— PiHSQ*CEL*ALSIN/CSQ*(-.5+.75*C(4,2)/CSQ) 

E( 2 „2 )=-. 75*P I HSQ*CEL*C(4, 1 )/(CSQ*CSQ) 

ETACN=-P I HOL*H* , 25*C (1 , 2 )*( 2 . +C ( 2 , 2 ) ) /CSQ**2 
6 S 1GMA«B*CEL 

E( 1,1 )«P 1 HOL*CEL*C( 3 , 2 )*ALS I N/C( 1 t 1 ) 

G(1 ,1 )=P!H0L*CEL*ALS!N*.5*C(4,2)/CSQ 
F(2,2)=»PIH0L*CEL*C(3 ,1 )/C( 1 ,1 ) 

G(2,1 )=0.0 
G(2,2)=0.0 

E(1,2)--SIGMA*F( 1,0*2, 

F(1,2)-S1GMA*E( 1,1) 

G(1 ,2)-0.0 

E( 2 , 1 )=-S 1 GMA*F (2,2) 

F( 2 „ ! )— S 1 GMA*E( 2 , 2 )*2 . 

HH(1 )-E(1 ,1 )+G(1 ,1 ) 

HH( 2 )=~F (2*2) 

X(1 )— ,44195*BL 
X(2)=-.26485*BL 
X(3)==.16195*BL 
X(4)=0.0 
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X(5)-“X(3) 

X(6)=~X(2) 

X(7)=*-X(1) 

TIME^WL/CEL 
DO 112 J-1 ,7 

D 1 AM=S 1 +S2*X ( J ) +S3*X ( J ) **2+S4*X ( J ) **3+S5*X ( J )**4 
D( J )=D3 AM+S6*X ( J )**5+S7*X ( J )**6 
DSQ( J )sD( J )*D( J ) 

DO 12 L=1 „ 2 
CMCNsHH( L)*TIME/D( J ) 

IF(CMCN=7. 5)30,30, 31 

30 CMCN1=2.22°7.2002670E°02*CMCN 

CMCN2-2 . 320 1 049E~02*CMCN**2~3 . 5200945E=03*CMCN**3 
GO TO 1! 2 

31 IF(CMCN-20. )32,32,33 

32 CMCND-8. 0717750=1 .6777438*CMCN+. 14323456*CMCN**2 
CMCN2«=5.4474432E=03*CMCN**3+7.9358277E-05*CMCN**4 
GO TO 112 

33 CMCN W . 36985EHD2+5 . 1 8527 1 5E=02*CMCN 

CMCN2==4. 3765884E=04*CMCN**2+1 . 5350526E-06*CMCN**3 
12 CM( «J » L )=CMCN 1 +CMCN2 

DO 2 1-1,11 

PP— 1 

T=.62832*(PP°1.) 

SMVF 8=0.0 
SMVFD=0.0 
SMVM-0.0 
SMHF=0.0 
SMHM=0 . 0 

ETA-. 5*(H*S I N(T )+ETACN*COS( 2 .*T) ) 

DO 3 J-1, 7 
YsB*X( J )*ALCOS+T 
DO 4 M-1,2 
DO 4 N-1,2 
P— N 

4 V,ip(M,N)-E(M,N)*SIN(P*Y)+F(M,N)*C0S(2.*Y/P)+G(M,N) 

DO 5 K=1,2 

RNCON=SQR( V I P(K,K)**2)*D( J )/Vl SC 
RN(K,K)«L0G(RNC0N)/2,303 
JF(RN(K,K)=>3. 301 )40,40,41 

40 CDCI^* 97 323833”. 6286544 1*RN(K,K)+. 1670641*RN(K,K)**2 
CDC2==1.5636585E=02*RN(K,K)**3“1 .8406664E~03*RN(K,K)**4 
GO TO 50 

41 IF(RN(K,K)-5. 699)42,42, 43 

42 CDC1— 1 . 2369606+3 . 6492456*RN( K ,K)-2. 27761 74*RN( K , K )**2 
CDC2».52382332*RN(K,K)**3“4.0687398E-02*RN(K,K)**4 

GO 7 0 50 

43 CDC 1=3=7 .7171679+2. 75921 1 3*RN(K„ K) 
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CDC2=“ . 34 1 63942*RN ( K , K )**2+1 . 395 87 82 E-02*RN ( K,K)**3 
50 P0P«CDC1+CDC2 
5 CD(K»K)=10.**( POP) 

VF l=CM( J , 2 )*. 7854*DSQ( J )*V I P(2 , 1 )*RH0*BL/7 . 

Vf D«=CD( 2 , 2 )*RHO*. 5*SQR( VI P(2 , 2 )**2)*V 1P(2,2)*D(J )*BL/7 . 
HFNRH0*CM(J,1 )*.7854*DSQ(J)*VIP(1,2)*BL/7.. 

HFD=CD(1 P 1 )*.5*D(J)*SQR(V!P(1 J )**2)*VIP(1 ,1 )*BL/7. 

HFsHF 1+HFD 

SMVF NSMVF 1 +VF 8 

SMVFDsSMVFD+VFD 

SHVMsSMVH+( VFD+VF 8 )*X( J ) 

SMHF^SMHF+HF 
3 SMHM-SMHM+HF*X ( J ) 

SHVFsSMVF I +SMVFD 
RAT 1 OsSMVFD/SMVF 
PRINT, SHVF , SMHF , SMVM , SMHM , RAT 1 0 
2 PRINT, ETA, CEL 
GO TO 10 
END 
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POTENTIAL FLOW SOLUTION FOR STOKIAN WAVE PROFILE 
FORCES AND MOMENTS ON A SUBMERGED BODY OF REVOLUTION 
R J ENGLERT J R THUNE THESIS 



DIMENSION A(4),C(4,2),E!(2),E2(2),E3(2),X(7),F(3) 

DIMENSION V(2),UT(3),Q1(2),Q2(2) 

10 PAUSE 

ACCEPT, WL, BL,H, DEPTH, Z,RHO, VI SC , ALPHA, SI , S2 , S3 , S4, S5 , S6, S7 
B=6.2832/WL 
A( 1 )-B*DEPTH 
A(2)-2.*A(1 ) 

A(3)=A(1)+B*Z 
A(4)«2**A(3 ) 

DO 1 1-1,4 
DO 1 J-1,2 

1 C(l # J)=.5*(EXP(A( I ))+((-1. )**J)*1./EXP(A( I ))) 
TANH=C(1 5 1I)/C(1,2) 

P1H0L-3. 1416*H/WL 
P 1 HSQsP I HOL*P 1 HOL 
ALS 6 N-S 1 N( ALPHA) 

ALCOS^COS( ALPHA) 

IF(C(1,1)-10.E10)8,8,9 
9 CEL=SQR(32.174*TANH/B) 

P3»0.0 
ETACN-0.0 
GO TO 6 

8 SS1-PIHSQ*(5*+2.*C(2,2)+2.*C(2,2)**2)/(8.*C( 1 , 1 )**4) 
CEL-SQR(32. D74*TANH/B*( 1 .+SS1 ) ) 

ETACN—P I HOL*H*. 25*C ( 1 , 2 )* ( 2 . +C ( 2 , 2 ) ) /C ( 1 , 1 )**4 
P3s=*P!HSQ*CEL*,75/C( H , 1 )**4 
6 SIGMA - B*CEL 

Pl-P I HOL*CEL/C( 1,1) 

P2-PIHSQ*CEL*,5/C(! ,1 )**2 
P4-P2 

EJ (1 )»S 1GMA*C(3 p 1 ) 

E2( 1 )=0.0 

E3( 1 )=2.*S 0GMA*C(4» II ) 

El (2)*C(3,2)*SIGMA 
E2( 2 )==S I GMA*2 . 

E3 ( 2 )-=C( 4 , 2 )*S 1 GMA*2 , 
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X( 1 )==.44!95*BL 
X ( 2 )— = * 2.6485*BL 
X(3)-“.16!95*BL 
X(4 )=0 . 0 
X(5)-~X(3) 

X(6)=~X(2) 

X(7)»“X(I) 

CC==1 2 . 5664*RH0 
DO 30 1-1, n 
PP=C 

T=.62832*(PP°K ) 

ETA=,5*(H*S IN(T)+ETACN*C0S(2.*T ) ) 

DO 3 M=»J „3 

3 F(M)=Q.O 
VM=0,0 
HM«0.0 

DO 20 J-1,7 
PSB«B*X(J)*ALCOS+T 
Q? ( 1 )=S ! N( PS ! ) 

Q2(1 )-C0S(2.*PSl ) 

Q1 ( 2)=C0S( PS I ) 

Q2(2 )=S 8 N( 2 . *PS i ) 

DO 4 K-1,2 

4 V(K)=P1*E1(K)*Q1(K)+(P2*E2(K)+P3*E3(K))*Q2(K) 

D I AM-S1 +S2*X( J )+S3*X( J )**2+S4*X( J )**3+S5*X( J )**4 

D=D6AM+S6*X(J)**5+S7*X(J)**6 

UT ( 1 )»-V(2)*ALC0S*D**2/16. 

UT( 2 )— V( 2 )*ALS I N*D**2/8 . 

UT(3) aes V( 1 )*D**2/8. 

DO 40 MM ,3 

40 F(M)=CC*UT(M)*BL/7.+F(M) 

VM=CC*UT(3 )*X( J )*BL/7 . +VM 
20 H^CC*UT(2)*X(J)*BL/7«+HM 
PR 0 NT » F ( 3 ) t VM , F ( 2 ) » HM , F ( 1 ) 

30 PRINT, ETA 
GO TO 10 
END 
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C R J ENGLERT J R THUNE THESIS PROJECT 

C POTENTIAL FLOW SOLUTION FOR SINE WAVE PROFILE 

DIMENSION X(7),D(7),DSQ(7) 

10 PAUSE 

ACCEPT, WL,BL,H, DEPTH, Z,RHO, VI SC, ALPHA, SI ,S2,S3,S4,S5,S6,S7 

B=6. 2832/WL 

ALS I N=S 8 N( ALPHA) 

ALCOS=COS( ALPHA) 

CEL-SQR(5.12062*WL) 

X(1 )=».44195*BL 
X(2)s=.26485*BL 
X(3)=“. 16195*BL 
X(4)-0.0 
X(5)*°X(3) 

X(6)==X(2 ) 

X(7)==X(1) 

T 8ME=WL/CEL 

PTC0N®B*RH0*25 . 26632*H*EXP( B*Z ) 

DO 12 J-1,7 

D ! AM=S 1 +S2*X ( J ) +S3*X ( J )**2+S4*X ( J ) **3+S5*X ( J ) **4 
D ( J )*D 8 AM+S6*X ( J )**5+S7*X ( J )**6 
12 DSQ(J)-D(J)**2 

DO 2 1-1,11 

PP-I 

T».62832*(PP~1.) 

SMPVF=0,0 

SMPVM-0,0 

SMPHF=0,0 

SMPHM=0.0 

ETA^.5*H*S8N(T) 

DO 3 J*1 » 7 
Y=B*X(J)*ALCOS+T 

$MPHF=SMPHF+PTC0N*ALSIN*DSQ(J)*C0S(Y)*BL/7. 

SMPHM«SMPHM+PTC0N*ALSIN*DSQ(J)*X(J)*C0S(Y)*BL/7. 

SMPVF=SMPVF+PTC0N*DSQ(J)*SIN(Y)*BL/7. 

3 SMPVM^SMPVM+PTCON*DSQ( J )*X( J)*SIN(Y )*BL/7 . 

PRINT, SMPVF , SMPHF , SMPVM , SMPHM 
2 PRINT, ETA 

0F( SENSE SWITCH 1)20,21 

20 PRINT, Y,PTCON,T 

21 PRINT, CEL 
GO TO 10 
END 
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APPENDIX Do DESIGN OF TORSION GAGE 



A. General 

A particular array of strain gages mounted on a thin 
walled cylindrical tube was found to be an effecient and 
simple design for the torsion gage (14). The arrangement 
of the strain gages is illustrated in Fig., 23 . 



Configuration of gages 




Figure 23 




Circuit Diagram 
Figure 2k 
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The configuration shown utilizes the fact that the 
principal strain axes of a shaft in torsion are 45° from 
the circumferential reference line. Note also that the 
magnitude of the principal compressive strain is equal to 
the magnitude of the principal tensile strain. By using 
the principle of the Wheatstone bridge which states that 
the output of opposite arms of a strain gage bridge are 
additive and adjacent arms subtractive, the gages were con- 
nected as indicated in Fig. 24 for maximum sensitivity. 

The characteristics of the torsion gage described above 
are as follows: 

1. Output is proportional to the 

applied torque. 

2. Output is independent of bending. 

3. Output is independent of axial 

thrust . 

4. As long as temperature changes take 

place uniformly throughout the 
tube, the output is independent of 
temperature changes . 

5. The total output is four times the 

output of a single gage. 



B. Calculation of Gage S ize and Sensitivity 

From the basic relationship of solid mechanics between 
stress and torsional moment applied to a shaft or cylindrical 
tube: 
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r p 



( 46 ) 



in which T 



I_ 



shear stress at the outer surface 
of the tube or shaft „ 

outside diameter. 

moment in the horizontal plane 
(torque applied) 



polar moment of inertia of the section, 
I ^ 



JL d 
52 



jl et- 

52 



for solid shaft 
,4 



i) 



Substituting the basic stress strain relationship: 



7 “ 



G 



in which if - shear strain at the surface 

G ** modulus of elasticity in shear 
(modulus of rigidity) 



Equation 46 then beomess 



ifG « M. 



y i 



o/2 



or 



M. 



I G 

-t 



P 



2 if 



(47) 



Prom an investigation of Mohrs circle for the torsional 



case 
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In which: = Principal 

tensile strain 



* Principal com- 
pressive strain 



it can be seen that <f 




€ 



2 



substituting X Into 47? 



- <€i - e 2 ) 

a 

o 



Also from the Mohrs circle: 




(48) 



then 



£ p » -£ 1 (Maximum compressive strain ip equal 
in magnitude but opposite in sign 
to the maximum tensile strain) 



M 



d o 



(2 q) 



(49) 



The above equation is applicable to two strain gages 
located on one side of the gage. By adding the output of 
the other two gages the torque equation becomes: 




Substituting 





G 



E 

2(1 + H) 



(50) 



into equation (50) results in 



2(1 + n) 



I 

d„ 



(4e x ) 



(51) 
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Because of difficulties encountered in bonding strain 
gages to small radius objects, a 1-inch outside diameter 



tubing was considered to be the minimum practical size. 

The Sanborn recorder instruction manual indicated 

the maximum sensitivity of a 4 active arm bridge to be 

about 10 micro inches/inch/ centimeter deflection of the 

Sanborn stylus. Therefore, by substituting this value 

for €, in equation (51) and the value of I corresponding 
— P 

to a one inch outside diameter, ,055 inch wall thickness 
aluminum tubing, (smallest commercial size available), the 
torque necessary to produce one centimeter of stylus de- 
flection can be calculated. Substituting in equation (50) 
the following values; 



E ■ 10 x 10° for aluminum 
£ 1 - 10 X 10~ 6 



d 



I 



P 



o 



1.0 

.0247 in 21 



then 




M = 5.71 inch - lbs/cm 



or 



M 



y 



.509 ft. -lbs/cm 
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This figure agrees very well with the actual calibration 
trace of the torsion gage which yielded .52 ft. -lbs/cm. 

The strain gages used for this instrument were Tatnall 
foil gages, type 121-R5A. 

Eastman 910 cement was used to bond the gages to the 
tubing. The area in which the gages were located was 
thoroughly painted with TMS-GW-1 type waterproofing in 
order to waterproof the gages and to hold the fragile 
wiring in place. The gage itself proved to be very ef- 
fective in measuring the moment in the horizontal plane 
about the vertical axis of the body. 

Although the torsion gage lowered the natural frequency 
of the overall gage system slightly in bending, the natural 
frequency of the system in torsion was about 50 cps, thus 
being well out of the range of the forcing frequencies used. 
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